Tag Archives: Código Aberto/Open Source

LibreOffice Suite: All your research needs, for free

Finally, after much delay due to my dissertation work, I can resume my work on the role free-software can play in archaeology, in this second entry of Open source archaeology, a brief users guide (II). This is a topic I really enjoy writing about, because there are a lot of tools out there created and updated by devoted groups of supporters that just want to provide us with programs that are not restrained by the market rules and, whether you want to believe or not, they can be as good, if not better, than firmware. In the case of LibreOffice, at least for me, this is absolutely true.

LibreOffice is a project of The Document Foundation, a forked version from OpenOffice, and comprises programs for general office work, such as word processing, the creation and editing of spreadsheets, slide-shows, diagrams and drawings, mathematical formulae and, most important for archaeologists, databases. You can work in the native formats or with other software formats such as .doc, and there is no issue, at least on my experience, to shift from one to the other.

Continue reading LibreOffice Suite: All your research needs, for free


Open source archaeology, a brief users guide (I)

A digital world like ours has many initiatives that perhaps we do not know, but that offer us open-source software for any task we can imagine. From whole operative systems (OS) to little programs for specific tasks, communities of coders work for free on alternatives to firmware software that in some cases can make us wonder why we still pay for programs when we can donate to this initiatives and enjoy great software. My work on GIS and GRASS was of course my first exposure to this technology but my best example would be LibreOffice: I have not used the Office Suit for more than a decade, and I have never encountered a problem on my research or on my publications.

App menu of ArcheOS. Look at all the archeological software ready to be used!

Continue reading Open source archaeology, a brief users guide (I)


Time has passed since my last post on GIS, sorry for the waiting. But here we are on track again, this time with an essential one… how to add the points from your GPS to your map. Remember you can download all tutorials from my Academia.edu page. Enjoy!

Vector points:

So now that we have our raster map of the region we want to study it is time to add some sites to it. In order to do this we will create a vector map from a list of points in a spreadsheet.

1.- Once you have a list of points you can organise them in a spreadsheet file like this:

GIS 5-1

The most important things to take into account are the order to the fields, which is not mandatory but strongly suggested, and that in order to avoid issues it is better to avoid spaces at least on the headers. This is a simplified version from my Ph.D. research, and it gives you the general structure. We need only three main fields, the two coordinates and the CAT, which is just a field from 1 to n, like the index on a database. Names, types of sites and many other fields can be added after this. It is up to you include as many as you want, but for the lesson of the day we will use this simplified version.

Continue reading INTRODUCTION TO GIS (V)

Introduction to GIS (IV):

Reprojecting to UTM:

So now we have a map that covers most of our working area, but the coordinates are in LatLong. The problem is that this is a sexagesimal system, and thus not very good for some of the operations we want to carry out. So, before we can finish our base map we need to covert it to UTM, a decimal system. It is a simple operation, once you get used to it, but in my experience it takes a couple of attempts before the reprojections are done with confidence. Leave any questions if some of the steps are not clear.

1.- Open your patched map. In Settings ~> Region ~> Set region we will set the region of the map. A region, in simple words, represents the boundaries of your maps, its mathematical limits. It is important to always know where the boundaries of your map are, since any calculation will end at the limits of the region, and many times those are not the same that the ones in the map. In the menu select Set region to match raster map and choose your patched map in the drop list. Then press Close.


2.- By doing this we are able now to create a vector map that will be our “template” to be able to reproject the map. In Vector ~> Generate area for current region. In the Required tag of the menu create a Name for the map (remember to avoid spaces). In Optional be sure the option of As area is selected. Run and Close. If you check this new map you will see that it is just a grey box exactly matching the limits of your map. That is what we need. We need to close GRASS and open it again.

3.- In order to wok with Utm we need to create a new location for our map. It is the same method we learnt in lesson one, but this time we need to Select coordinate system parameters from a list. In the next window, look for UTM in Choose Projection.


4.- What we have done with all those steps was to create one location that fits our needs, rather than extracting the information from a file. But we still need to “translate” our LatLong map to the UTM location. For this we access the permanent mapset and we will use the vector map we have created before to establish the boundaries of our map.

5.- In Vector ~> Develop raster map ~> Reproject vector map from different GRASS location. In Required select the location in which you have made the vector map. In Source select first Mapset and then name of input map. Press Run.



6.- Now you have reprojected the vector map to UTM. Perhaps you have notice that it is not a rectangle but a trapezoid. That is actually what you are looking for, because it indicates that the map was reprojected properly.

7.- Set the region as in step 1, but this time Set region to match vector map.


8.- We are almost there. The last step is to reproject the raster from the Latlong location to the UTM location. To do this, go Raster ~> Develop raster map ~> Reproject raster map from different GRASS location. The Required and Source steps are the same that those for the vector. The difference in a raster is that we need to choose a method of interpolation (lanczos_f) and indicate the resolution of output raster map (30). In simple words the method is the way the software calculates the value of the cells of the new map. Other methods such as nearest or bicubic are faster, but they affect the results of the output map. Lanczos_f allows us to do the task without loosing information. The reason for selecting 30 in the resolution is because it is the same that our source map (ignore bicubic on the picture, choose lanczos_f option).

And that is all! I know this task can be a little bit confusing, but it is essentiall for your work on maps. Now you have your map in UTM units, ready to work and to be used on some operations we will exploring soon. If you have any doubt please leave a comment!

Previous post

Introduction to GIS (III):

Importing and Patching Maps:

In the last lesson we learnt how to import our first map in GRASS. However, as it happens many times, our working area needs of more than one map. In this lesson we will import several new DEM maps to our location and we will learn to patch them all together in a single, larger map. I have used several .tif files that I have uploaded here. Remember to change the extension to .rar (follow the instructions of lesson II).

1.- In the Layer Manager window, select File ~> Import raster data ~> Common formats import.

2.- In the menu, select File if you need just one map, or Directory if you need more. If you choose the second option be sure the Format selected is GeoTiff and the Extension .tif. Be sure all the boxes on the left side of your map names are selected. Finally, select Extend region extents based on new dataset. Press Run. When finished, press Close.

It is time now to introduce a very important concept: Region. In the previous lesson we decided to set up the region based on the DEM file. Imagine a Region is like the boundaries of your location. It has an N, S, E, and W limit, right now just the ones of the map we used in lesson 2. But if we add new maps without changing that region, those boundaries of our location, we will just pile up one on top of the other. That is why we select Extend region extents based on new dataset, so the new maps will push those boundaries to fit the new information we are adding.


3.- Now you should be able to see the maps you have added. If you include the map we used in lesson 2 you will have a rather decent portion of central Greece and northern Peloponnese in your display. But this is not yet a single map but just all the ones we have imported shown together. In order to patch different maps into a single one go to Raster ~> Overlay rasters ~> Patch raster maps.


4.- In Required select the maps to be added and a new name for your map. When adding several maps you need to press the arrow on the right to select them one by one, so the software will create a large string of names in that text box. It is a little bit counterintuitive, but just go one by one carefully and you should not have any problem. When ready press Run.

Note the long string of names in the [multiple] section.

And done. From several maps now we have a single one that represents our working area. Change the colours as we have seen in lesson 2 and try to play a little bit on the display window exploring the options and zooming in and out. As you can see in the bottom left of the screen, the map is expressed in latitude and longitude, which is a sexagesimal system and more complicated to operate with. So next lesson we will learn how to reproject your map from latlong to UTM.

DEM map of the Peloponnese

Previous post

Next post