Computing tile origin in dst_srs failed. Starting datacube failed. Starting datacube(s) failed.

66 views
Skip to first unread message

moh713...@gmail.com

unread,
Nov 13, 2020, 11:58:54 AM11/13/20
to FORCE
Dear David, I have gone through this documentation : 
trying to understand the issue regarding this issue :
Computing tile origin in dst_srs failed. 
Starting datacube failed. 
Starting datacube(s) failed.
arising due to using non-built-in projection.
in this document, it has been mentioned that :
Note that coordinates may be undefined if they are too far away from the origin of the coordinate system (see Transverse Mercator example below); in this case the algorithm will fail. It is good practice to use a point of origin that is relatively close to the study area. The allowed distance varies greatly with different output projections.
the documentation also referred to to see figures 13,14,15 which are not available.

I am asking again how can I set the projection without getting this error?
thanks in advance

thie...@gmail.com

unread,
Nov 25, 2020, 4:24:13 AM11/25/20
to FORCE
Hi Mohammed,
looks like an issue with gdaltransform.
Try running gdaltransform manually in the shell with your desired settings. Does is work? So for example
gdaltransform -s_srs EPSG:28992 -t_srs EPSG:31370

david frantz

unread,
Dec 4, 2020, 7:56:51 AM12/4/20
to FORCE
Do you happen to use anaconda?
Best,
David

moh713...@gmail.com

unread,
Dec 4, 2020, 9:42:03 AM12/4/20
to FORCE
Hi David,
Yes, I have anaconda installed on my pc.
thanks

david frantz

unread,
Dec 4, 2020, 9:47:23 AM12/4/20
to FORCE
I have recently realized that anaconda messes with the system paths and thus corrupts the libraries linked with FORCE...
Please disable your anaconda environment, and try again. 
Thanks,
David

moh713...@gmail.com

unread,
Dec 5, 2020, 2:38:06 AM12/5/20
to FORCE
Hi David,
I have disable anaconda but the issue is still arising.
Do you think that I should reinstall the required open-source dependencies?
Regards
Mohammed

david frantz

unread,
Dec 10, 2020, 3:01:33 AM12/10/20
to FORCE
Hi Mohammed.
if you suspect that this is a dependency problem, please try using the docker version. 
If the problem persists, it is very likely that there is a small issue in your parameterization.
Best,
David

moh713...@gmail.com

unread,
Dec 10, 2020, 2:55:28 PM12/10/20
to FORCE
Dear David,
My question was based on your suggestion to disable anaconda.
AS I already mentioned that the problem arising when I try to use a non-built-in projection where the study area located in more than 3 UTM zones.
could you show me the parameterization of the data cube which you already used successfully non-built-in projection similar to my case?
in other words, how the parameter "PROJECTION " should be written down? Do you have a similar case?
DO_REPROJ = TRUE
DO_TILE = TRUE
FILE_TILE = NULL #/home/localuser/force/mis/allow-list_landsat_equal_area_6000.txt
TILE_SIZE = 6000
BLOCK_SIZE = 3000
RESOLUTION_LANDSAT = 30
RESOLUTION_SENTINEL2 = 10
ORIGIN_LON = 40.9781 
ORIGIN_LAT = 19.1660
PROJECTION = PROJCRS["ProjWiz_Custom_Albers",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Albers Equal Area",ID["EPSG",9822]],PARAMETER["Latitude of false origin",15.5248447,
ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",48.1695557,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",13.2047435,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",17.844946,ANGLEUNIT["Degree",0.0174532925199433],
ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,
LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,
ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]

RESAMPLING = CC

moh713...@gmail.com

unread,
Dec 10, 2020, 3:48:35 PM12/10/20
to FORCE
Dear David,

these are two different parameterizations of the data cube  : the first one works well with force-level2

# DATA CUBES
# ------------------------------------------------------------------------
DO_REPROJ = TRUE
DO_TILE = TRUE
FILE_TILE = /home/localuser/force/mis/allow_list/socotra/allow-list_landsat_equal_area_6000.txt
TILE_SIZE = 60000
BLOCK_SIZE = 6000
RESOLUTION_LANDSAT = 30
RESOLUTION_SENTINEL2 = 10
ORIGIN_LON =  53.25506500842682     
ORIGIN_LAT =  12.75565804078176  
PROJECTION = PROJCRS["ProjWiz_Custom_Cylindrical_Equal_Area",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Lambert Cylindrical Equal Area (Spherical)",ID["EPSG",9834]],PARAMETER["Latitude of 1st standard parallel",6.3630421,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Longitude of natural origin",53.9222717,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]

RESAMPLING = CC



