1

Topic: a problem of orthorectifition and resampleing

In  order  to  orthorectifite  a  ASTER  raw  image  with  a  shadea  DEM,  I  selected  15  tie  points  between  them  and  then  optimized  GCPS.  the  optimized  result  shows  that  the  mean  of  their  misregistration  are  close  to  zero.  All  the  above  processes  are  no  problem.
>>>  when  the  step  of  orthorectification/resampling,the  programe  popup  a  dialog  box  titled  "system  error".the  content  of  the  dialog  box  is  as  follows
>>>      "  SAT_ORTHORECTIFICATION--Array  subscript  for  SARR  must  have  same  size  as  source  expression".
>>>
>>>  I  am  eager  to  know  what  is  the  reason  of  this  error,and  how  to  resolve  it  ??
>>>  I  hope  someone  may  help  me!
>>>
>>>  the  dem  I  used  do  not  fully  coverage  the  raw  image.  the  overlap  between  the  two  images  are  most.
>>>  the  resolution  of  the  dem  is  30m.the  dem  is  ASTER  GDEM.
>>>  to  define  the  map  grid,  I  used  "from  the  raw  image"

>>> the optimized GCPS file is as follows:

;File generated the : Sun Aug 29 11:05:14 2010
;GCPS generated from :
;- Tied point file : D:\ѧϰ\ASTER\yushu\yushu_post\shaded_image1763.pts
;- GeoReferenced image : D:\ѧϰ\ASTER\yushu\yushu_dem\shaded-relif
;- DEM : D:\ѧϰ\ASTER\yushu\yushu_dem\ASTGTM_N33E096\ASTGTM_N33E096_dem.tif
;Values : Easting, Northing, Altitude, X, Y, Weight, Opti, dE, dN, dA
;Note : Opti, dE, dN, and dA, are only used in GCPS Optimization
96.7155733960     33.4282174676    3779.6907    2818    3020   0.36574   1     0.0000     0.0000     0.0000
96.7505207920     33.3900270203    3690.4805    3110    3249   0.55679   1     0.0000     0.0000     0.0000
96.7059880087     33.4664363556    3750.2556    2689    2726   0.38150   1     0.0000     0.0000     0.0000
96.6388924214     33.4657709787    3816.7316    2302    2823   0.33411   1     0.0000     0.0000     0.0000
96.6081568456     33.3843767205    4528.4871    2264    3436   0.60826   1     0.0000     0.0000     0.0000
96.5705455785     33.4911762643    4035.9469    1833    2704   0.39087   1     0.0000     0.0000     0.0000
96.4920590754     33.6089986546    3894.0418    1176    1922   0.65289   1     0.0000     0.0000     0.0000
96.4450994470     33.6361118811    4855.1261     868    1783   0.46216   1     0.0000     0.0000     0.0000
96.4504679934     33.6771976855    3827.4429     807    1485   0.52156   1     0.0000     0.0000     0.0000
96.5393179584     33.6745738287    4798.0535    1348    1422   0.51217   1     0.0000     0.0000     0.0000
96.6633594447     33.6411341096    3894.2910    2160    1520   0.51807   1     0.0000     0.0000     0.0000
96.7246479414     33.5820308705    4006.8265    2625    1873   0.61792   1     0.0000     0.0000     0.0000
96.5001474073     33.6983131962    4193.1315    1059    1286   0.52653   1     0.0000     0.0000     0.0000
96.4442981657     33.7699986303    3931.4453     603     815   0.53378   1     0.0000     0.0000     0.0000
96.5115787257     33.8127157263    4018.6572     947     421   0.74615   1     0.0000     0.0000     0.0000

2

Re: a problem of orthorectifition and resampleing

Hello,

The error you're mentioning seems to come from a pb with reading the GCP file. The file you posted is the initial gcp file. Would you mind posting the optimized GCP file?

Also, you may have a pb if the DEM doesn't cover the whole raw image.

Sebastien

3

Re: a problem of orthorectifition and resampleing

