This paper proposes an approach for handling multivariate data in an archaeological Geographical Information System (GIS), providing a new tool to archaeologists and historians. Our method extracts potential objects of known shapes in a geographical database (GDB) devoted to archaeological excavations. In this work, archaeological information is organized according to three components: location, date and a shape parameter, in a context where data are imprecise and lacunar. To manage these aspects, a three-step methodology was developed using fuzzy sets modeling and adapting the fuzzy Hough transform. This methodology is applied in order to define the appropriate tool for a GDB of Roman street remains in Reims, France. The defined queries return an estimation of the possible presence of streets during a fuzzy time interval given by experts on the Roman period in Reims.