Skip to content

Instantly share code, notes, and snippets.

View dbuscombe-usgs's full-sized avatar
🎯
Focusing

Daniel Buscombe dbuscombe-usgs

🎯
Focusing
186 followers · 109 following
View GitHub Profile
@dbuscombe-usgs
dbuscombe-usgs / EndPoint_bearing_Distance.m
Last active December 20, 2015 06:09
Destination point given distance and bearing from start point You have starting points: xcoordinate, a list of x-coordinates (longitudes) ycoordinate, a list of y-coordinates (latitudes) num_samples, the number of samples in the plane towards the destination point bearings: heading, a list of headings (degrees) and distances: range, a list of di…
R = 6371 % mean radius of the Earth, in km
% pre-allocate for results
x=cell(1,num_points)
y=cell(1,num_points)
for k=1:num_points
% COORDINATES OF STARTING POINT
lon1=xcoordinate(k) ;
lat1=ycoordinate(k) ;
@dbuscombe-usgs
dbuscombe-usgs / bash memory aid.sh
Last active December 20, 2015 06:09
Some common things I do in BASH
#Do something if a file exists
if [ -f $target_dir"/$target_file" ]
then do something here
fi
#Calculate pi using bc
pi=`echo "4*a(1)" | bc -l`
#Convert degs to radians
rad=`echo "$deg*($pi/180)" | bc -l`
@dbuscombe-usgs
dbuscombe-usgs / cow flip.sh
Last active March 30, 2021 16:55
How to start a BASH script (flip a coin and use cowsay)
# flip a coin
FLIP=$(($(($RANDOM%10))%2))
# if heads, use cowsay to present your splash
if [ $FLIP -eq 1 ]
then
# start with some wisdom
fortune -s | cowsay -n | zenity --text-info --title="Your message here" --width 500 --height 500
else
@dbuscombe-usgs
dbuscombe-usgs / cs2cs 1.sh
Last active December 20, 2015 06:19
Convert coordinates in Arizona Central State Plane to WGS84 Latitude/Longitude using cs2cs
cs2cs -f "%.6f" +proj=tmerc +lat_0=31 +lon_0=-111.9166666666667 +k=0.9999 +x_0=213360 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs infile > outfile
@dbuscombe-usgs
dbuscombe-usgs / awk trim.awk
Last active December 20, 2015 06:19
Trim file into 2 columns separated by a tab
awk '{print $(NF-2),"\t",$(NF-1)}' outfile > outfile2
@dbuscombe-usgs
dbuscombe-usgs / zernicke.py
Last active December 20, 2015 06:19
Classify images based on Zernicke moments
features = np.zeros([len(imlist),25])
for i,f in enumerate(imlist):
im = np.array(Image.open(f).convert('L'))
h=mahotas.features.zernike_moments(im, 50, 8)
features[i] = h.flatten()
@dbuscombe-usgs
dbuscombe-usgs / haralick.py
Last active December 20, 2015 06:19
Classify images based on Haralick descriptors
features = np.zeros([len(imlist),13])
for i,f in enumerate(imlist):
im = np.array(Image.open(f))
h = mahotas.texture.haralick(im)[0]
features[i] = h.flatten()
@dbuscombe-usgs
dbuscombe-usgs / class moments.py
Last active December 20, 2015 06:19
Classify images based on moments
features = np.zeros([len(imlist),5])
for i,f in enumerate(imlist):
im = np.array(Image.open(f).convert('L'))
Z=np.array([],np.float)
for k in range(5):
Z=np.append(Z,mahotas.moments(im, k, k))
features[i]=Z
@dbuscombe-usgs
dbuscombe-usgs / tas class.py
Last active December 20, 2015 06:19
Classify images based on TAS features
features = np.zeros([len(imlist),162])
for i,f in enumerate(imlist):
im = np.array(Image.open(f))
h=mahotas.features.tas(im)
features[i] = h.flatten()
@dbuscombe-usgs
dbuscombe-usgs / pftas class.py
Last active December 20, 2015 06:19
Classify images based on PFTAS features
features = np.zeros([len(imlist),162])
for i,f in enumerate(imlist):
im = np.array(Image.open(f))
h=mahotas.features.pftas(im)
features[i] = h.flatten()