I figured out it was some simple maths to do. Using a scaling matrix [4,0,0,0,4,0,0,0,1]*[x,y,1]. So after multiplying we get[4x,4y,1]. So if we multiply all points the image is scaled.Learnt this in school. Also the map was inverted. So my program finally looks like this
from OpenGL.GL import *
from OpenGL.GLUT import *
from OpenGL.GLU import *
import math
def initFun():
glClearColor(1.0,1.0,1.0,0.0)
glColor3f(0.0,0.0, 0.0)
glMatrixMode(GL_PROJECTION)
glLoadIdentity()
gluOrtho2D(0,700,0,700)
def displayFun():
glClear(GL_COLOR_BUFFER_BIT)
glColor3f(0.0,0.0,1.0)
xpts=[]
ypts=[]
N=len(newlist)
for i in range(0,N):
xpts.append(700+100-float(newlist[i][4])*20)
ypts.append(700+1300-float(newlist[i][5])*20)
glBegin(GL_POINTS)
for j in range(0,N):
glVertex2f(xpts[j],ypts[j])
glEnd()
glFlush()
if __name__ == '__main__':
glutInit()
glutInitWindowSize(700,700)
glutCreateWindow("My Display")
glutInitDisplayMode(GLUT_SINGLE | GLUT_RGB)
glutDisplayFunc(displayFun)
initFun()
glutMainLoop()