hello,
as you told me,the GCP file I used is wrong.I retried it with a optimized GCP file which is as follows:
the program run correctly.and I completed the orthorectifition and resampling
Thank you very much!
optimized GCP file

;File generated the : Sun Aug 29 11:33:41 2010
;GCPS optimization report:
;GCP / Easting misregistration (meter) / Northing Misregistration (meter)
;     1    -0.000243    -0.000046
;     2    -0.000042    -0.000022
;     3    -0.000289    -0.000375
;     4    -0.000320     0.000009
;     5     0.000270     0.000311
;     6    -0.000412    -0.000103
;     7     0.000212    -0.000507
;     8     0.000549    -0.000099
;     9    -0.000007    -0.000157
;    10    -0.000233     0.000446
;    11     0.000064     0.000193
;    12     0.000040    -0.000111
;    13    -0.000556     0.000192
;    14     0.000072     0.000012
;    15     0.000075     0.000015
;Easting mis-registration in meters (Average / Standard Deviation):     -0.0000256     0.0002701
;Northing mis-registration in meters (Average / Standard Deviation):     -0.0000129     0.0002413
;Norm mis-registration in meters (Average / Standard Deviation):      0.0000287     0.0003622
;
;Files:
;- Image: D:\ѧϰ\ASTER\yushu\yushu_post\AST_L1A_00304302010042653_20100526224204_1763.hdf
;- Reference Image: D:\ѧϰ\ASTER\yushu\yushu_dem\shaded-relif
;- DEM: D:\ѧϰ\ASTER\yushu\yushu_dem\ASTGTM_N33E096\ASTGTM_N33E096_dem.tif
;- Ancillary file: D:\ѧϰ\ASTER\yushu\yushu_post\AST_L1A_00304302010042653_20100526224204_1763.anc
;- GCPS/Tied Points/ICP file: D:\ѧϰ\ASTER\yushu\yushu_post\shaded_image1763.pts
;
;Resampling: Sinc
;Correlator: Statistic
;- Window Size X: 128
;- Window Size Y: 128
;- Search Range X: 20
;- Search Range Y: 20
;Easting ouput pixel size (meter):   0.00
;Northing ouput pixel size (meter):   0.00
;
;Convergence quality report:
;Loop number, Avg x, Avg y, Norm Avg xy, Stdev x, Stdev y, Norm Stdev xy
;     1    -0.000416     0.000125     0.000434     0.000461     0.000524     0.000698
;     2    -0.000082     0.000147     0.000168     0.000271     0.000242     0.000363
;     3    -0.000031     0.000032     0.000045     0.000265     0.000219     0.000344
;     4    -0.000010    -0.000003     0.000010     0.000286     0.000218     0.000359
;     5    -0.000026    -0.000013     0.000029     0.000270     0.000241     0.000362
;
;Original GCPS:
;GCP number, Longitude (degree), Latitude (degree), Altitude (meter), X (pix), Y (pix), SNR, Opti, dE (meter), dN (meter), dA (meter)
;     1     96.7134722211     33.4281944456    3751.0000    2818    3020   1.00000   1     0.0000     0.0000     0.0000
;     2     96.7493055545     33.3898611122    3687.0000    3110    3249   1.00000   1     0.0000     0.0000     0.0000
;     3     96.7037499989     33.4681944456    3735.0000    2689    2726   1.00000   1     0.0000     0.0000     0.0000
;     4     96.6365277767     33.4651388900    3818.0000    2302    2823   1.00000   1     0.0000     0.0000     0.0000
;     5     96.6084722211     33.3818055567    4534.0001    2264    3436   1.00000   1     0.0000     0.0000     0.0000
;     6     96.5679166656     33.4906944456    3937.9999    1833    2704   1.00000   1     0.0000     0.0000     0.0000
;     7     96.4923611100     33.6104166678    3854.0000    1176    1922   1.00000   1     0.0000     0.0000     0.0000
;     8     96.4476388878     33.6356944456    4988.0000     868    1783   1.00000   1     0.0000     0.0000     0.0000
;     9     96.4498611100     33.6770833344    3833.0000     807    1485   1.00000   1     0.0000     0.0000     0.0000
;    10     96.5379166656     33.6720833344    4926.9999    1348    1422   1.00000   1     0.0000     0.0000     0.0000
;    11     96.6637499989     33.6406944456    3883.0000    2160    1520   1.00000   1     0.0000     0.0000     0.0000
;    12     96.7245833322     33.5831944456    3975.0000    2625    1873   1.00000   1     0.0000     0.0000     0.0000
;    13     96.4973611100     33.6970833344    4040.0000    1059    1286   1.00000   1     0.0000     0.0000     0.0000
;    14     96.4445833322     33.7695833344    3928.0000     603     815   1.00000   1     0.0000     0.0000     0.0000
;    15     96.5120833322     33.8129166678    4014.0000     947     421   1.00000   1     0.0000     0.0000     0.0000
;
;Optimized GCPS:
96.7155733960     33.4282174676    3779.6907    2818    3020   0.36574   1     0.0000     0.0000     0.0000
96.7505207920     33.3900270203    3690.4805    3110    3249   0.55679   1     0.0000     0.0000     0.0000
96.7059880087     33.4664363556    3750.2556    2689    2726   0.38150   1     0.0000     0.0000     0.0000
96.6388924214     33.4657709787    3816.7316    2302    2823   0.33411   1     0.0000     0.0000     0.0000
96.6081568456     33.3843767205    4528.4871    2264    3436   0.60826   1     0.0000     0.0000     0.0000
96.5705455785     33.4911762643    4035.9469    1833    2704   0.39087   1     0.0000     0.0000     0.0000
96.4920590754     33.6089986546    3894.0418    1176    1922   0.65289   1     0.0000     0.0000     0.0000
96.4450994470     33.6361118811    4855.1261     868    1783   0.46216   1     0.0000     0.0000     0.0000
96.4504679934     33.6771976855    3827.4429     807    1485   0.52156   1     0.0000     0.0000     0.0000
96.5393179584     33.6745738287    4798.0535    1348    1422   0.51217   1     0.0000     0.0000     0.0000
96.6633594447     33.6411341096    3894.2910    2160    1520   0.51807   1     0.0000     0.0000     0.0000
96.7246479414     33.5820308705    4006.8265    2625    1873   0.61792   1     0.0000     0.0000     0.0000
96.5001474073     33.6983131962    4193.1315    1059    1286   0.52653   1     0.0000     0.0000     0.0000
96.4442981657     33.7699986303    3931.4453     603     815   0.53378   1     0.0000     0.0000     0.0000
96.5115787257     33.8127157263    4018.6572     947     421   0.74615   1     0.0000     0.0000     0.0000

