[gdal-dev] RADARSAT2 NITF RPC georeferencing
Even Rouault
even.rouault at mines-paris.org
Sun Jun 28 13:29:38 EDT 2009
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
More information about the gdal-dev
mailing list