while this one arises  with this issue : 
Computing tile origin in dst_srs failed. 
Starting datacube failed. 
Starting datacube(s) failed.

# DATA CUBES
# ------------------------------------------------------------------------
DO_REPROJ = TRUE
DO_TILE = TRUE
FILE_TILE = NULL 
TILE_SIZE = 6000
BLOCK_SIZE = 3000
RESOLUTION_LANDSAT = 30
RESOLUTION_SENTINEL2 = 10
ORIGIN_LON = 40.9781 
ORIGIN_LAT = 19.1660
PROJECTION = PROJCRS["ProjWiz_Custom_Albers",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Albers Equal Area",ID["EPSG",9822]],PARAMETER["Latitude of false origin",15.5248447,
ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",48.1695557,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",13.2047435,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",17.844946,ANGLEUNIT["Degree",0.0174532925199433],
ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,
LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,
ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]

RESAMPLING = CC

Do you see any difference between the two parameterizations?
the only difference is the projection !!!

david frantz

unread,
Dec 11, 2020, 2:41:36 AM12/11/20
to FORCE
Mohammed..

First of all, the 2nd projection is not written in one line. This can't work.
Secondly, and as said multiple times before, you need to validate your projection with gdalsrsinfo! Both of your projections are not valid (see gdalsrsinfo test below)!
I cannot follow why the 1st one works for you. They both don't work for me (see gdaltransform test below).
How did you come up with this projection? I highly recommend to use one from here: https://spatialreference.org/ref/epsg/
However, you still need to validate it to properly work with gdal.

Best,
David


gdalsrsinfo -V 'PROJCRS["ProjWiz_Custom_Cylindrical_Equal_Area",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Lambert Cylindrical Equal Area (Spherical)",ID["EPSG",9834]],PARAMETER["Latitude of 1st standard parallel",6.3630421,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Longitude of natural origin",53.9222717,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'

ERROR 1: ERROR - failed to load SRS definition from PROJCRS["ProjWiz_Custom_Cylindrical_Equal_Area",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Lambert Cylindrical Equal Area (Spherical)",ID["EPSG",9834]],PARAMETER["Latitude of 1st standard parallel",6.3630421,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Longitude of natural origin",53.9222717,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]


gdalsrsinfo -V 'PROJCRS["ProjWiz_Custom_Albers",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Albers Equal Area",ID["EPSG",9822]],PARAMETER["Latitude of false origin",15.5248447,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",48.1695557,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",13.2047435,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",17.844946,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'

ERROR 1: ERROR - failed to load SRS definition from PROJCRS["ProjWiz_Custom_Albers",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Albers Equal Area",ID["EPSG",9822]],PARAMETER["Latitude of false origin",15.5248447,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",48.1695557,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",13.2047435,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",17.844946,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]


gdaltransform -s_srs 'EPSG:4326' -t_srs 'PROJCRS["ProjWiz_Custom_Cylindrical_Equal_Area",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Lambert Cylindrical Equal Area (Spherical)",ID["EPSG",9834]],PARAMETER["Latitude of 1st standard parallel",6.3630421,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Longitude of natural origin",53.9222717,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'

ERROR 1: Translating source or target SRS failed:
PROJCRS["ProjWiz_Custom_Cylindrical_Equal_Area",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Lambert Cylindrical Equal Area (Spherical)",ID["EPSG",9834]],PARAMETER["Latitude of 1st standard parallel",6.3630421,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Longitude of natural origin",53.9222717,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]


gdaltransform -s_srs 'EPSG:4326' -t_srs 'PROJCRS["ProjWiz_Custom_Albers",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Albers Equal Area",ID["EPSG",9822]],PARAMETER["Latitude of false origin",15.5248447,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",48.1695557,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",13.2047435,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",17.844946,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'

ERROR 1: Translating source or target SRS failed:
PROJCRS["ProjWiz_Custom_Albers",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Albers Equal Area",ID["EPSG",9822]],PARAMETER["Latitude of false origin",15.5248447,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",48.1695557,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",13.2047435,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",17.844946,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]

moh713...@gmail.com

unread,
Dec 11, 2020, 5:49:00 AM12/11/20
to FORCE
Hi David,

First of all, the 2nd projection is not written in one line. This can't work
Could you tell me please how it can be in one line?? how do you discover this as I followed the same way with the first one !!

Secondly, and as said multiple times before, you need to validate your projection with gdalsrsinfo! Both of your projections are not valid (see gdalsrsinfo test below)!
I have already validated them.
It is strange !! for validation using the same command you use, I get successful results as you can see below 