4

Re: a problem of orthorectifition and resampleing

Hello,

There seems to be a pb because the pixel size of your reference image (the shaded DEM) is zero! look at the field:
;Easting ouput pixel size (meter):   0.00
;Northing ouput pixel size (meter):   0.00

Hence the very small mis-registration, which is representative of the reality.

Sebastien

5

Re: a problem of orthorectifition and resampleing

although this error ,i still got an orthorectified images. is it wrong?will it bring more error to the next steps?
and what's the method to resolve it?
many thanks!

6

Re: a problem of orthorectifition and resampleing

Yes, this might be wrong. You have no way of assessing the quality of your GCP under these circumstances and the
GCP are therefore not optimized. Except for the numerical precision errors, your optimized GCP will be identical to your initial GCP.

Sebastien

7

Re: a problem of orthorectifition and resampleing

Thank you very much for your patience.
I  don't know if I understand correctly.to solve the pb,is it necessary to use the DEM which cover the whole image?

8

Re: a problem of orthorectifition and resampleing

Put simply, the answer is yes. There are ways around it if the DEM doesn't cover the whole area, but it's more complicated. For simplicity and to make sure things work well, I would suggest using a DEM that entirely encompasses the footprint of your images.

Sebastien

9

Re: a problem of orthorectifition and resampleing

