TILE_SIZE must be a multiple of RESOLUTION

53 views
Skip to first unread message

moh713...@gmail.com

unread,
Nov 5, 2020, 3:50:31 AM11/5/20
to FORCE
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

unread,
Nov 5, 2020, 6:26:38 AM11/5/20
to FORCE
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

unread,
Nov 5, 2020, 8:15:35 AM11/5/20
to FORCE
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

unread,
Nov 5, 2020, 8:48:47 AM11/5/20
to FORCE
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

unread,
Nov 9, 2020, 9:27:47 AM11/9/20
to FORCE
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

unread,
Nov 10, 2020, 9:04:22 AM11/10/20
to FORCE
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
Reply all
Reply to author
Forward
0 new messages