gdalsrsinfo -V 'PROJCRS["ProjWiz_Custom_Cylindrical_Equal_Area",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Lambert Cylindrical Equal Area (Spherical)",ID["EPSG",9834]],PARAMETER["Latitude of 1st standard parallel",6.3630421,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Longitude of natural origin",53.9222717,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'

Validate Succeeds

PROJ.4 : +proj=cea +lat_ts=6.3630421 +lon_0=53.9222717 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs

OGC WKT2:2018 :
PROJCRS["ProjWiz_Custom_Cylindrical_Equal_Area",
    BASEGEOGCRS["WGS 84",
        DATUM["World Geodetic System 1984",
            ELLIPSOID["WGS 84",6378137,298.257223563,
                LENGTHUNIT["metre",1]],
            ID["EPSG",6326]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["Degree",0.0174532925199433]]],
    CONVERSION["unnamed",
        METHOD["Lambert Cylindrical Equal Area (Spherical)",
            ID["EPSG",9834]],
        PARAMETER["Latitude of 1st standard parallel",6.3630421,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8823]],
        PARAMETER["Longitude of natural origin",53.9222717,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["False easting",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["(E)",east,
            ORDER[1],
             LENGTHUNIT["metre",1,
                ID["EPSG",9001]]],
        AXIS["(N)",north,
            ORDER[2],
            LENGTHUNIT["metre",1,
                ID["EPSG",9001]]]]

gdalsrsinfo -V 'PROJCRS["ProjWiz_Custom_Albers",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Albers Equal Area",ID["EPSG",9822]],PARAMETER["Latitude of false origin",15.5248447,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",48.1695557,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",13.2047435,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",17.844946,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'

Validate Succeeds

PROJ.4 : +proj=aea +lat_0=15.5248447 +lon_0=48.1695557 +lat_1=13.2047435 +lat_2=17.844946 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs

OGC WKT2:2018 :
PROJCRS["ProjWiz_Custom_Albers",
    BASEGEOGCRS["WGS 84",
        DATUM["World Geodetic System 1984",
            ELLIPSOID["WGS 84",6378137,298.257223563,
                LENGTHUNIT["metre",1]],
            ID["EPSG",6326]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["Degree",0.0174532925199433]]],
    CONVERSION["unnamed",
        METHOD["Albers Equal Area",
            ID["EPSG",9822]],
        PARAMETER["Latitude of false origin",15.5248447,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8821]],
        PARAMETER["Longitude of false origin",48.1695557,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8822]],
        PARAMETER["Latitude of 1st standard parallel",13.2047435,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8823]],
        PARAMETER["Latitude of 2nd standard parallel",17.844946,
            ANGLEUNIT["Degree",0.0174532925199433],
            ID["EPSG",8824]],
        PARAMETER["Easting at false origin",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8826]],
        PARAMETER["Northing at false origin",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8827]]],
    CS[Cartesian,2],
        AXIS["(E)",east,
            ORDER[1],
            LENGTHUNIT["metre",1,
                ID["EPSG",9001]]],
        AXIS["(N)",north,
            ORDER[2],
            LENGTHUNIT["metre",1,
                ID["EPSG",9001]]]]

I cannot follow why the 1st one works for you. They both don't work for me (see gdaltransform test below).

Regarding the  "gdaltransform test"
I do not see any input coordinates you enter to test the command as this is one of the requirements of testing !!
when  I used the same command you already used, it gives me the choice to put some value to check as it is documented here: https://gdal.org/programs/gdaltransform.html

Could you tell me what are the values that you enter to get that error after these two commands :

gdaltransform -s_srs 'EPSG:4326' -t_srs 'PROJCRS["ProjWiz_Custom_Cylindrical_Equal_Area",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Lambert Cylindrical Equal Area (Spherical)",ID["EPSG",9834]],PARAMETER["Latitude of 1st standard parallel",6.3630421,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Longitude of natural origin",53.9222717,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8802]],PARAMETER["False easting",0,LENGTHUNIT["metre",1],ID["EPSG",8806]],PARAMETER["False northing",0,LENGTHUNIT["metre",1],ID["EPSG",8807]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHhttps://gdal.org/programs/gdaltransform.htmlUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'

gdaltransform -s_srs 'EPSG:4326' -t_srs 'PROJCRS["ProjWiz_Custom_Albers",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Albers Equal Area",ID["EPSG",9822]],PARAMETER["Latitude of false origin",15.5248447,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",48.1695557,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",13.2047435,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",17.844946,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'

How did you come up with this projection? I highly recommend to use one from here: https://spatialreference.org/ref/epsg/
I used an online tool that can give the best equal-area projection for any area and then I validated it using gdal as documented in FORCE  totorial and as I showed that how it has been successfully validated.
 using the appropriate projection is crucial in getting the less distortion and this is very clear in FORCE documentation here :
and this is why I want the most appropriate projection as can as possible.

Hopefully, this issue be solved and you forgive me for all my questions and issues >
Regards

david frantz

unread,
Dec 11, 2020, 6:04:33 AM12/11/20
to FORCE
Well, it doesn't work for me. This could be related to different GDAL versions, I don't know. 

On my machine, the gdaltransform fails before I can input any coordinates. If it works for you, try your  ORIGIN_LON / ORIGIN_LAT coordinates and see whether the transformation works for these points. 
My best bet is: your online tool optimized the projection on a very, very, small area, thus the origin coordinates are outside of the area where the projection is valid.

If you need to make this projection work, I suggest you inquire at GDAL developers for help on this.

If the key aspect of your projection requirements is "Equal Area" I suggest to simply use "GLANCE" as projection.  There is no distortion in area when using this projection.

Please note: since I cannot reproduce your problem, I am not able to help you further with this...  

Best,
David

moh713...@gmail.com

unread,
Dec 11, 2020, 7:51:03 AM12/11/20
to FORCE
Well, it doesn't work for me. This could be related to different GDAL versions, I don't know. 
this is mine :
GDAL 3.0.4, released 2020/01/28
What about you ? is there is a difference between the two versions that  can lead to this issue??!

If it works for you, try your  ORIGIN_LON / ORIGIN_LAT coordinates and see whether the transformation works for these points. 
I tried as follow :
gdaltransform -s_srs 'EPSG:4326' -t_srs 'PROJCRS["ProjWiz_Custom_Albers",BASEGEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]],ID["EPSG",6326]],PRIMEM["Greenwich",0,ANGLEUNIT["Degree",0.0174532925199433]]],CONVERSION["unnamed",METHOD["Albers Equal Area",ID["EPSG",9822]],PARAMETER["Latitude of false origin",15.5248447,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",48.1695557,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",13.2047435,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",17.844946,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8824]],PARAMETER["Easting at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8826]],PARAMETER["Northing at false origin",0,LENGTHUNIT["metre",1],ID["EPSG",8827]]],CS[Cartesian,2],AXIS["(E)",east,ORDER[1],LENGTHUNIT["metre",1,ID["EPSG",9001]]],AXIS["(N)",north,ORDER[2],LENGTHUNIT["metre",1,ID["EPSG",9001]]]]'
40.9781 19.1660

and it works to give this result : 
-757226.018925644 415747.304410907 0

My best bet is: your online tool optimized the projection on a very, very, small area, thus the origin coordinates are outside of the area where the projection is valid.
If you need to make this projection work, I suggest you inquire at GDAL developers for help on this.
I will check and try this again and report this here as can as possible.

If the key aspect of your projection requirements is "Equal Area" I suggest to simply use "GLANCE" as projection.  There is no distortion in area when using this projection.
 I will try to compare GLANCE  projection for a KNOWN AREA AND SEE THAT and compare the results with the customised projection I want to work with and see what happens.

Please note: since I cannot reproduce your problem, I am not able to help you further with this...  
ok, can you please tell me what is your comments about my last inquiery regarding this point in the last conversation :

First of all, the 2nd projection is not written in one line. This can't work
Could you tell me please how it can be in one line?? how do you discover this as I followed the same way with the first one !! could you rewrite it to be in one line please like the first one which worked normally, please?

this is because I have faced the same problem with the first projection before and I rewrote the projection and it worked normally but  I do not know how that happened and you told me today the second projection is not in one line.

Thanks again 

david frantz

unread,
Dec 11, 2020, 7:54:08 AM12/11/20
to FORCE
I am really sorry Mohammed, this thread takes too much of my time. Please use one of the built-in projections. Best, David

moh713...@gmail.com

unread,
Dec 11, 2020, 8:30:31 AM12/11/20
to FORCE
I really appreciate this and when I put an issue here, this is only to report something that is related to the purpose of the creation of this group only and  I do not expect to get the results urgently.
 I also thought that my final offer regarding what written below will not take a lot of time :
First of all, the 2nd projection is not written in one line. This can't work
Could you tell me please how it can be in one line?? how do you discover this as I followed the same way with the first one !! could you rewrite it to be in one line please like the first one which worked normally, please?
this is because I have faced the same problem with the first projection before and I rewrote the projection and it worked normally but  I do not know how that happened and you told me today the second projection is not in one line.

Thanks again Dr. David for all your efforts and your time
Reply all
Reply to author
Forward
0 new messages