已有未经几何纠正卫星影像的四个角的WGS空间直角坐标,如何通过gdal完成几何精纠正?

5 views
Skip to first unread message

277480729

unread,
Mar 3, 2009, 8:50:43 AM3/3/09
to geosings
   
 
我有一张未经几何纠正的原始卫星影像,现在想给它做几何纠正。通过一系列数学计算,已计算出影像的四个角坐标(其实可以算出任一像元对应的地面坐标---WGS-84坐标系统)。
用gdalwarp可以直接把这幅原始卫星影像进行几何粗纠正(把图像从无投影变成wgs-84下的投影坐标),请问我是不是可以通过更改这幅影像的gcp信息,把高精
度计算出的图像四个角的坐标直接当作gcp来纠正图像(暂时把这个过程当作几何精纠正)?
接下来我该怎么做才能完成这个几何精纠正呢?我是不是可以更改geotif图像的tag文件?进行自定义gcp,以通过gdalwarp完成纠正。
谢谢先!
这是一幅未经任何几何纠正的原始卫星影像的gdalinfo 信息,它没有投影,但是自带gcp有投影信息,我能否改动它自带的gcp?
C:\Dwocuments and Settings\dell>gdalinfo f:/scene01/imagery.tif
Driver: GTiff/GeoTIFF
Size is 24000, 24000
Coordinate System is `'
GCP Projection = GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.2572235630016,AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]]
GCP[  0]: Id=1, Info=
          (1,1) -> (102.808469142837,25.3221393577761,0)
GCP[  1]: Id=2, Info=
          (24000,1) -> (103.38543742346,25.1997490622572,0)
GCP[  2]: Id=3, Info=
          (24000,24000) -> (103.247943888031,24.6709097848961,0)
GCP[  3]: Id=4, Info=
          (1,24000) -> (102.673419826765,24.7930519134803,0)
Metadata:
  AREA_OR_POINT=Point
  TIFFTAG_IMAGEDESCRIPTION=PAN
  TIFFTAG_DATETIME=2004:04:19 11:50:37
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,24000.0)
Upper Right (24000.0,    0.0)
Lower Right (24000.0,24000.0)
Center      (12000.0,12000.0)
Band 1 Block=24000x128 Type=Byte, ColorInterp=Gray
  Overviews: 6000x6000, 3000x3000, 1500x1500, 750x750, 375x375, 188x188, 94x94, 47x47
 
 
     

linux23...@gmail.com

unread,
Mar 8, 2009, 10:59:47 AM3/8/09
to gdal+python+GIS+geosings论坛

呵呵,不好意思,没有非常认真得玩过GCPS,不敢乱说……

玩过的同学可以说说滴嘛

Reply all
Reply to author
Forward
0 new messages