Skip to content

redplanet.Crust.dichotomy.get_coords 

get_coords() -> ndarray

Get a list of dichotomy boundary coordinates.

The origin of the dataset is not fully clear. We use the file "dichotomy_coordinates-JAH-0-360.txt" (0.046 MiB) downloaded from Wieczorek (2022). They attribute the data to "Andrews-Hanna et al. (2008)" which is ambiguous (to me atleast, I could be missing something obvious) — my best guess is Andrews-Hanna et al. (2008).

Returns:

Type Description
np.ndarray

A numpy array of shape (n, 2) where n is the number of coordinates, and the two columns are longitude (0->360) and latitude respectively.
References
  1. Andrews-Hanna, J. C., Zuber, M. T., & Banerdt, W. B. (2008). The Borealis basin and the origin of the martian crustal dichotomy. Nature, 453(7199), 1212-1215. https://doi.org/10.1038/nature07011
  2. Wieczorek, M. A. (2022). InSight Crustal Thickness Archive [Dataset]. Zenodo. https://doi.org/10.5281/zenodo.6477509 [Data from "dichotomy/" is mirrored to https://rutgers.box.com/v/redplanet-data]
Source code in src/redplanet/Crust/dichotomy.py 
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
@substitute_docstrings
def get_coords() -> np.ndarray:
    """
    Get a list of dichotomy boundary coordinates.

    The origin of the dataset is not fully clear. We use the file "dichotomy_coordinates-JAH-0-360.txt" (0.046 MiB) downloaded from {@Wieczorek2022_icta.n}. They attribute the data to *"Andrews-Hanna et al. (2008)"* which is ambiguous (to me atleast, I could be missing something obvious) — my best guess is {@dichotomy_paper.n}.

    Returns
    -------
    np.ndarray
        A numpy array of shape `(n, 2)` where `n` is the number of coordinates, and the two columns are longitude (0->360) and latitude respectively.
    """
    fpath = _get_fpath_dataset('dichotomy_coords')
    dat_dichotomy_coords = np.loadtxt(fpath)
    return dat_dichotomy_coords