def METHOD_NAME(tmpdir, gdf):
    # Create a geodataset and an ogr compliant version of it
    ds = GeoDataset.from_gdf(gdf)
    oc = ds.vector.ogr_compliant()
    # Assert some ogr compliant stuff
    assert oc.ogr_layer_type.upper() == "MULTIPOLYGON"
    assert list(oc.dims)[0] == "index"
    assert len(oc.Roman) == 2
    # Write and load
    fn = str(tmpdir.join("dummy_ogr.nc"))
    ds.vector.to_netcdf(fn, ogr_compliant=True)
    ds1 = GeoDataset.from_netcdf(fn)
    assert np.all(ds.vector.geometry == ds1.vector.geometry)