Re: [gdal-dev] feet converting to meters unexpectedly

2023-01-09 Thread Even Rouault 


Le 09/01/2023 à 20:22, Rahkonen Jukka a écrit :


Hi,

I noticed that the survey foot is now deemed obsolete 
https://www.nist.gov/pml/us-surveyfoot since 1^st of January. 
Should/could there be warnings in GDAL/Proj about that? Or will also 
the corresponding entries in the EPSG database become obsolete in time?


The US government rightly wants to fix the future, but it can't fix the 
past. I bet that the use of CRS using US Survey Foot will outlive all of 
us. IMHO, I don't think it is GDAL/Proj responsibility to do anything 
special regarding that. Probably that tagging those CRS in the EPSG 
database will indeed be the next action, but I believe this should be 
with care as it could have surprising effects. Probably the supersession 
mechanism would be more appropriate than deprecation (deprecation in 
EPSG is often used to tag broken entries. there is nothing broken with 
US Survey Foot)


Even


-Jukka Rahkonen-

*Lähettäjä:*gdal-dev  *Puolesta 
*Even Rouault

*Lähetetty:* maanantai 9. tammikuuta 2023 21.10
*Vastaanottaja:* Kirk Waters - NOAA Federal ; 
gdal-dev@lists.osgeo.org

*Aihe:* Re: [gdal-dev] feet converting to meters unexpectedly

Hi,

This is indeed https://github.com/OSGeo/gdal/issues/6839 
<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2FOSGeo%2Fgdal%2Fissues%2F6839&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332787798%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=50fmqgZrGZQQdpKoc1fOlTL0yhI2z1Wf9Dabs0KwExs%3D&reserved=0> 
which has been fixed in GDAL 3.6.1.


Even

Le 09/01/2023 à 19:48, Kirk Waters - NOAA Federal via gdal-dev a écrit :

I think I've gotten this solved for my purposes now, but I'm
reporting it here in case anyone else might run across issues
using gdal in a conda environment. If I make a conda environment
containing:

python

pdal=2.4.2

gdal

(which ends up with gdal version 3.5.1, pdal version 2.4.2), I get
a weird behavior where

gdalwarp input.tif output.tif

will change the values in the raster as if it is converting from
feet to meters. My example raster where I see this has
georeferencing of EPSG:3521+6360 and is a 32-bit floating point
DEM. My use case would have more arguments to gdalwarp, but I
simplified to the minimum that showed the problem.

If I change the environment package list for conda to not specify
the pdal version, I'll get pdal version 2.4.3 and gdal version
3.6.2 and the raster will convert correctly.

I did note that in a setup where it converts the values, I see
this warning:

Warning 1: Unknown units=US survey foot


Kirk Waters, PhD

NOAA Office for Coastal Management

Applied Sciences Program

‪Phone: 843-284-6962 (New as of 2022)

coast.noaa.gov/digitalcoast

<https://eur06.safelinks.protection.outlook.com/?url=http%3A%2F%2Fcoast.noaa.gov%2Fdigitalcoast&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332787798%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=mFBiuSZzlkEop5b4qQLEcKr1Xmmmu1sGw0NxxOE0dj4%3D&reserved=0>



___

gdal-dev mailing list

gdal-dev@lists.osgeo.org

https://lists.osgeo.org/mailman/listinfo/gdal-dev  
<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.osgeo.org%2Fmailman%2Flistinfo%2Fgdal-dev&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332787798%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=3ZnHOEI50I0QfHU5qO6IIOZlEMB4WFJRaBhAaB0vW%2Bk%3D&reserved=0>

--
http://www.spatialys.com  
<https://eur06.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.spatialys.com%2F&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332944041%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=3XoD5Rxng%2BozCi%2FyIAClJLUZHwuJnAhXM%2BAjTn4QtrM%3D&reserved=0>
My software is free, but my time generally not.


--
http://www.spatialys.com
My software is free, but my time generally not.
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] feet converting to meters unexpectedly

2023-01-09 Thread Kirk Waters - NOAA Federal via gdal-dev 
It's not exactly that it is obsolete. A better way to think of it might be
that it will not be part of the state plane 2022 definitions, which won't
be out until about 2025. It will always be part of the state plane 1983
definitions.

Kirk


On Mon, Jan 9, 2023 at 2:23 PM Rahkonen Jukka <
jukka.rahko...@maanmittauslaitos.fi> wrote:

