hi, I'm doing a land cover map for Serbia and I'm stuck at 15:34 mark, where my min and max values are 1.000000 and 1.000008 respectively. I've been using chatgpt to troubleshoot my issue for 2 hours now but with no luck at all. what could be the issue?
Hi @@KatarinaaBlaa Did you manage to solve this problem? I did as @milos-makes-maps suggested but no change. And my first map shows a single blue color (terra::plotRGB(land_cover_bosnia).
@@Ismael-WASH I dont know if youtube or milos directly remove my replies , but I've commented 3 times before and they all appear to be gone. I have a solution on a pastebin link but I do not know if it is going to be taken down again
Thank you for the fantastic presentation! I only have one issue when I get to #10 Render Object: whenever I run rayshader::render_highquality code, I get this error: "Error in path_material_df$type : $ operator is invalid for atomic vectors" I tried running it in RStudio and VS Code, but I get the same error. Any suggestions on how to fix it? Thanks in advance!
I am truly in love with this map style! I'd love to use it to show sample locations for a project in Ecuador, though I wondered whether you have some aesthetic tips for showing which part of the map borders the ocean and which part borders a land mass? By cutting the country out entirely this differentiation is lost
Hi, and thanks for your kind words! If I'm not mistaken, you would like to show the surrounding water , too. For this purpose I would carve out a bounding box based on your country polygon. If you're not sure how to do it, check out one of my previous tutorials here: th-cam.com/video/kGadI6_ZIR4/w-d-xo.html
@@milos-makes-maps Another question in this regard: would it also be possible to create similar maps using Python actually? I can imagine you hold a strong preference for R, but aren't there packages in Python with similar functionalities?
@@milos-makes-maps Hi Milos, I have tried running the code a few times, and even after applying some modifications I was not able to render the map. Since you are generating these maps on a rather frequent rate, could you include one for Ecuador (like this one that shows climate details)?
Thank you for this lovely video.. My requirement is same but I want to use the existing LULC file of my area of interest, not the living atlas or anything like that... Any help will be highly appreciated
Excellent contribution. One question, would it be possible to extract percentage data for each characteristic? That is, percentage values for Water, Trees, Built area, etc. Thks
Thanks! Yes, it's definitely possible in R and I have a tutorial that shows you how to calculate zonal statistics for multiple classes here: th-cam.com/video/3cSDYKQg1xw/w-d-xo.html
i have tried 2 of your scripts to run on my computer, non of them works, i think you should also tell us under which version of Rstudio and packages you create this scripts
Hi Milos, please guide me i am getting an error which is like this Error in eval(ei, envir) : object 'terra' not found, i am stuck on this step please guide how to resolve this issue?
I finally managed to get the final result! Although my Bosnia borders are sharp and not detailed like yours. I realized it was the resolution. Only 20 works for me, 10, 03 and 01 says its corrupt. Is it because my pc can't handle such high resolution?
Hi Milos! Thank you for your tutorial, it's wonderful. I would like to ask, is there a way to restrict the borders down to certain states or provinces only and not the entire country? Thanks so much
Thank you, and yes, you can! You can use geodata package to filter administrative levels 1, 2 or 3 for any country in the world. Here is a simple line of code that will filter Lombardy, a region in Italy: libs
i'm just finiching this masterpeace thank you so much , my question is why you dont use rstudio as an editor instead vcstudio what is the différence? keep going.
Hi, that's a smart question. RStudio and VCStudio are both editors that can help you write and run code, but they have some differences. RStudio is mainly designed for the R programming language, which is widely used for data analysis and visualization. VCStudio, or Visual Studio Code, is a more general-purpose editor that supports many programming languages, including R. VCStudio also has more features and extensions that can enhance your coding experience. You can use either editor depending on your preference and needs. I hope this helps. Thank you for your support!
That's a great question but I haven't tried it yet. It should be possible via ground_material = rayrender::diffuse_material(), which is passed to render_highquality
Hey Milos! I have been trying this for a whole day but the datas i am downloading from Esri Living atlas, are showing to be single banded. Why this? I have tried like 100 times! But it says the same. Thats why didnt get any result at last! It was just a plain blank map! Please help me.
Can i know which code editor has been used here? I am just initially going to start learn R. So it would be helpful if i get to know which code editor or studio will be good for a beginner.
Hi Dr. Thank you for the interesting videos. Is it possible to replace giscoR with shapefiles, since I would like to work on a province or District in Zimbabwe.
Thanks! You can use geodata package to retrieve subnational boundaries for any country in the world. For example, this code will fetch the administrative level 1 shapefile of Zimbabwe in sf format. You can then filter NAME_1 column to get the boundaries for your region of interest zimbabwe_admin1 sf::st_as_sf()
Hey, sorry to hear that! What does it mean 3d view didn't create? The window didn't open at all or it closed when you tried to render the scene or something else?
Hey! Love your channel and tutorials. I have tried to create this one following along exactly, but my map comes out as all blue for water. My terra version is 1.7-78 and rayshader is 0.35.7. I'm on MacOS that makes any difference. Any ideas?
Same here too! I tried to figure out for weeks using different platforms from Windows to MacOS, Rstudio and Kaggle. Still trying to figure it out. Any follow ups?
truly big fan of your works and amazed by the way you present them. I am an architecture student from India and have zero knowledge of coding and I can’t afford to take a course to learn it, neither will I be able to learn it quick. I need similar maps for my project presentation. Do I simply copy the codes you put and change the city name that I want to generate a 3D map for? How do I go about it
Hi nice to meet you ☺️! You would need to create a bounding box, turn it into a spatial object, and pass it into my code in place of the country polygon. You can check how to create a bounding box in my tutorial here: th-cam.com/video/qTDf5VVnjMM/w-d-xo.htmlsi=AOFmoJc9poaZqeNA
Hello, Thank you for the video. Got this error on step 4: Error in h(simpleError(msg, call)) : error in evaluating the argument 'x' in selecting a method for function 'project': error in evaluating the argument 'x' in selecting a method for function 'aggregate': [crop] file exists. You can use 'overwrite=TRUE' to overwrite it Used same code as you have.
@@milos-makes-maps windows, fresh R installation. Also - absolute beginner in R, so this is all a bit overwhelming atm. Tried retyping your code, and it had error from the get go with, while your copied code hadn't had that error. If you have any good book or YT series to recommend to start with, would appreciate the advice.
Thank you so much for your kind words! I'm thrilled to hear that you found the tutorial helpful. If you're experiencing issues with packages unloading, it might be related to your system's memory or the specific version of R you're using. Don't worry, though; there's no need to throw anything away. You can always reinstall the packages or try updating them to the latest version. Keep exploring and experimenting, and you'll master 3D land cover mapping in no time! If you have any questions or need further assistance, feel free to reach out.
I want to use a shp with a custom projection (WGS84 UTM Zone 19S) I assume i have to change something in #4 load tiles because thats where i specify epsg:4326. What do I have to add or change to make it work? Because currently I get a land_cover_map full black. Thanks!
That's amazing. I have an error when i run the code. It says "Error in st_transform.sfc(st_geometry(x), crs, ...) : cannot transform sfc object with missing crs".
This error means that you are trying to transform a spatial feature object in R, but the object does not have a coordinate reference system (CRS) defined. Where does this error occur?
Thank you so much Dr. This error occurs when I try to Load my downloaded tiles. I have already defined the CRS. But still this error occurs.@@milos-makes-maps
what I am facing is merging 4 tiffs together. I tried all possible solution for this. Another issue is on coloring the plot. The plot getting only one color which is blue. I am not getting expected result. My query is why it getting one single color on the map plotting
It sounds like you’re encountering a couple of challenges with your 3D land cover map. Merging TIFF files can be tricky, and it’s essential to ensure they’re aligned and projected correctly before merging. For the coloring issue, you might be using the latest rayshader version and users have already reported the issue to the package owner here: github.com/tylermorganwall/rayshader/issues/308. Alternatively, you can install an earlier version of rayshader that I used in my tutorial by running this chunk in R: remotes::install_version( "rayshader", version = "0.35.7", repos = "cran.us.r-project.org" )
I ran into the same problem, half of my map of France is blue. I think it is the coordinate system that is wrong, as I used the same as you did for Bosnia. Could that be the issue?
Hi, thanks for the video ! I have one question, is it normal that in my elmat matrix i have a lot of NaN values. Then when i run the render_highquality function i have this error message : Error in rayvertex::validate_mesh(mesh) : !any(is.na(normals)) is not TRUE. I tried the exact same code as you. If someone can help me
Hi, you are receiving this error because of some issues with the current version of rayshader. You should install an earlier version using this chunk: require(remotes) install_version("rayshader", version = "0.35", repos = "cran.us.r-project.org")
Hi, that's a relevant and interesting question. GiscoR doesn't use India's official boundaries. But you can check one of my tutorials that shows you how to fetch official borders of India. Here is the video: th-cam.com/video/hEcKZAS7gK4/w-d-xo.html
Hi! amazing tutorial man. I have this Issue with the colors when i plot the resampled version. Im doing a region in Chile so I reprojected the maps that might be it? Land_cover_RM (before resampled) gives me the colors just fine. Thanks man! > terra::plotRGB(land_cover_RM_resampled) Error in grDevices::rgb(RGB[, 1], RGB[, 2], RGB[, 3], alpha = out$alpha, : invalid value of 'maxColorValue'
Hey! The error might occur because there is an invalid value of 'maxColorValue' in the `rgb()` function from the `grDevices` package in R. It might be due to incorrect or unexpected input values when plotting the resampled version of your land cover data. According to the R documentation, the `rgb()` function expects values between 0 and 1 for the red, green, and blue components of the color. However, you are passing values between 1 and 255, which are out of range. This causes the `rgb()` function to return an error message. To fix this error, you can either divide your values by 255 to get them in the 0-1 range, or use the `maxColorValue` argument to specify the maximum value of your color components. I hope this helps you resolve your error. If you have any other questions, feel free to ask me. 😊
@@milos-makes-maps Hi! Thank you for sharing! I have the exact same problem. What is the best way to normalize our values? Not sure if is a problem on resampling, but all values on RGB are NaN after resampling. Thank you!
Hi Dr, I'm actually doing it for one of the countries, and I got this message saying: ''Your request will download approximately 17438Mb. That's probably too big. If you really want to do this, set override_size_check = TRUE. Note that the OpenTopography API Limit will likely be exceeded.", will API limit affect the download because I can see it will take time to download the elevation ?
That's a bummer! Can you give me more info about the place where the error occurs? There are several ways to identify errors in R, such as using traceback (), browser (), debug (), or trace () functions. You can also use options (error = recover) and options (show.error.locations = TRUE) to launch a debugging session and locate the source of the error
@@milos-makes-maps I just started learning R, i come from ArcGIS Pro so im having way more difficulties than I thought. The error keeps happening on Country Borders. I tried what u said but maybe Im doing it wrong. I keep getting the same error.
@@lucas6419 I appreciate your patience and the willingness to learn R. Sometimes it can be tedious though 😁. It looks like you're having an issue with giscoR. Can you add these two lines to that function: cache = TRUE, update_cache = TRUE
I tried doing Austria. Changed the Country Code, the land cover data and names but I keep getting this error in 8. Digital Elevation Model. Error in h(simpleError(msg, call)) : error in evaluating the argument 'x' in selecting a method for function 'project': error in evaluating the argument 'x' in selecting a method for function 'aggregate': [crop] extents do not overlap
@@milos-makes-maps Maybe the issue is that there's 4-5 .tif files. I tried in aggregate changing the fact number but I don't know what else to do. Thanks!
Thank you Milos for another amazing tutorial. It exceeded all my expectations.
My pleasure, Luis! Thanks for watching it!
Thank you for this great tutorial and providing all your knowledge for free! ☺
Another banger video! 💯🔥
Thank you for your generous support, Aleks!
This is great! Could you please share your PC's specifications? I have tried several times on my computer, but the final output is not satisfactory.
This is very insightful and i will apply that to my own understanding of R and GIS data
WOW, that was amazing! Thanks!
hi, I'm doing a land cover map for Serbia and I'm stuck at 15:34 mark, where my min and max values are 1.000000 and 1.000008 respectively. I've been using chatgpt to troubleshoot my issue for 2 hours now but with no luck at all. what could be the issue?
Hi, try updating the terra package to the latest version, then restart your R session and try again. This has already helped a few users
@@milos-makes-maps Thank you!
Hi @@KatarinaaBlaa
Did you manage to solve this problem? I did as @milos-makes-maps suggested but no change. And my first map shows a single blue color (terra::plotRGB(land_cover_bosnia).
@@Ismael-WASH I dont know if youtube or milos directly remove my replies , but I've commented 3 times before and they all appear to be gone. I have a solution on a pastebin link but I do not know if it is going to be taken down again
Thank you for the fantastic presentation!
I only have one issue when I get to #10 Render Object: whenever I run rayshader::render_highquality code, I get this error:
"Error in path_material_df$type : $ operator is invalid for atomic vectors"
I tried running it in RStudio and VS Code, but I get the same error. Any suggestions on how to fix it? Thanks in advance!
Thanks for this powerful presentation
I am truly in love with this map style! I'd love to use it to show sample locations for a project in Ecuador, though I wondered whether you have some aesthetic tips for showing which part of the map borders the ocean and which part borders a land mass? By cutting the country out entirely this differentiation is lost
Hi, and thanks for your kind words! If I'm not mistaken, you would like to show the surrounding water , too. For this purpose I would carve out a bounding box based on your country polygon. If you're not sure how to do it, check out one of my previous tutorials here: th-cam.com/video/kGadI6_ZIR4/w-d-xo.html
@@milos-makes-maps Another question in this regard: would it also be possible to create similar maps using Python actually? I can imagine you hold a strong preference for R, but aren't there packages in Python with similar functionalities?
@@Appelouz waiting for his reply
@@milos-makes-maps Hi Milos, I have tried running the code a few times, and even after applying some modifications I was not able to render the map. Since you are generating these maps on a rather frequent rate, could you include one for Ecuador (like this one that shows climate details)?
@@Appelouz Sure!
Thank you for this lovely video.. My requirement is same but I want to use the existing LULC file of my area of interest, not the living atlas or anything like that... Any help will be highly appreciated
Excellent contribution. One question, would it be possible to extract percentage data for each characteristic? That is, percentage values for Water, Trees, Built area, etc. Thks
Thanks! Yes, it's definitely possible in R and I have a tutorial that shows you how to calculate zonal statistics for multiple classes here: th-cam.com/video/3cSDYKQg1xw/w-d-xo.html
I know python, I have started learning R just to try these masterpieces.
i have tried 2 of your scripts to run on my computer, non of them works, i think you should also tell us under which version of Rstudio and packages you create this scripts
Hi Milo , how to get the URL for the tiles in the map? thank you for the tutorial!
Hey! You can find the data here: livingatlas.arcgis.com/landcoverexplorer/#mapCenter=122.354%2C63.854%2C11&mode=step&timeExtent=2017%2C2022&year=2022
Thanks for this nice resource. I wonder how long the rendering take in your case? for me, it was quite slow in the case of Bosnia
Object rendering can be a tedious process and it took me around 30 minutes
Is it possible to remove the "black" line that borders each color in the raster? Thanks again.
Hello, can you specify what you mean by "black lines"? And is this related to the creation of the PNG file or rendering the 3D object?
Hi Milos, please guide me i am getting an error which is like this Error in eval(ei, envir) : object 'terra' not found, i am stuck on this step please guide how to resolve this issue?
I finally managed to get the final result! Although my Bosnia borders are sharp and not detailed like yours. I realized it was the resolution. Only 20 works for me, 10, 03 and 01 says its corrupt. Is it because my pc can't handle such high resolution?
Hmmmm, if you are getting the corrupt file error try setting cache = TRUE and update_cache = TRUE within gisco_get_countries()
Hi Milos! Thank you for your tutorial, it's wonderful.
I would like to ask, is there a way to restrict the borders down to certain states or provinces only and not the entire country? Thanks so much
Thank you, and yes, you can! You can use geodata package to filter administrative levels 1, 2 or 3 for any country in the world. Here is a simple line of code that will filter Lombardy, a region in Italy:
libs
@@milos-makes-maps Thank you! I'll try this
i'm just finiching this masterpeace thank you so much , my question is why you dont use rstudio as an editor instead vcstudio what is the différence? keep going.
Hi, that's a smart question. RStudio and VCStudio are both editors that can help you write and run code, but they have some differences. RStudio is mainly designed for the R programming language, which is widely used for data analysis and visualization. VCStudio, or Visual Studio Code, is a more general-purpose editor that supports many programming languages, including R. VCStudio also has more features and extensions that can enhance your coding experience. You can use either editor depending on your preference and needs. I hope this helps. Thank you for your support!
Great video, but Im having a problem where it says command not found when I run pattern and full.names
Thank you for the video. Do you know if you can make a background a gradient colour with shadows mixed smoothly? In rayshader
That's a great question but I haven't tried it yet. It should be possible via ground_material = rayrender::diffuse_material(), which is passed to render_highquality
@@milos-makes-maps Thank you. I will let you know if I try one
Hey Milos! I have been trying this for a whole day but the datas i am downloading from Esri Living atlas, are showing to be single banded. Why this? I have tried like 100 times! But it says the same. Thats why didnt get any result at last! It was just a plain blank map! Please help me.
Can i know which code editor has been used here? I am just initially going to start learn R. So it would be helpful if i get to know which code editor or studio will be good for a beginner.
Hey, I use R in VSCode and you can find the instructions on how to set it up here: th-cam.com/video/k79H8EeR5Jo/w-d-xo.html
Hi Dr. Thank you for the interesting videos. Is it possible to replace giscoR with shapefiles, since I would like to work on a province or District in Zimbabwe.
Thanks! You can use geodata package to retrieve subnational boundaries for any country in the world. For example, this code will fetch the administrative level 1 shapefile of Zimbabwe in sf format. You can then filter NAME_1 column to get the boundaries for your region of interest
zimbabwe_admin1
sf::st_as_sf()
milos amazing work cong. but i use your code but 3d wiew didnt created. what is the system reqiument vga card or system ?
Hey, sorry to hear that! What does it mean 3d view didn't create? The window didn't open at all or it closed when you tried to render the scene or something else?
Hey Milos, how come your R looks different then mine? Mine is white and looks kind of old/out dated. Are you in R Studio with a blackout filter?
Hi George, I use VSCode and I prepared a tutorial that shows you how to set up R in it: th-cam.com/video/k79H8EeR5Jo/w-d-xo.htmlsi=ef0Uk7WscxXQCtaW
Hey! Love your channel and tutorials. I have tried to create this one following along exactly, but my map comes out as all blue for water. My terra version is 1.7-78 and rayshader is 0.35.7. I'm on MacOS that makes any difference. Any ideas?
Same here too! I tried to figure out for weeks using different platforms from Windows to MacOS, Rstudio and Kaggle. Still trying to figure it out. Any follow ups?
truly big fan of your works and amazed by the way you present them. I am an architecture student from India and have zero knowledge of coding and I can’t afford to take a course to learn it, neither will I be able to learn it quick. I need similar maps for my project presentation. Do I simply copy the codes you put and change the city name that I want to generate a 3D map for? How do I go about it
Hi nice to meet you ☺️! You would need to create a bounding box, turn it into a spatial object, and pass it into my code in place of the country polygon. You can check how to create a bounding box in my tutorial here: th-cam.com/video/qTDf5VVnjMM/w-d-xo.htmlsi=AOFmoJc9poaZqeNA
thank you
Hello,
Thank you for the video.
Got this error on step 4:
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'x' in selecting a method for function 'project': error in evaluating the argument 'x' in selecting a method for function 'aggregate': [crop] file exists. You can use 'overwrite=TRUE' to overwrite it
Used same code as you have.
Hi, are you running this in Mac? I would advise you to update terra. This error is associated with terra
@@milos-makes-maps windows, fresh R installation.
Also - absolute beginner in R, so this is all a bit overwhelming atm. Tried retyping your code, and it had error from the get go with, while your copied code hadn't had that error.
If you have any good book or YT series to recommend to start with, would appreciate the advice.
Hello. Can we use USGS dataset tile instead of ESRI??
Good question, you could theoretically use it but you would need to adjust for the different composition of the land classes in the USGS data
Wow, that was amazing! Thanks for the videos. I have gone through this tutorial, and some packages are unloading. can I learn this from the trash?
Thank you so much for your kind words! I'm thrilled to hear that you found the tutorial helpful. If you're experiencing issues with packages unloading, it might be related to your system's memory or the specific version of R you're using. Don't worry, though; there's no need to throw anything away. You can always reinstall the packages or try updating them to the latest version. Keep exploring and experimenting, and you'll master 3D land cover mapping in no time! If you have any questions or need further assistance, feel free to reach out.
I want to use a shp with a custom projection (WGS84 UTM Zone 19S) I assume i have to change something in #4 load tiles because thats where i specify epsg:4326. What do I have to add or change to make it work? Because currently I get a land_cover_map full black.
Thanks!
Hi Lucas 👋🏻, you will need to change line 67 in main.r in the repo from crs
Can you add coordinates to your maps?
Hey Milos can i do this for state level or district level.
Hey Gaurav, absolutely, you would just need to load your custom shapefile using st_read from sf package
@@milos-makes-maps ok Thank you
That's amazing. I have an error when i run the code. It says "Error in st_transform.sfc(st_geometry(x), crs, ...) :
cannot transform sfc object with missing crs".
This error means that you are trying to transform a spatial feature object in R, but the object does not have a coordinate reference system (CRS) defined. Where does this error occur?
Thank you so much Dr.
This error occurs when I try to Load my downloaded tiles. I have already defined the CRS. But still this error occurs.@@milos-makes-maps
Dear Milos, I got " Error: unexpected symbol in:
" #-----------------
rayshader" " errors what will be the reason?
Hey, it looks like you may unintentionally ran the comment part: #----------------
Can you rerun the code again without that part?
what I am facing is merging 4 tiffs together. I tried all possible solution for this. Another issue is on coloring the plot. The plot getting only one color which is blue. I am not getting expected result. My query is why it getting one single color on the map plotting
It sounds like you’re encountering a couple of challenges with your 3D land cover map. Merging TIFF files can be tricky, and it’s essential to ensure they’re aligned and projected correctly before merging. For the coloring issue, you might be using the latest rayshader version and users have already reported the issue to the package owner here: github.com/tylermorganwall/rayshader/issues/308. Alternatively, you can install an earlier version of rayshader that I used in my tutorial by running this chunk in R:
remotes::install_version(
"rayshader",
version = "0.35.7",
repos = "cran.us.r-project.org"
)
Did it work?
did you find any solution ?
I ran into the same problem, half of my map of France is blue. I think it is the coordinate system that is wrong, as I used the same as you did for Bosnia. Could that be the issue?
@@geologic987 no, it might be because if the terra version you are using. Please install the latest version of terra and restart your R session
Hi, thanks for the video !
I have one question, is it normal that in my elmat matrix i have a lot of NaN values. Then when i run the render_highquality function i have this error message : Error in rayvertex::validate_mesh(mesh) :
!any(is.na(normals)) is not TRUE.
I tried the exact same code as you.
If someone can help me
Hi, you are receiving this error because of some issues with the current version of rayshader. You should install an earlier version using this chunk:
require(remotes)
install_version("rayshader", version = "0.35", repos = "cran.us.r-project.org")
Hi Milos, I want to create a map for a state within India. Is it possible to do that via giscoR?
Hi, that's a relevant and interesting question. GiscoR doesn't use India's official boundaries. But you can check one of my tutorials that shows you how to fetch official borders of India. Here is the video: th-cam.com/video/hEcKZAS7gK4/w-d-xo.html
hello
i have a province shp
can i crop it based on my province shp instead using country?
Yes
@@kocheng-pd4bo country_sf
Is there any tutorial for this excercise with Python is available?
While you are here, feel free to try out the R solution and let me know how it goes
Hi! amazing tutorial man. I have this Issue with the colors when i plot the resampled version. Im doing a region in Chile so I reprojected the maps that might be it? Land_cover_RM (before resampled) gives me the colors just fine. Thanks man!
> terra::plotRGB(land_cover_RM_resampled)
Error in grDevices::rgb(RGB[, 1], RGB[, 2], RGB[, 3], alpha = out$alpha, :
invalid value of 'maxColorValue'
Hey! The error might occur because there is an invalid value of 'maxColorValue' in the `rgb()` function from the `grDevices` package in R. It might be due to incorrect or unexpected input values when plotting the resampled version of your land cover data.
According to the R documentation, the `rgb()` function expects values between 0 and 1 for the red, green, and blue components of the color. However, you are passing values between 1 and 255, which are out of range. This causes the `rgb()` function to return an error message.
To fix this error, you can either divide your values by 255 to get them in the 0-1 range, or use the `maxColorValue` argument to specify the maximum value of your color components.
I hope this helps you resolve your error. If you have any other questions, feel free to ask me. 😊
@@milos-makes-maps Hi! Thank you for sharing!
I have the exact same problem. What is the best way to normalize our values?
Not sure if is a problem on resampling, but all values on RGB are NaN after resampling.
Thank you!
@@SofiaOliveira-os9xi Now I'm curious why this is happening 🤔. Can you share the link to your code in a GitHub repo?
Did you find a solution ?@@SofiaOliveira-os9xi
Any solution ?
Hi Dr, I'm actually doing it for one of the countries, and I got this message saying: ''Your request will download approximately 17438Mb. That's probably too big. If you
really want to do this, set override_size_check = TRUE. Note
that the OpenTopography API Limit will likely be exceeded.", will API limit affect the download because I can see it will take time to download the elevation ?
Hi Sofia, it looks like you are trying to download a huge DEM file vis elevate? You should try to decrease the zoom so that the output is smaller
@@milos-makes-maps Thanks Dr, I will do as you said, and give it a try!
> terra::plotRGB(land_cover_bosnia_resampled)
Error in grDevices::rgb(RGB[, 1], RGB[, 2], RGB[, 3], alpha = out$alpha, :
invalid value of 'maxColorValue'
I keep getting this error:
Skipping unit = BA
Corrupted file
Error: No results for c('BA')
D:!
That's a bummer! Can you give me more info about the place where the error occurs? There are several ways to identify errors in R, such as using traceback (), browser (), debug (), or trace () functions. You can also use options (error = recover) and options (show.error.locations = TRUE) to launch a debugging session and locate the source of the error
@@milos-makes-maps I just started learning R, i come from ArcGIS Pro so im having way more difficulties than I thought. The error keeps happening on Country Borders. I tried what u said but maybe Im doing it wrong. I keep getting the same error.
@@lucas6419 I appreciate your patience and the willingness to learn R. Sometimes it can be tedious though 😁. It looks like you're having an issue with giscoR. Can you add these two lines to that function: cache = TRUE, update_cache = TRUE
@@milos-makes-maps when remove *resolution = "1"* then the code run successfully
@@surajbhagat2672 Hm, that's interesting 🤔. Thanks for sharing 🙏🏼
I am trying to run the R program but it's getting stuck ...after 9th section
I'm sorry to hear that. Can you share the error message and the offending line
I tried doing Austria. Changed the Country Code, the land cover data and names but I keep getting this error in 8. Digital Elevation Model.
Error in h(simpleError(msg, call)) :
error in evaluating the argument 'x' in selecting a method for function 'project': error in evaluating the argument 'x' in selecting a method for function 'aggregate': [crop] extents do not overlap
Did you change the land cover tile URLs? Because they are not the same for Bosnia and Austria
@@milos-makes-maps I did Indeed.
@@lGaaspouhMc can you tell me the names of the land cover tiles that you downloaded?
@@milos-makes-maps lc2022/47N_20220101-20230101.tif, 48P, 47P, 48O, 47Q. Im trying with Thailand now btw. I get the same error.
@@milos-makes-maps Maybe the issue is that there's 4-5 .tif files. I tried in aggregate changing the fact number but I don't know what else to do.
Thanks!
Sir, how can i use my own shp file in this project??
Also how can I use my own LULC tif imagery??
You can use st_read from sf package to load your shapefile like this:
country_sf
@@milos-makes-maps thanks a lot sir, and how to use any other LULC tif image instead of ESRI live atlas
which editor is this ?
VSCode
👏👏👏
♥️
What is the best way to learn Python at this level?
R is great!
Sorry, I never used Python at this level
Can I have codes : TIA
ESA world cover data is more detailed than ESRI .
Thanks! Is there a study that compares the two datasets? I'm keen on making a more detailed tutorial that uses ESA dataset
@@milos-makes-maps this video compares the datasets : th-cam.com/video/6CVJEtl3CZo/w-d-xo.htmlsi=9D3CqJXet-RaSbN9
crs_lambert
Hello, you can find CRS for any country here: epsg.io/