Hi Frank,

I only want to note that the default is bilinear interpolation of the DEM.

Part of the code: const char *pszDEMInterpolation = CSLFetchNameValueDef( papszOptions, "RPC_DEMINTERPOLATION", "bilinear" );

Best regards,
    Dmitriy

16.01.2013 22:47, Frank Warmerdam ?????:
Yehiyam,

I do not see why you think the AREA_OR_POINT value ought to have any influence on how DEM values are interpolated.

Note that a PixelIsPoint GeoTIFF file will already result in a GDAL GeoTransform that is offset by half a pixel. The GDAL GeoTransform always treats pixels as areas, and GDAL interpolators - presumably including the for RPCs - should treat the pixel value as being at the center of the pixel area for the purpose of interpolation (regardless of whether that value is point or area).

On the other hand, I am denying the possibility that the RPC DEM interpolation is always off by half a pixel. I haven't actually looked closely at that code lately and the RPC code is not so very well tested and validated.

Best regards,
Frank


On Wed, Jan 16, 2013 at 4:40 AM, Livneh Yehiyam <ye...@rafael.co.il <mailto:ye...@rafael.co.il>> wrote:


    Hi

    I have a question about the RPC transform implemented in gdal_rpc.cpp.

    If a dem file is specified in the transform options (RPC_DEM) the
    file is sampled to retrieve the elevation of the transformed point.

    In the case where the dem file is a DTED file (which defines
    AREA_OR_POINT to POINT), shouldn't the pixel/line values in the
    dem file be offset by half a pixel to get the correct height?

    Thanks

    Yehiyam


    ________________________________________
    This message (including any attachments) issued by RAFAEL-
    ADVANCED DEFENSE SYSTEMS LTD. (hereinafter "RAFAEL") contains
    confidential information intended for a specific individual and
    purpose, may constitute information that is privileged or
    confidential or otherwise protected from disclosure. If you are
    not the intended recipient, you should contact us immediately and
    thereafter delete this message from your system. You are hereby
    notified that any disclosure, copying, dissemination, distribution
    or forwarding of this message, or the taking of any action based
    on it, is strictly prohibited. If you have received this e-mail in
    error, please notify us immediately by e-mail
    mailto:law...@rafael.co.il <mailto:law...@rafael.co.il> and
    completely delete or destroy any and all electronic or other
    copies of the original message and any attachments thereof.

    _______________________________________________
    gdal-dev mailing list
    gdal-dev@lists.osgeo.org <mailto:gdal-dev@lists.osgeo.org>
    http://lists.osgeo.org/mailman/listinfo/gdal-dev




--
---------------------------------------+--------------------------------------
I set the clouds in motion - turn up | Frank Warmerdam, warmer...@pobox.com <mailto:warmer...@pobox.com> light and sound - activate the windows | http://pobox.com/~warmerdam <http://pobox.com/%7Ewarmerdam>
and watch the world go round - Rush    | Geospatial Software Developer


_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/gdal-dev

Reply via email to