> Hi,
>
>
>
> I noticed that the survey foot is now deemed obsolete
> https://www.nist.gov/pml/us-surveyfoot since 1st of January. Should/could
> there be warnings in GDAL/Proj about that? Or will also the corresponding
> entries in the EPSG database become obsolete in time?
>
>
>
> -Jukka Rahkonen-
>
>
>
>
>
> *Lähettäjä:* gdal-dev  *Puolesta *Even
> Rouault
> *Lähetetty:* maanantai 9. tammikuuta 2023 21.10
> *Vastaanottaja:* Kirk Waters - NOAA Federal ;
> gdal-dev@lists.osgeo.org
> *Aihe:* Re: [gdal-dev] feet converting to meters unexpectedly
>
>
>
> Hi,
>
> This is indeed https://github.com/OSGeo/gdal/issues/6839
> <https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2FOSGeo%2Fgdal%2Fissues%2F6839&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332787798%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=50fmqgZrGZQQdpKoc1fOlTL0yhI2z1Wf9Dabs0KwExs%3D&reserved=0>
> which has been fixed in GDAL 3.6.1.
>
> Even
>
> Le 09/01/2023 à 19:48, Kirk Waters - NOAA Federal via gdal-dev a écrit :
>
> I think I've gotten this solved for my purposes now, but I'm reporting it
> here in case anyone else might run across issues using gdal in a conda
> environment. If I make a conda environment containing:
>
> python
>
> pdal=2.4.2
>
> gdal
>
> (which ends up with gdal version 3.5.1, pdal version 2.4.2), I get a weird
> behavior where
>
> gdalwarp input.tif output.tif
>
> will change the values in the raster as if it is converting from feet to
> meters. My example raster where I see this has georeferencing of
> EPSG:3521+6360 and is a 32-bit floating point DEM. My use case would have
> more arguments to gdalwarp, but I simplified to the minimum that showed the
> problem.
>
>
>
> If I change the environment package list for conda to not specify the pdal
> version, I'll get pdal version 2.4.3 and gdal version 3.6.2 and the raster
> will convert correctly.
>
>
>
> I did note that in a setup where it converts the values, I see this
> warning:
>
> Warning 1: Unknown units=US survey foot
>
>
> Kirk Waters, PhD
>
> NOAA Office for Coastal Management
>
> Applied Sciences Program
>
> ‪Phone: 843-284-6962 (New as of 2022)
>
> coast.noaa.gov/digitalcoast
> <https://eur06.safelinks.protection.outlook.com/?url=http%3A%2F%2Fcoast.noaa.gov%2Fdigitalcoast&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332787798%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=mFBiuSZzlkEop5b4qQLEcKr1Xmmmu1sGw0NxxOE0dj4%3D&reserved=0>
>
>
>
>
>
> ___
>
> gdal-dev mailing list
>
> gdal-dev@lists.osgeo.org
>
> https://lists.osgeo.org/mailman/listinfo/gdal-dev 
> <https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.osgeo.org%2Fmailman%2Flistinfo%2Fgdal-dev&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332787798%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=3ZnHOEI50I0QfHU5qO6IIOZlEMB4WFJRaBhAaB0vW%2Bk%3D&reserved=0>
>
> --
>
> http://www.spatialys.com 
> <https://eur06.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.spatialys.com%2F&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332944041%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=3XoD5Rxng%2BozCi%2FyIAClJLUZHwuJnAhXM%2BAjTn4QtrM%3D&reserved=0>
>
> My software is free, but my time generally not.
>
>
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] feet converting to meters unexpectedly

2023-01-09 Thread Rahkonen Jukka 
Hi,

I noticed that the survey foot is now deemed obsolete 
https://www.nist.gov/pml/us-surveyfoot since 1st of January. Should/could there 
be warnings in GDAL/Proj about that? Or will also the corresponding entries in 
the EPSG database become obsolete in time?

-Jukka Rahkonen-


Lähettäjä: gdal-dev  Puolesta Even Rouault
Lähetetty: maanantai 9. tammikuuta 2023 21.10
Vastaanottaja: Kirk Waters - NOAA Federal ; 
gdal-dev@lists.osgeo.org
Aihe: Re: [gdal-dev] feet converting to meters unexpectedly


Hi,

This is indeed 
https://github.com/OSGeo/gdal/issues/6839<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2FOSGeo%2Fgdal%2Fissues%2F6839&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332787798%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=50fmqgZrGZQQdpKoc1fOlTL0yhI2z1Wf9Dabs0KwExs%3D&reserved=0>
 which has been fixed in GDAL 3.6.1.

Even
Le 09/01/2023 à 19:48, Kirk Waters - NOAA Federal via gdal-dev a écrit :
I think I've gotten this solved for my purposes now, but I'm reporting it here 
in case anyone else might run across issues using gdal in a conda environment. 
If I make a conda environment containing:
python
pdal=2.4.2
gdal
(which ends up with gdal version 3.5.1, pdal version 2.4.2), I get a weird 
behavior where
gdalwarp input.tif output.tif
will change the values in the raster as if it is converting from feet to 
meters. My example raster where I see this has georeferencing of EPSG:3521+6360 
and is a 32-bit floating point DEM. My use case would have more arguments to 
gdalwarp, but I simplified to the minimum that showed the problem.

If I change the environment package list for conda to not specify the pdal 
version, I'll get pdal version 2.4.3 and gdal version 3.6.2 and the raster will 
convert correctly.

