Register
Login
Resources
Docs Blog Datasets Glossary Case Studies Tutorials & Webinars
Product
Data Engine LLMs Platform Enterprise
Pricing Explore
Connect to our Discord channel
Resources  
Docs Blog Datasets Glossary Case Studies Tutorials & Webinars
Explore Pricing
Product  
Data Engine LLMs Platform Enterprise
Connect to Dagshub Discord channel
Register
Login

maxrjones
/
gmt-fork
mirror of https://github.com/meghanrjones/gmt.git

Watch 3
Issues   Discussions 0

hexbinning.sh 2.0 KB
Permalink History Raw

You have to be logged in to leave a comment. Sign In
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46 
  1. #!/usr/bin/env bash
    2. 

  2. # Test hexagon tile binning on counting values in gmtbinstats
    3. 

  3. #
    4. 

  4. # Hardwire some values to draw the bottom plot so we avoid lots of shell math
    5. 

  5. # The commented stuff shows how the script and data were developed
    6. 

  6. #e=5
    7. 

  7. #n=3
    8. 

  8. #dy=1
    9. 

  9. #ny=$(gmt math -Q $n SUB $dy DIV =)
    10. 

  10. #r=$(gmt math -Q $dy 0.5 MUL 30 COSD DIV =)
    11. 

  11. #dx=$(gmt math -Q 3 $r MUL =)
    12. 

  12. #nx=$(gmt math -Q $e $dx DIV RINT 1 ADD =)
    13. 

  13. #e2=$(gmt math -Q 0 $nx 1 SUB $dx MUL ADD =)
    14. 

  14. #r=$(gmt math -Q $dy 30 COSD DIV =)
    15. 

  15. #echo "Revised: 0/$e2/0/$n -I$dx/$dy for a pseudo $nx by $ny grid with hexagon radius $r"
    16. 

  16. #dx2=$(gmt math -Q $dx 2 DIV =)
    17. 

  17. #dy2=$(gmt math -Q $dy 2 DIV =)
    18. 

  18. #wb=$(gmt math -Q 0 $dx2 ADD =)
    19. 

  19. #eb=$(gmt math -Q $e2 $dx2 SUB =)
    20. 

  20. #sb=$(gmt math -Q 0 $dy2 ADD =)
    21. 

  21. #nb=$(gmt math -Q $n $dy2 SUB =)
    22. 

  22. #rh=$(gmt math -Q 4 30 COSD DIV =)
    23. 

  23. #gmt math -T1/100/1 0 $e2 RAND = x
    24. 

  24. #gmt math -T1/100/1 0 $n RAND = y
    25. 

  25. #gmt math -T1/100/1 50 25 NRAND = z
    26. 

  26. #gmt convert -A x y z -o1,3,5 > hex_data.txt
    27. 

  27. 

  28. # Create grids to find nodes of hexagons
    29. 

  29. gmt grdmath -R0/5.19615242271/0/3 -I1.73205080757/1 0 = b.grd
    30. 

  30. gmt grdmath -R0.866025403785/4.33012701893/0.5/2.5 -I1.73205080757/1 0 = r.grd
    31. 

  31. gmt grdmath -R0/5.19615242271/0/3 -I0.866025403785/0.5 0 = t.grd
    32. 

  32. gmt begin hexbinning ps
    33. 

  33. 	gmt grd2xyz b.grd | gmt plot -Sh3.46410161514c -Glightblue@50 -W1p -N -R0/5.19615242271/0/3 -Jx3c -X3.2c
    34. 

  34. 	gmt grd2xyz b.grd | gmt plot -Sh3.46410161514c -Glightblue -W1p
    35. 

  35. 	gmt grd2xyz r.grd | gmt plot -Sh3.46410161514c -Glightred@50 -W1p -N
    36. 

  36. 	gmt grd2xyz r.grd | gmt plot -Sh3.46410161514c -Glightred -W1p
    37. 

  37. 	gmt grdmath -R0/5.19615242271/0/3 -I0.866025403785/0.5 0 = t.grd
    38. 

  38. 	gmt grd2xyz t.grd | gmt plot -Sc0.2c -Gwhite -Wfaint -N
    39. 

  39. 	gmt grd2xyz r.grd b.grd | gmt plot -Sc0.2c -Gblack -N
    40. 

  40. 	gmt basemap -Bxafg0.866025403785 -Byafg0.5
    41. 

  41. 	gmt plot hex_data.txt -Ss0.2c -Gyellow -Wfaint
    42. 

  42. 	gmt binstats hex_data.txt -R0/5/0/3 -I1 -Th -Cn > bin.txt
    43. 

  43. 	gmt makecpt -Cjet -T1/10/1 -A50
    44. 

  44. 	gmt plot -Baf -Y12c bin.txt -C -Sh3.46410161514c -W1p -B+t"Hexagonal Binning"
    45. 

  45. 	gmt text bin.txt -F+f14p,Helvetica-Bold -Gwhite -W0.25p -N
    46. 

  46. gmt end show
    47.
Prev
Next
Tip!

Press p or to see the previous file or, n or to see the next file

Comments

Write Preview
Loading...
Integrate AWS S3
Use S3 remote
Select bucket
Access key
Finish

Use AWS S3 as storage!

Browsing data directories saved to S3 is possible with DAGsHub. Let's configure your repository to easily display your data in the context of any commit!

Specify your S3 bucket

Access key (If needed)

Congratulations!

gmt-fork is now integrated with AWS S3!

Delete Storage Key

Are you sure you want to delete this access key?

No
Yes
Integrate Google Cloud Storage
Use Google Storage
Select bucket
Upload key
Finish

Use Google Cloud Storage!

Browsing data directories saved to Google Cloud Storage is possible with DAGsHub. Let's configure your repository to easily display your data in the context of any commit!

Specify your Google Storage bucket

Service Account Key

Congratulations!

gmt-fork is now integrated with Google Cloud Storage!

Delete Storage Key

Are you sure you want to delete this access key?

No
Yes
Integrate Azure Cloud Storage
Use Azure Storage
Select bucket
Set key
Finish

Use Azure Cloud Storage!

Browsing data directories saved to Azure Cloud Storage is possible with DAGsHub. Let's configure your repository to easily display your data in the context of any commit!

Specify your Azure Storage bucket

Access key (If needed)

Congratulations!

gmt-fork is now integrated with Azure Cloud Storage!

Delete Storage Key

Are you sure you want to delete this access key?

No
Yes
Integrate S3 compatible storage
Use S3 like remote
Select bucket
Access key
Finish

Use any S3 compatible storage!

Browsing data directories saved to S3 compatible storage is possible with DAGsHub. Let's configure your repository to easily display your data in the context of any commit!

Specify your S3 bucket

Bucket name cannot be the same as the repository name. Please change one of them.

Access key (If needed)

Congratulations!

gmt-fork is now integrated with your S3 compatible storage!

Delete Storage Key

Are you sure you want to delete this access key?

No
Yes