Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1415 1657 (N49W21 L= 45)   305     77        0 (1381/1646-1452/1669)   8.1 W17-W24 N48-N50
      2 1030 1548 (N33E09 L= 15)  1667    340        5 ( 870/1435-1233/1635) -10.1 E22-W06 N28-N36
      3  615 1552 (N26E38 L=346)   898    205        5 ( 573/1470- 662/1622) -78.7 E41-E35 N20-N30
      4 1587 1466 (N37W35 L= 59)   392    100        2 (1536/1451-1630/1481)   7.1 W30-W38 N36-N38
      5  569 1397 (N15E36 L=348)   369    163        0 ( 532/1344- 624/1446) -29.1 E39-E30 N12-N18
      6 1085 1376 (N22E02 L= 22)   388    139        0 (1010/1348-1150/1396)   7.1 E06-W02 N20-N23
      7 1511 1287 (N23W28 L= 52)   373     64        2 (1485/1264-1544/1315) -26.6 W26-W30 N21-N25
      8  593  915 (S13E25 L=359)   484    153        1 ( 548/ 862- 652/ 976) -26.6 E29-E20 S16-S10
      9 1330  422 (S33W34 L= 58)   498    118        0 (1289/ 401-1377/ 451)  -8.1 W30-W37 S34-S31
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y