anta-database 0.3.3__tar.gz → 0.3.5__tar.gz

{anta_database-0.3.3 → anta_database-0.3.5}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: anta_database
-Version: 0.3.3
+Version: 0.3.5
 Summary: SQLite database for the AntArchitecture Community Data
 Author: Antoine Hermant
 Author-email: antoine.hermant@etik.com
@@ -48,12 +48,12 @@ Visit the [Home Page](https://antoinehermant.github.io/anta_database/intro) for
 
 The AntADatabase contains all traced and dated Internal Reflective Horizons (isochrones) that have been published at that date. All the data is organized by flight transect for each dataset. The data is encoded in HDF5 (Hierarchical Data Format) which provides convenient data structure and performant read speeds. 
 Variables per file (when exist):
-- x(point)
-- y(point)
-- distance(point)
+- PSX(point)
+- PSY(point)
+- Distance(point)
 - IRH_DEPTH(point, age)
-- IRH_DENS(point, age)
-- ICE_THCK(point, age)
+- IRH_NUM(point, age)
+- ICE_THK(point)
 - SURF_ELEV(point, age)
 - BED_ELEV(point, age)
 
@@ -66,7 +66,7 @@ This Python module provides SQL indexing for the AntADatabase, as well as quick
 - age
 - acquisition_year (year of radar acquisition)
 - region
-- basin
+- IMBIE_basin
 - var (variable)
 - flight id
 

{anta_database-0.3.3 → anta_database-0.3.5}/README.md

@@ -16,12 +16,12 @@ Visit the [Home Page](https://antoinehermant.github.io/anta_database/intro) for
 
 The AntADatabase contains all traced and dated Internal Reflective Horizons (isochrones) that have been published at that date. All the data is organized by flight transect for each dataset. The data is encoded in HDF5 (Hierarchical Data Format) which provides convenient data structure and performant read speeds. 
 Variables per file (when exist):
-- x(point)
-- y(point)
-- distance(point)
+- PSX(point)
+- PSY(point)
+- Distance(point)
 - IRH_DEPTH(point, age)
-- IRH_DENS(point, age)
-- ICE_THCK(point, age)
+- IRH_NUM(point, age)
+- ICE_THK(point)
 - SURF_ELEV(point, age)
 - BED_ELEV(point, age)
 
@@ -34,7 +34,7 @@ This Python module provides SQL indexing for the AntADatabase, as well as quick
 - age
 - acquisition_year (year of radar acquisition)
 - region
-- basin
+- IMBIE_basin
 - var (variable)
 - flight id
 

anta_database-0.3.5/ROADMAP.md

@@ -0,0 +1,12 @@
+# Roadmap
+
+Here is a todo list of features that we want to include in the AntADatabase.
+
+## Important
+- [x] Implement HDF5 format (In process, for next release)
+- [x] Filter by radar system
+- [x] Create shapefiles, e.g. for QAntarctica (In process)
+
+## To discuss
+- [ ] A netCDF version 
+

{anta_database-0.3.3 → anta_database-0.3.5}/anta_database/compiling/compile.py

@@ -294,7 +294,7 @@ class CompileDatabase:
                             header=0,
                             sep=sep,
                             # usecols=original_new_columns.columns,
-                            na_values=['-9999', '-9999.0', 'NaN', 'nan', ''],
+                            # na_values=['-9999', '-9999.0', 'NaN', 'nan', ''],
                             dtype=str
                             )
 
@@ -304,6 +304,26 @@ class CompileDatabase:
         ds = ds[ds.columns.intersection(original_new_columns.columns)]
         ds.columns = original_new_columns[ds.columns].iloc[0].values  # renaming the columns
 
+        values_to_replace = ['-9999', '-9999.0', 'NaN', 'nan', '']
+
+        if 'PSX' in ds.columns and 'PSY' in ds.columns:
+            # Create a mask for rows where both PSX and PSY are in values_to_replace
+            mask = (
+                ds['PSX'].isin(values_to_replace) &
+                ds['PSY'].isin(values_to_replace)
+            )
+            if mask.any():
+                problematic_indices = ds.index[mask].tolist()
+                warnings.warn(
+                    f"Found {len(problematic_indices)} rows where both PSX and PSY are in {values_to_replace}. "
+                    f"Indices: {problematic_indices}. Dropping these rows."
+                )
+
+            ds = ds[~mask].copy()
+
+        cols_to_replace = [col for col in ds.columns if col not in ['PSX', 'PSY']]
+        ds[cols_to_replace] = ds[cols_to_replace].replace(values_to_replace, pd.NA)
+
         pattern_values = original_new_columns[1:]
         pattern_header = original_new_columns.iloc[0]
         pattern_values.columns = pattern_header
@@ -873,9 +893,9 @@ class CompileDatabase:
 
                     df = gpd.GeoDataFrame({
                         "geometry": [Point(x, y) for x, y in zip(ds.PSX.values, ds.PSY.values)],
-                        "irh_num": N.astype('int16'),
+                        "irh_num": N.astype('int8'),
                         "max_age": age_max_depth.astype('int32'),
-                        "max_depth": max_depth.astype('int16'),
+                        "max_depth": np.rint(max_depth).astype('int16'),
                         # "flight_id": [flight_id] * len(N),
                         # "doi_data": [doi_data] * len(N),
                     }, crs="EPSG:3031")

{anta_database-0.3.3 → anta_database-0.3.5}/anta_database/database/database.py

@@ -715,8 +715,9 @@ class Database:
         self,
         metadata: Union[None, Dict, 'MetadataResult'] = None,
         data_dir: Optional[str] = None,
-        downscale_factor: Optional[str] = None,
+        downsampling_factor: Optional[str] = None,
         disable_tqdm: bool = False,
+        fraction_depth: Optional[bool] = False,
     ) -> Generator[Tuple[pd.DataFrame, Dict]]:
         """
         Generates xarray Datasets from HDF5 files, one at a time, with lazy loading.
@@ -744,35 +745,51 @@ class Database:
         file_paths = self._get_file_paths_from_metadata(metadata=md)
         file_paths = np.unique(file_paths) # Be carefull as many pointers point to the same file
 
+        if disable_tqdm or self._disable_tqdm:
+            disable_tqdm = True
+
         for file_path in tqdm(file_paths, desc='Generating dataframes', total=len(file_paths), unit='file', disable=disable_tqdm):
             full_path = os.path.join(data_dir, file_path)
             file_md = self._get_file_metadata(file_path)
             ds = h5py.File(full_path, 'r')
-            df = pd.DataFrame({'PSX': ds['PSX'][::downscale_factor],
-                                'PSY': ds['PSY'][::downscale_factor]})
+            df = pd.DataFrame({'PSX': ds['PSX'][::downsampling_factor],
+                                'PSY': ds['PSY'][::downsampling_factor]})
 
             if 'Distance' in ds.keys():
-                df['Distance'] = ds['Distance'][::downscale_factor]
+                df['Distance'] = ds['Distance'][::downsampling_factor]
 
+            var_impl = []
+            age_impl = []
             for var in md['var']:
                 if var == 'IRH_DEPTH':
+                    var_impl.append(var)
                     for age in md['age']:
                         irh_values = ds['IRH_AGE'][:]
                         irh_index = np.where(irh_values == int(age))[0]
                         if len(irh_index) == 0:
                             continue
                         irh_index = irh_index[0]
+                        age_impl.append(age)
 
-                        df[int(age)] = ds[var][::downscale_factor, irh_index]
+                        df[age] = ds[var][::downsampling_factor, irh_index]
+                        if fraction_depth:
+                            if 'ICE_THK' in ds.keys():
+                                df[age] *= 100/ds['ICE_THK'][::downsampling_factor]
+                            else:
+                                df[age] = np.nan
+                                # print(f'WARNING: flight line {file_md['flight_id']} in {file_md['dataset']} does not have ICE_THK, cannot compute fraction depth, will return NaN')
 
                 else:
                     if var in ds.keys():
-                        df[var] = ds[var][::downscale_factor]
+                        var_impl.append(var)
+                        df[var] = ds[var][::downsampling_factor]
 
+            if 0 < len(age_impl) <= 1:
+                age_impl = list(age_impl)[0]
             metadata = {
                 'dataset': file_md['dataset'],
-                'var': md['var'],
-                'age': md['age'],
+                'var': var_impl,
+                'age': age_impl,
                 'flight_id': file_md['flight_id'],
                 'institute': file_md['institute'],
                 'project': file_md['project'],
@@ -789,105 +806,6 @@ class Database:
 
             yield df, metadata
 
-    # def return_dataset(
-    #     self,
-    #     metadata: Union[None, Dict, 'MetadataResult'] = None,
-    #     data_dir: Optional[str] = None,
-    #     downscale_factor: Optional[str] = None,
-    # ) -> Generator[Tuple[pd.DataFrame, Dict]]:
-    #     """
-    #     Generates xarray Datasets from HDF5 files, one at a time, with lazy loading.
-
-    #     Args:
-    #         metadata: Metadata for filtering files.
-    #         data_dir: Directory containing the data files.
-    #         vars_to_load: List of variables to load from each file.
-
-    #     Yields:
-    #         Tuple[xr.Dataset, Dict]: A lazy-loaded xarray Dataset and its metadata.
-    #     """
-    #     # Resolve metadata
-    #     md = metadata or self._md
-    #     if not md:
-    #         print('Please provide metadata of the files you want to generate the data from. Exiting...')
-    #         return
-
-    #     # Resolve data directory
-    #     data_dir = data_dir or self._db_dir
-    #     if not data_dir:
-    #         print('No data directory provided. Exiting...')
-    #         return
-
-    #     if len(md['var']) > 1:
-    #         print(f'WARNING: you requested multiple variables ({md['var']}), this will return individual dataframes for each variable in each transect.')
-
-    #     file_paths = self._get_file_paths_from_metadata(metadata=md)
-    #     file_paths = np.unique(file_paths) # Be carefull as many pointers point to the same file
-
-    #     for file_path in file_paths:
-    #         full_path = os.path.join(data_dir, file_path)
-    #         file_md = self._get_file_metadata(file_path)
-    #         for var in md['var']:
-    #             if var == 'IRH_DEPTH':
-    #                 for age in md['age']:
-    #                     with h5py.File(full_path, 'r') as ds:
-    #                         irh_values = ds['age'][:]
-    #                         irh_index = np.where(irh_values == int(age))[0]
-    #                         if len(irh_index) == 0:
-    #                             continue
-    #                         irh_index = irh_index[0]
-    #                         df = pd.DataFrame({'x': ds['x'][::downscale_factor],
-    #                                         'y': ds['y'][::downscale_factor],
-    #                                         'distance': ds['distance'][::downscale_factor],
-    #                                         var: ds[var][::downscale_factor, irh_index]})
-
-    #                         metadata = {
-    #                             'dataset': file_md['dataset'],
-    #                             'var': var,
-    #                             'age': age,
-    #                             'flight_id': file_md['flight_id'],
-    #                             'institute': file_md['institute'],
-    #                             'project': file_md['project'],
-    #                             'acquisition_year': file_md['acquisition_year'],
-    #                             'age_unc': file_md['age'],
-    #                             'reference': file_md['reference'],
-    #                             'DOI_dataset': file_md['DOI_dataset'],
-    #                             'DOI_publication': file_md['DOI_publication'],
-    #                             'flight_id': file_md['flight_id'],
-    #                             'region': file_md['region'],
-    #                             'IMBIE_basin': file_md['IMBIE_basin'],
-    #                             'radar_instrument': file_md['radar_instrument'],
-    #                         }
-
-    #                         yield df, metadata
-
-    #                 else:
-    #                     ds = h5py.File(full_path, 'r')
-    #                     df = pd.DataFrame({'PSX': ds['PSX'][::downscale_factor],
-    #                                     'PSY': ds['PSY'][::downscale_factor],
-    #                                     'Distance': ds['Distance'][::downscale_factor],
-    #                                     var: ds[var][::downscale_factor]})
-
-    #                     metadata = {
-    #                         'dataset': file_md['dataset'],
-    #                         'var': var,
-    #                         'age': None,
-    #                         'flight_id': file_md['flight_id'],
-    #                         'institute': file_md['institute'],
-    #                         'project': file_md['project'],
-    #                         'acquisition_year': file_md['acquisition_year'],
-    #                         'age_unc': file_md['age'],
-    #                         'reference': file_md['reference'],
-    #                         'DOI_dataset': file_md['DOI_dataset'],
-    #                         'DOI_publication': file_md['DOI_publication'],
-    #                         'flight_id': file_md['flight_id'],
-    #                         'region': file_md['region'],
-    #                         'IMBIE_basin': file_md['IMBIE_basin'],
-    #                         'radar_instrument': file_md['radar_instrument'],
-    #                     }
-
-    #                     yield df, metadata
-
     @property
     def plot(self):
         if self._plotting is None:

{anta_database-0.3.3 → anta_database-0.3.5}/anta_database/plotting/plotting.py

@@ -62,7 +62,7 @@ class Plotting:
     def dataset(
             self,
             metadata: Union[None, Dict, 'MetadataResult'] = None,
-            downscale_factor: Optional[int] = None,
+            downsampling_factor: Optional[int] = None,
             title: Optional[str] = None,
             xlim: Optional[tuple] = (None, None),
             ylim: Optional[tuple] = (None, None),
@@ -82,7 +82,7 @@ class Plotting:
         self._base_plot(
             color_by='dataset',
             metadata=metadata,
-            downscale_factor=downscale_factor,
+            downsampling_factor=downsampling_factor,
             marker_size=marker_size,
             title=title,
             xlim=xlim,
@@ -100,7 +100,7 @@ class Plotting:
     def institute(
             self,
             metadata: Union[None, Dict, 'MetadataResult'] = None,
-            downscale_factor: Optional[int] = None,
+            downsampling_factor: Optional[int] = None,
             title: Optional[str] = None,
             xlim: Optional[tuple] = (None, None),
             ylim: Optional[tuple] = (None, None),
@@ -120,7 +120,7 @@ class Plotting:
         self._base_plot(
             color_by='institute',
             metadata=metadata,
-            downscale_factor=downscale_factor,
+            downsampling_factor=downsampling_factor,
             marker_size=marker_size,
             title=title,
             xlim=xlim,
@@ -138,7 +138,7 @@ class Plotting:
     def flight_id(
             self,
             metadata: Union[None, Dict, 'MetadataResult'] = None,
-            downscale_factor: Optional[int] = None,
+            downsampling_factor: Optional[int] = None,
             title: Optional[str] = None,
             xlim: Optional[tuple] = (None, None),
             ylim: Optional[tuple] = (None, None),
@@ -158,7 +158,7 @@ class Plotting:
         self._base_plot(
             color_by='flight_id',
             metadata=metadata,
-            downscale_factor=downscale_factor,
+            downsampling_factor=downsampling_factor,
             scale_factor=scale_factor,
             marker_size=marker_size,
             title=title,
@@ -176,7 +176,8 @@ class Plotting:
     def var(
             self,
             metadata: Union[None, Dict, 'MetadataResult'] = None,
-            downscale_factor: Optional[int] = None,
+            downsampling_factor: Optional[int] = None,
+            fraction_depth: Optional[bool] = False,
             title: Optional[str] = None,
             xlim: Optional[tuple] = (None, None),
             ylim: Optional[tuple] = (None, None),
@@ -196,7 +197,8 @@ class Plotting:
         self._base_plot(
             color_by='var',
             metadata=metadata,
-            downscale_factor=downscale_factor,
+            downsampling_factor=downsampling_factor,
+            fraction_depth=fraction_depth,
             title=title,
             xlim=xlim,
             ylim=ylim,
@@ -215,7 +217,7 @@ class Plotting:
             self,
             metadata: Union[None, Dict, 'MetadataResult'] = None,
             elevation: Optional[bool] = False,
-            downscale_factor: Optional[int] = None,
+            downsampling_factor: Optional[int] = None,
             title: Optional[str] = None,
             xlim: Optional[tuple] = (None, None),
             ylim: Optional[tuple] = (None, None),
@@ -234,7 +236,7 @@ class Plotting:
             color_by='transect_1D',
             elevation=elevation,
             metadata=metadata,
-            downscale_factor=downscale_factor,
+            downsampling_factor=downsampling_factor,
             title=title,
             xlim=xlim,
             ylim=ylim,
@@ -261,7 +263,8 @@ class Plotting:
             self,
             metadata: Union[None, Dict, 'MetadataResult'] = None,
             elevation: Optional[bool] = False,
-            downscale_factor: Optional[int] = None,
+            fraction_depth: Optional[bool] = False,
+            downsampling_factor: Optional[int] = None,
             title: Optional[str] = None,
             xlim: Optional[tuple] = (None, None),
             ylim: Optional[tuple] = (None, None),
@@ -314,8 +317,8 @@ class Plotting:
         if color_by == 'dataset':
             if not title:
                 title = f'AntADatabase by datasets'
-            if downscale_factor == None:
-                downscale_factor = 1
+            if downsampling_factor == None:
+                downsampling_factor = 1
             datasets = list(metadata['dataset'])
             bedmap_entries = {'BEDMAP1', 'BEDMAP2', 'BEDMAP3'}
             bedmap_colors = {
@@ -341,8 +344,8 @@ class Plotting:
                 for f in file_paths:
                     full_path = os.path.join(self._db._db_dir, f)
                     with h5py.File(full_path, 'r') as ds:
-                        all_x.append(ds['PSX'][::downscale_factor])
-                        all_y.append(ds['PSY'][::downscale_factor])
+                        all_x.append(ds['PSX'][::downsampling_factor])
+                        all_y.append(ds['PSY'][::downsampling_factor])
                 df = pd.DataFrame({'PSX': np.concatenate(all_x),
                                 'PSY': np.concatenate(all_y)})
                 plt.scatter(df['PSX']/1000, df['PSY']/1000, color=colors[dataset], s=marker_size, zorder=zorder, linewidths=0)
@@ -359,8 +362,8 @@ class Plotting:
         if color_by == 'institute':
             if not title:
                 title = f'AntADatabase by institutes'
-            if downscale_factor == None:
-                downscale_factor = 1
+            if downsampling_factor == None:
+                downsampling_factor = 1
             institutes = list(metadata['institute'])
 
             color_indices = np.linspace(0.1, 0.9, len(institutes))
@@ -378,8 +381,8 @@ class Plotting:
                 for f in file_paths:
                     full_path = os.path.join(self._db._db_dir, f)
                     with h5py.File(full_path, 'r') as ds:
-                        all_x.append(ds['PSX'][::downscale_factor])
-                        all_y.append(ds['PSY'][::downscale_factor])
+                        all_x.append(ds['PSX'][::downsampling_factor])
+                        all_y.append(ds['PSY'][::downsampling_factor])
                 df = pd.DataFrame({'PSX': np.concatenate(all_x),
                                 'PSY': np.concatenate(all_y)})
                 plt.scatter(df['PSX']/1000, df['PSY']/1000, color=colors[institute], s=marker_size, linewidths=0)
@@ -402,9 +405,8 @@ class Plotting:
             else:
                 var = var[0]
 
-
             all_dfs = []
-            for ds, _ in tqdm(self._db.data_generator(metadata, downscale_factor=downscale_factor, disable_tqdm=True), desc="Plotting", total=total_traces, unit='trace', disable=self._disable_tqdm):
+            for ds, _ in tqdm(self._db.data_generator(metadata, downsampling_factor=downsampling_factor, disable_tqdm=True, fraction_depth=fraction_depth), desc="Plotting", total=total_traces, unit='trace', disable=self._disable_tqdm):
                 all_dfs.append(ds)
             df = pd.concat(all_dfs)
 
@@ -448,67 +450,61 @@ class Plotting:
                     if cmap == None:
                         cmap = cmaps.bukavu
                     extend = 'both'
-                    if not vmin: vmin = -1000
-                    if not vmax: vmax = 1000
+                    if vmin is None: vmin = -1000
+                    if vmax is None: vmax = 1000
                     if not title:
                         title = f'AntADatabase Bed Elevation'
                 elif var == 'ICE_THK':
                     if cmap == None:
                         cmap = cmaps.torch_r
                     extend = 'max'
-                    if not vmin: vmin = 0
-                    if not vmax: vmax = 4000
+                    if vmin is None: vmin = 0
+                    if vmax is None: vmax = 4000
                     if not title:
                         title = f'AntADatabase Ice Thickness'
                 elif var == 'SURF_ELEV':
                     if cmap == None:
-                        cmap = cmaps.torch_r
+                        cmap = cmaps.ice_r
                     extend = 'max'
-                    if not vmin: vmin = 2000
-                    if not vmax: vmax = 4000
+                    if vmin is None: vmin = 1000
+                    if vmax is None: vmax = 4000
                     if not title:
                         title = f'AntADatabase Surface Elevation'
                 elif var == 'BASAL_UNIT':
                     if cmap == None:
                         cmap = cmaps.torch_r
                     extend = 'both'
-                    if not vmin: vmin = 2000
-                    if not vmax: vmax = 400
+                    if vmin is None: vmin = 2000
+                    if vmax is None: vmax = 400
                     if not title:
                         title = f'AntADatabase Basal Unit'
                 scatter = plt.scatter(df['PSX']/1000, df['PSY']/1000, c=df[var], cmap=cmap, s=marker_size, vmin=vmin, vmax=vmax, linewidths=0, rasterized=True)
 
-            elif var in ['IRH_DEPTH', 'IRH_FRAC_DEPTH']:
-                if var == 'IRH_DEPTH':
-                    label = f'IRH Depth [m]'
-                    if cmap == None:
-                        cmap = cmaps.torch_r
-                    extend = 'both'
-                    age = list(metadata['age'])
-                    if not title:
+            elif var in ['IRH_DEPTH']:
+                if cmap == None:
+                    cmap = cmaps.torch_r
+                extend = 'both'
+                age = list(metadata['age'])
+                if not title:
+                    if fraction_depth:
+                        title = f'AntADatabase IRH Fraction Depth'
+                    else:
                         title = f'AntADatabase IRH Depth'
-                    if len(age) > 1:
-                        print('WARNING: Multiple layers provided: ', age,
-                            '\nSelect a unique age for better results')
-                    elif len(age) == 0:
-                        print('No layer provided, please provide a valid age.')
-                        return
-                elif var == 'IRH_FRAC_DEPTH':
-                    if cmap == None:
-                        cmap = cmaps.torch_r
-                    extend = 'both'
-                    age = list(metadata['age'])
-                    if len(age) > 1:
-                        print('WARNING: Multiple layers provided: ', age,
-                            '\nSelect a unique age for better results')
-                    elif len(age) == 0:
-                        print('No layer provided, please provide a valid age.')
-                        return
-                    age = age[0]
+                if len(age) > 1:
+                    print('WARNING: Multiple layers provided: ', age,
+                        '\nSelect a unique age for better results')
+                elif len(age) == 0:
+                    print('No layer provided, please provide a valid age.')
+                    return
+                if fraction_depth:
                     label = r'IRH Fractional Depth [\%]'
+                    if vmin is None: vmin = 0
+                    if vmax is None: vmax = 100
+                else:
+                    label = f'IRH Depth [m]'
 
                 for age in metadata['age']:
-                    scatter = plt.scatter(df['PSX']/1000, df['PSY']/1000, c=df[int(age)], cmap=cmap, s=marker_size, vmin=vmin, vmax=vmax, linewidths=0, rasterized=True)
+                    scatter = plt.scatter(df['PSX']/1000, df['PSY']/1000, c=df[age], cmap=cmap, s=marker_size, vmin=vmin, vmax=vmax, linewidths=0, rasterized=True)
 
         if color_by == 'transect_1D':
             flight_id = list(metadata['flight_id'])
@@ -605,13 +601,14 @@ class Plotting:
                         all_y.append(ds['PSY'][:])
                 df = pd.DataFrame({'PSX': np.concatenate(all_x),
                                 'PSY': np.concatenate(all_y)})
-                plt.scatter(df['PSX'][::downscale_factor]/1000, df['PSY'][::downscale_factor]/1000, color=colors[flight_id], s=marker_size, linewidths=0)
+                plt.scatter(df['PSX'][::downsampling_factor]/1000, df['PSY'][::downsampling_factor]/1000, color=colors[flight_id], s=marker_size, linewidths=0)
 
                 plt.plot([], [], color=colors[flight_id], label=flight_id, linewidth=3)
             ncol = 2 if len(flight_ids) > 40 else 1
 
         # --- Format Figure ---
-        print('Formatting ...')
+        if not self._disable_tqdm:
+            print('Formatting ...')
 
         ax.set_xlim(*xlim)
         ax.set_ylim(*ylim)
@@ -710,7 +707,7 @@ class Plotting:
             site_coords = pd.read_pickle(self._site_coords_path)
             for i in site_coords.index:
                 site = site_coords.loc[i]
-                ax.scatter(site['x']/1000, site['y']/1000, color='red', s=50, marker='^', edgecolor='black', linewidth=1.5)
+                ax.scatter(site['x']/1000, site['y']/1000, color='red', s=50, marker='^', edgecolor='black', linewidth=1.5, zorder=50)
         # --- Save/Show ---
         with self._plot_context():
             if save:

{anta_database-0.3.3 → anta_database-0.3.5}/anta_database.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: anta_database
-Version: 0.3.3
+Version: 0.3.5
 Summary: SQLite database for the AntArchitecture Community Data
 Author: Antoine Hermant
 Author-email: antoine.hermant@etik.com
@@ -48,12 +48,12 @@ Visit the [Home Page](https://antoinehermant.github.io/anta_database/intro) for
 
 The AntADatabase contains all traced and dated Internal Reflective Horizons (isochrones) that have been published at that date. All the data is organized by flight transect for each dataset. The data is encoded in HDF5 (Hierarchical Data Format) which provides convenient data structure and performant read speeds. 
 Variables per file (when exist):
-- x(point)
-- y(point)
-- distance(point)
+- PSX(point)
+- PSY(point)
+- Distance(point)
 - IRH_DEPTH(point, age)
-- IRH_DENS(point, age)
-- ICE_THCK(point, age)
+- IRH_NUM(point, age)
+- ICE_THK(point)
 - SURF_ELEV(point, age)
 - BED_ELEV(point, age)
 
@@ -66,7 +66,7 @@ This Python module provides SQL indexing for the AntADatabase, as well as quick
 - age
 - acquisition_year (year of radar acquisition)
 - region
-- basin
+- IMBIE_basin
 - var (variable)
 - flight id
 

{anta_database-0.3.3 → anta_database-0.3.5}/anta_database.egg-info/SOURCES.txt

@@ -41,6 +41,7 @@ book/contact.md
 book/contribute.md
 book/intro.md
 book/logo.png
+book/pism_example.ipynb
 book/quick_start.ipynb
 book/references.bib
 book/requirements.txt

{anta_database-0.3.3 → anta_database-0.3.5}/book/_config.yml

@@ -9,6 +9,7 @@ logo: logo.png
 # See https://jupyterbook.org/content/execute.html
 execute:
   execute_notebooks: force
+  timeout: 600  # Timeout in seconds (e.g., 600 = 10 minutes)
 
 # Define the name of the latex output file for PDF builds
 latex:

{anta_database-0.3.3 → anta_database-0.3.5}/book/_toc.yml

@@ -5,6 +5,7 @@ format: jb-book
 root: intro
 chapters:
 - file: quick_start
+- file: pism_example
 - file: advanced
 - file: contact
 - file: contribute

{anta_database-0.3.3 → anta_database-0.3.5}/book/contribute.md

@@ -6,3 +6,5 @@ Maybe first contact me so I know that people are interested and on which feature
 
 The project is hosted on [GitHub](https://github.com/antoinehermant/anta_database). First you should fork the repository, create your own features branches and submit pull request.
 
+Check the [Roadmap](https://github.com/antoinehermant/anta_database/blob/main/ROADMAP.md).
+