Le Thursday 25 June 2009 22:04:58 Gong, Shawn (Contractor), vous avez écrit : > hi list, > > RSAT2's NITF product comes in "imagery_##.ntf" and "product.xml" > gdalinfo output for these two are listed below. > > Georeferencing is off between the two files (their corners don't match > using OpenEV). > Q1: Is it true that imagery_##.ntf uses affine transform (see > GeoTransform tag), and product.xml uses 143 GCPs?
For what I know, NITF files can contain only 4 GCPs and the NITF driver tries to compute a GeoTransform from those 4 ones. > Q2: I see that NITF RPC georeferencing has been implemented, > http://trac.osgeo.org/gdal/wiki/rfc22_rpc How do I use the RPC in > RADARSAT2's NITF file? Will the RPC provide better georeferencing? You can use the RPC info by using the '-rpc' option of gdalwarp. > > > > thanks, > Shawn > > > ========================================= > imagery_##.ntf > Driver: NITF/National Imagery Transmission Format > Files: imagery_HH.ntf > Size is 8792, 10070 > Coordinate System is: > GEOGCS["WGS 84", > DATUM["WGS_1984", > SPHEROID["WGS 84",6378137,298.257223563, > AUTHORITY["EPSG","7030"]], > TOWGS84[0,0,0,0,0,0,0], > AUTHORITY["EPSG","6326"]], > PRIMEM["Greenwich",0, > AUTHORITY["EPSG","8901"]], > UNIT["degree",0.0174532925199433, > AUTHORITY["EPSG","9108"]], > AUTHORITY["EPSG","4326"]] > GeoTransform = > 50.03718817626492, 5.656037108659124e-005, 1.153154346122985e-005 > 13.76523037807838, 1.126467725830654e-005, -5.535416735635244e-005 > Metadata: > NITF_FHDR=NITF02.10 > NITF_CLEVEL=06 > NITF_STYPE=BF01 > NITF_OSTAID=RSI > NITF_FDT=20081019203042 > NITF_FTITLE=RS2_RECTIFIED:2008-10-19T14:51:25.306614Z > ... > ... > RPC Metadata: > LINE_OFF=5034 > LINE_SCALE=5036 > SAMP_OFF=4388 > SAMP_SCALE=4538 > LONG_OFF=50.3435 > LONG_SCALE=0.3068 > LAT_OFF=13.5359 > LAT_SCALE=0.3283 > HEIGHT_OFF=0 > HEIGHT_SCALE=700 > LINE_NUM_COEFF=0.0006728793 0.2105967 -1.130644 -1.223845e-006 > -0.0003698973 -1.002738e-007 -3.118422e-007 -0.0006155271 -0.000504439 0 > -1.174964e-008 2.143919e-007 -4.856652e-006 7.659651e-008 -6.164503e-006 > -8.210502000000001e-007 -4.267789e-007 -4.321238e-008 0 -2.272615e-008 > LINE_DEN_COEFF=1 -0.0004301605 0.0002676128 -3.302012e-007 > -3.970334e-006 0 1.027667e-008 -6.566892e-006 2.541771e-006 > 3.784513e-007 0 0 -1.749914e-006 0 3.37938e-007 3.157953e-006 0 0 0 0 > SAMP_NUM_COEFF=0.0003633695 1.147341 0.2557103 -0.03114342 > -0.0006508754 0.0008477953 0.0001658026 -5.388591e-005 -0.0002000257 > 2.006948e-005 -1.95455e-005 -2.907722e-006 -9.035467999999999e-006 > -1.276494e-006 1.372415e-006 -1.699765e-006 -5.031646e-007 > -6.214799e-005 -1.609981e-006 3.722058e-008 > SAMP_DEN_COEFF=1 -0.0001622142 0.0007102185 -0.001001374 > -8.652518e-007 4.728022e-005 1.732678e-005 -2.196369e-006 8.936224e-006 > -1.514486e-006 -2.812931e-006 0 1.072109e-007 1.4604e-007 -1.900576e-008 > -2.564915e-007 0 -6.090468e-006 -3.221172e-007 0 > MIN_LONG=50.1901 > MAX_LONG=50.4969 > MIN_LAT=13.37175 > MAX_LAT=13.70005 > Corner Coordinates: > Upper Left ( 50.0371882, 13.7652304) ( 50d 2'13.88"E, 13d45'54.83"N) > Lower Left ( 50.1533108, 13.2078139) ( 50d 9'11.92"E, 13d12'28.13"N) > Upper Right ( 50.5344670, 13.8642694) ( 50d32'4.08"E, 13d51'51.37"N) > Lower Right ( 50.6505896, 13.3068530) ( 50d39'2.12"E, 13d18'24.67"N) > Center ( 50.3438889, 13.5360417) ( 50d20'38.00"E, 13d32'9.75"N) > Band 1 Block=1024x1024 Type=UInt16, ColorInterp=Gray > > > ========================================= > product.xml > Driver: RS2/RadarSat 2 XML Product > Files: product.xml > lutBeta.xml > lutSigma.xml > lutGamma.xml > imagery_HH.ntf > imagery_HV.ntf > Size is 8792, 10070 > Coordinate System is `' > GCP Projection = GEOGCS[,DATUM[,SPHEROID["WGS > 1984",6378137,298.2572235629972]],PRIMEM["Greenwich",0],UNIT["degree",0. > 0174532925199433]] > GCP[ 0]: Id=1, Info= > (0,0) -> (50.0368241447173,13.7652233850433,0) > GCP[ 1]: Id=2, Info= > (879.099976,0) -> (50.0866031135979,13.775165475226,0) > GCP[ 2]: Id=3, Info= > (1758.19995,0) -> (50.1363862930026,13.7850974953681,0) > GCP[ 3]: Id=4, Info= > (2637.30005,0) -> (50.1861736889333,13.7950194384848,0) > GCP[ 4]: Id=5, Info= > (3516.3999,0) -> (50.2359652797982,13.8049312921043,0) > GCP[ 5]: Id=6, Info= > (4395.5,0) -> (50.2857610923955,13.8148330533982,0) > ... > ... > GCP[139]: Id=140, Info= > (6153.7002,10069) -> (50.5012171146149,13.2772541435474,0) > GCP[140]: Id=141, Info= > (7032.7998,10069) -> (50.5509029539819,13.2871595011868,0) > GCP[141]: Id=142, Info= > (7911.8999,10069) -> (50.6005928456264,13.297055157266,0) > GCP[142]: Id=143, Info= > (8791,10069) -> (50.6502867615454,13.3069410980946,0) > Metadata: > PRODUCT_TYPE=SGF > PIXEL_SPACING=6.25000000e+00 > LINE_SPACING=6.25000000e+00 > BETA_NOUGHT_LUT=lutBeta.xml > SIGMA_NOUGHT_LUT=lutSigma.xml > GAMMA_LUT=lutGamma.xml > SATELLITE_IDENTIFIER=RADARSAT-2 > SENSOR_IDENTIFIER=SAR > BEAM_MODE=F1 > ACQUISITION_START_TIME=2008-10-19T14:51:24.966569Z > NEAR_RANGE_INCIDENCE_ANGLE=3.68685226e+01 > FAR_RANGE_INCIDENCE_ANGLE=3.98800888e+01 > SLANT_RANGE_NEAR_EDGE=9.631139343954010e+05 > Subdatasets: > SUBDATASET_3_NAME=RADARSAT_2_CALIB:BETA0:product.xml > SUBDATASET_3_DESC=Beta Nought calibrated > SUBDATASET_2_NAME=RADARSAT_2_CALIB:SIGMA0:product.xml > SUBDATASET_2_DESC=Sigma Nought calibrated > SUBDATASET_4_NAME=RADARSAT_2_CALIB:GAMMA:product.xml > SUBDATASET_4_DESC=Gamma calibrated > SUBDATASET_1_NAME=RADARSAT_2_CALIB:UNCALIB:product.xml > SUBDATASET_1_DESC=Uncalibrated digital numbers > Corner Coordinates: > Upper Left ( 0.0, 0.0) > Lower Left ( 0.0,10070.0) > Upper Right ( 8792.0, 0.0) > Lower Right ( 8792.0,10070.0) > Center ( 4396.0, 5035.0) > Band 1 Block=1024x1024 Type=UInt16, ColorInterp=Undefined > Metadata: > POLARIMETRIC_INTERP=HH > Band 2 Block=1024x1024 Type=UInt16, ColorInterp=Undefined > Metadata: > POLARIMETRIC_INTERP=HV _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