I did note that in a setup where it converts the values, I see this warning:
Warning 1: Unknown units=US survey foot

Kirk Waters, PhD
NOAA Office for Coastal Management
Applied Sciences Program
‪Phone: 843-284-6962 (New as of 2022)
coast.noaa.gov/digitalcoast<https://eur06.safelinks.protection.outlook.com/?url=http%3A%2F%2Fcoast.noaa.gov%2Fdigitalcoast&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332787798%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=mFBiuSZzlkEop5b4qQLEcKr1Xmmmu1sGw0NxxOE0dj4%3D&reserved=0>




___

gdal-dev mailing list

gdal-dev@lists.osgeo.org<mailto:gdal-dev@lists.osgeo.org>

https://lists.osgeo.org/mailman/listinfo/gdal-dev<https://eur06.safelinks.protection.outlook.com/?url=https%3A%2F%2Flists.osgeo.org%2Fmailman%2Flistinfo%2Fgdal-dev&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332787798%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=3ZnHOEI50I0QfHU5qO6IIOZlEMB4WFJRaBhAaB0vW%2Bk%3D&reserved=0>

--

http://www.spatialys.com<https://eur06.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.spatialys.com%2F&data=05%7C01%7Cjukka.rahkonen%40maanmittauslaitos.fi%7Cff99c3d359784e7f45ce08daf2752dad%7Cc4f8a63255804a1c92371d5a571b71fa%7C0%7C1%7C63802332944041%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=3XoD5Rxng%2BozCi%2FyIAClJLUZHwuJnAhXM%2BAjTn4QtrM%3D&reserved=0>

My software is free, but my time generally not.
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] feet converting to meters unexpectedly

2023-01-09 Thread Even Rouault 
I should also mention that you can probably workaround the issue by 
adding the -novshift switch to your gdalwarp command line: 
https://gdal.org/programs/gdalwarp.html#cmdoption-gdalwarp-novshift


Le 09/01/2023 à 20:10, Even Rouault a écrit :


Hi,

This is indeed https://github.com/OSGeo/gdal/issues/6839 which has 
been fixed in GDAL 3.6.1.


Even

Le 09/01/2023 à 19:48, Kirk Waters - NOAA Federal via gdal-dev a écrit :
I think I've gotten this solved for my purposes now, but I'm 
reporting it here in case anyone else might run across issues using 
gdal in a conda environment. If I make a conda environment containing:

python
pdal=2.4.2
gdal
(which ends up with gdal version 3.5.1, pdal version 2.4.2), I get a 
weird behavior where

gdalwarp input.tif output.tif
will change the values in the raster as if it is converting from feet 
to meters. My example raster where I see this has georeferencing of 
EPSG:3521+6360 and is a 32-bit floating point DEM. My use case would 
have more arguments to gdalwarp, but I simplified to the minimum that 
showed the problem.


If I change the environment package list for conda to not specify the 
pdal version, I'll get pdal version 2.4.3 and gdal version 3.6.2 and 
the raster will convert correctly.


I did note that in a setup where it converts the values, I see this 
warning:

Warning 1: Unknown units=US survey foot

Kirk Waters, PhD
NOAA Office for Coastal Management
Applied Sciences Program
‪Phone: 843-284-6962‬ (New as of 2022)
coast.noaa.gov/digitalcoast 



___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

--
http://www.spatialys.com
My software is free, but my time generally not.

___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev


--
http://www.spatialys.com
My software is free, but my time generally not.
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev

Re: [gdal-dev] feet converting to meters unexpectedly

2023-01-09 Thread Even Rouault 

Hi,

This is indeed https://github.com/OSGeo/gdal/issues/6839 which has been 
fixed in GDAL 3.6.1.


Even

Le 09/01/2023 à 19:48, Kirk Waters - NOAA Federal via gdal-dev a écrit :
I think I've gotten this solved for my purposes now, but I'm reporting 
it here in case anyone else might run across issues using gdal in a 
conda environment. If I make a conda environment containing:

python
pdal=2.4.2
gdal
(which ends up with gdal version 3.5.1, pdal version 2.4.2), I get a 
weird behavior where

gdalwarp input.tif output.tif
will change the values in the raster as if it is converting from feet 
to meters. My example raster where I see this has georeferencing of 
EPSG:3521+6360 and is a 32-bit floating point DEM. My use case would 
have more arguments to gdalwarp, but I simplified to the minimum that 
showed the problem.


If I change the environment package list for conda to not specify the 
pdal version, I'll get pdal version 2.4.3 and gdal version 3.6.2 and 
the raster will convert correctly.


I did note that in a setup where it converts the values, I see this 
warning:

Warning 1: Unknown units=US survey foot

Kirk Waters, PhD
NOAA Office for Coastal Management
Applied Sciences Program
‪Phone: 843-284-6962‬ (New as of 2022)
coast.noaa.gov/digitalcoast 



___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev


--
http://www.spatialys.com
My software is free, but my time generally not.
___
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev