TILE_SIZE must be a multiple of RESOLUTION

[moh713...@gmail.com]

Dear David,

I am running the  : force-level2 /home/localuser/force/param/yem.prm with  FORCE v. 3.5.2  on UBUNTU 18.04 to produce ARD for landsat with A custom projection >

the log file give me this : 

TILE_SIZE must be a multiple of RESOLUTION. If DO_REPROJ = FALSE, the image resolution overrides the RESOLUTION parameter. Starting datacube failed.

Starting datacube(s) failed.

the data cube parameters as this :

# DATA CUBES

# ------------------------------------------------------------------------

DO_REPROJ = TRUE

DO_TILE = TRUE

FILE_TILE = NULL #/home/localuser/force/mis/allow-list_landsat_equal_area_6000.txt

TILE_SIZE = 5000

BLOCK_SIZE = 2500

RESOLUTION_LANDSAT = 30

RESOLUTION_SENTINEL2 = 10

ORIGIN_LON = 40.97805253919763

ORIGIN_LAT = 19.16599361925848  

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.5729207,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8821]],PARAMETER["Longitude of false origin",48.0541992,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8822]],PARAMETER["Latitude of 1st standard parallel",13.274483,ANGLEUNIT["Degree",0.0174532925199433],ID["EPSG",8823]],PARAMETER["Latitude of 2nd standard parallel",17.8713583,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

thanks in advance

[moh713...@gmail.com]

sorry this is as a result of writing  :TILE_SIZE = 5000

this should be a multiple of landsat resolution  :30
when I changed it to  TILE_SIZE =6000,  the error disappeared

[Stefan]

TILE_SIZE = 5000

RESOLUTION = 30

5000 / 30 = 166,666667

Error:  TILE_SIZE must be a multiple of RESOLUTION.

You are trying to fit 30 m pixels into a 5000 m tile. Imagine you have a kitchen that is 5000 cm wide and want redo the floor with 30 cm tiles. You won't be able to fit the 30 cm tiles without cutting the last one.
We can't just cut pixels, so the raster size needs to be divisible by the pixel size.

[moh713...@gmail.com]

Thank you very much for your clarification Stefan.

I have read the documentation and reported this here .

Now I am facing a new kind of problem which I already solved. 

this problem is reported in the log file as :

Transformation failed.

Computing tile origin in dst_srs failed. Starting datacube failed.

Starting datacube(s) failed.

I tried to solve it following the same steps when I already faced the problem and reported that here in a previous post, unfortunately no progress happened.

this issue always arise when using a non-built-in projection.

this projection is  : 

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]]]]

the all parameters in the data cube are :

# DATA CUBES

# ------------------------------------------------------------------------

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

[david frantz]

Dear Mohammed,

"Now I am facing a new kind of problem which I already solved.   "

I am not sure if I understand you correctly, but please do not use this group as a personal log for any problem you run into..

Also, I am wodering. Several times now, you arrived at a working solution. But still you continue to change your parameterization. Why not stick to a parameterization that works for you?

Best,

David

[moh713...@gmail.com]

Dear David,

Also, I am wodering. Several times now, you arrived at a working solution. But still you continue to change your parameterization. Why not stick to a parameterization that works for you?

I was doing testing for a small area and now I changed the projection for a new large area .

 As I followed the same solution when I faced the same problem for the small  testing area, the error is still arising for the new large area !!

hopefully this is very clear now.

 thanks