I’m so happy my parents got to cross northern lights off their bucket list!
Continue reading..AuroraI already want to revisit this neat lakeside park to snorkel!
Continue reading..Esrey ParkToday, I kicked off a long weekend vacation with my parents! We decided to take a little road trip somewhere none of us had been. Apparently, pipe vesicles can be found on the bluffs just north of Bessemer. “The traps from Bessemer to North Bessemer appear to be older than any beds exposed north…Elongate, so-called […]
Continue reading..KeweenawThis is my job submission script for the raster extraction. #!/bin/bash -l #PBS -l walltime=48:00:00,nodes=1:ppn=24,mem=61gb #PBS -m abe #PBS -M myaddress@email.com module load gdal/2.0.1 module load proj module load R module load ompi module load udunits module load geos export OMP_NUM_THREADS=1 mpirun -v -hostfile $PBS_NODEFILE -np 1 R –slave CMD BATCH ./extract.R /dev/stdout This is […]
Continue reading..Parallel Raster Extractionhttps://bookdown.org/yihui/rmarkdown/word-document.html Writing a MS-Word document using R (with as little overhead as possible) Click to access knitr.pdf Tables from R into Word https://rdrr.io/github/crsh/papaja/man/apa_table.html
Continue reading..R MarkdownSee part 1 for the relevant shape file and description of what’s going on! In the last post, we identified the map codes of interest to subset the shape file. Here’s some code for working with it in R. library(rgdal) bedrock <- readOGR(dsn=”~/Arrowhead_South_Bedrock_Geology”,”Arrowhead_South_Bedrock_Geology”) basalt = subset(bedrock,map_label==”Mnd” | map_label==”Mna”) writeOGR(basalt,dsn=”.”,layer=”basalt”,driver=”ESRI Shapefile”) So now we have a […]
Continue reading..Weekend GIS Fun: Part 2The stunning geology of the north shore is hard to ignore! So, my beach walking has led me to a rock collecting hobby, in order to learn more about the brilliant colors & textures along the shoreline. I’m interested in finding pipe vesicles, which form at the base of a lava flow. It’s kind of […]
Continue reading..Weekend GIS Fun: Part 1gdal_merge.py -init -1 -o all.tif n -1 -a_nodata -1 new_FPWF.tif new_MH.tif new_FD.tif
Continue reading..GDAL merge examplehttps://dl.acm.org/citation.cfm?id=3274985 Click to access a9ba03ba65cc9e9e0a0c1b73a8359e6772ac.pdf https://www.perrygeo.com/ Click to access postgis_zonal.pdf Click to access zonalstat_tr.pdf http://www.guru-gis.net/efficient-zonal-statistics-using-r-and-gdal/ https://robbymarrotte.weebly.com/blog/faster-zonal-statistics
Continue reading..Efficient Zonal StatisticsToday, I surveyed the plots on this WMA, and also at Split Rock Creek state park. I surprised myself by making it through Smiley Barn WMA before the morning was out! This was my record-breaking day for number of sites in a morning, which calmed my nerves about not finishing within the survey period. Now, […]
Continue reading..Eden WMA