Thank you very much for all your help!I wil take a futher experiment with a whole DEM.
thank you one more time!

10

Re: a problem of orthorectifition and resampleing

Dear Sebastien,
I have retried a experiment,and this time I used a DEM which cover the raw image fully enough.
but it still has the same problem.
the optimized GCPS file is:
;File generated the : Wed Sep 01 15:00:19 2010
;GCPS optimization report:
;GCP / Easting misregistration (meter) / Northing Misregistration (meter)
;     1     0.000137     0.000244
;     2    -0.000124    -0.000225
;     3     0.000489     0.000027
;     4    -0.000366    -0.000085
;     5    -0.000287     0.000096
;     6    -0.000157    -0.000332
;     7     0.000356    -0.000193
;     8    -0.000120     0.000016
;     9    -0.000050    -0.000096
;    10    -0.000134     0.000343
;    11    -0.000210     0.000229
;    12    -0.000102    -0.000106
;    13     0.000153    -0.000040
;    14     0.000652    -0.000900
;    15     0.000058    -0.000018
;    16     0.000117     0.000063
;Easting mis-registration in meters (Average / Standard Deviation):      0.0000147     0.0002455
;Northing mis-registration in meters (Average / Standard Deviation):     -0.0000257     0.0002386
;Norm mis-registration in meters (Average / Standard Deviation):      0.0000296     0.0003423
;
;Files:
;- Image: D:\ѧϰ\ASTER\yushu\yushu_post\AST_L1A_00304302010042653_20100526224204_1763.hdf
;- Reference Image: D:\ѧϰ\ASTER\yushu\yushu_dem\shaded-relif_2
;- DEM: D:\ѧϰ\ASTER\yushu\yushu_dem\yushu.mos
;- Ancillary file: D:\ѧϰ\ASTER\yushu\yushu_post\AST_L1A_00304302010042653_20100526224204_1763.anc
;- GCPS/Tied Points/ICP file: D:\ѧϰ\ASTER\yushu\yushu_dem\shaded_image1763_2.pts
;
;Resampling: Sinc
;Correlator: Statistic
;- Window Size X: 128
;- Window Size Y: 128
;- Search Range X: 30
;- Search Range Y: 30
;Easting ouput pixel size (meter):   0.00
;Northing ouput pixel size (meter):   0.00
;
;Convergence quality report:
;Loop number, Avg x, Avg y, Norm Avg xy, Stdev x, Stdev y, Norm Stdev xy
;     1    -0.000015    -0.000186     0.000187     0.000567     0.000301     0.000642
;     2    -0.000131    -0.000042     0.000138     0.000259     0.000276     0.000378
;     3     0.000010    -0.000026     0.000028     0.000246     0.000237     0.000341
;     4     0.000000    -0.000018     0.000018     0.000255     0.000247     0.000355
;     5     0.000015    -0.000026     0.000030     0.000246     0.000239     0.000342
;
;Original GCPS:
;GCP number, Longitude (degree), Latitude (degree), Altitude (meter), X (pix), Y (pix), SNR, Opti, dE (meter), dN (meter), dA (meter)
;     1     96.7493055556     33.3884722222    3695.0000    3105    3253   1.00000   1     0.0000     0.0000     0.0000
;     2     96.7040277778     33.4684722222    3729.0000    2687    2721   1.00000   1     0.0000     0.0000     0.0000
;     3     96.6906944444     33.5651388889    3762.0000    2454    2047   1.00000   1     0.0000     0.0000     0.0000
;     4     96.7551388889     33.5079166667    4080.0000    2933    2385   1.00000   1     0.0000     0.0000     0.0000
;     5     96.8104166667     33.4606944444    3974.0000    3345    2669   1.00000   1     0.0000     0.0000     0.0000
;     6     96.5412500000     33.6190277778    3814.0000    1458    1807   1.00000   1     0.0000     0.0000     0.0000
;     7     96.4645833333     33.6190277778    3758.0000     977    1883   1.00000   1     0.0000     0.0000     0.0000
;     8     96.4056944444     33.5904166667    4048.0000     684    2167   1.00000   1     0.0000     0.0000     0.0000
;     9     96.4962500000     33.5018055556    3941.0000    1375    2717   1.00000   1     0.0000     0.0000     0.0000
;    10     96.5148611111     33.6829166667    4142.0000    1190    1364   1.00000   1     0.0000     0.0000     0.0000
;    11     96.6404166667     33.7237500000    3985.0000    1883     942   1.00000   1     0.0000     0.0000     0.0000
;    12     96.7329166667     33.7709722222    4075.0000    2358     491   1.00000   1     0.0000     0.0000     0.0000
;    13     96.7570833333     33.6559722222    4037.0000    2690    1312   1.00000   1     0.0000     0.0000     0.0000
;    14     96.7773611111     33.5801388889    4328.0000    2942    1827   1.00000   1     0.0000     0.0000     0.0000
;    15     96.5120833333     33.8129166667    4014.0000     949     418   1.00000   1     0.0000     0.0000     0.0000
;    16     96.4456944444     33.7706944444    3957.0000     604     816   1.00000   1     0.0000     0.0000     0.0000
;
;Optimized GCPS:
96.7489716261     33.3890856141    3688.5072    3105    3253   0.57133   1     0.0000     0.0000     0.0000
96.7049669440     33.4669360381    3753.5044    2687    2721   0.37997   1     0.0000     0.0000     0.0000
96.6888257485     33.5650906997    3814.3581    2454    2047   0.51261   1     0.0000     0.0000     0.0000
96.7573455354     33.5071743865    4143.1276    2933    2385   0.33517   1     0.0000     0.0000     0.0000
96.8127294257     33.4608591644    3994.7928    3345    2669   0.52043   1     0.0000     0.0000     0.0000
96.5418314812     33.6171083413    3822.5573    1458    1807   0.47585   1     0.0000     0.0000     0.0000
96.4622066273     33.6175926218    3940.3882     977    1883   0.45361   1     0.0000     0.0000     0.0000
96.4054110504     33.5901088578    4050.1611     684    2167   0.57521   1     0.0000     0.0000     0.0000
96.4956151699     33.5007104201    3958.8370    1375    2717   0.57023   1     0.0000     0.0000     0.0000
96.5157300304     33.6847224622    4129.2684    1190    1364   0.49936   1     0.0000     0.0000     0.0000
96.6420104284     33.7249331396    4004.0932    1883     942   0.59293   1     0.0000     0.0000     0.0000
96.7345132122     33.7707098002    4080.3284    2358     491   0.54786   1     0.0000     0.0000     0.0000
96.7573181735     33.6559081215    4038.5099    2690    1312   0.48745   1     0.0000     0.0000     0.0000
96.7749310844     33.5756559315    4503.6835    2942    1827   0.28518   1     0.0000     0.0000     0.0000
96.5119765610     33.8130538386    4012.2125     949     418   0.74263   1     0.0000     0.0000     0.0000
96.4448680958     33.7709363519    3946.0909     604     816   0.53667   1     0.0000     0.0000     0.0000

the first step, I selected 16 tie points between the raw image and the shaded DEM manually.
here,the shaded DEM is not in UTM projection,but the geographic lat/long projection the same as ASTER GDEM.
the second step  optimized the tie points with the anc file, the reference image(shaded DEM) and the DEM file which is used to generate the shaded DEM.the correlation method is statistic


dxy

11

Re: a problem of orthorectifition and resampleing

Hi,

The pb has nothing to do with you. The pb comes from the map information of your DEM, which is weird. Make sure your DEM is in UTM projection with equal resolution in Easting and Northing.

Sebastien

12

Re: a problem of orthorectifition and resampleing

Thank you!
I have made DEM in UTM projection.The result shows correct.
Thanks a lot.

dxy