Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1450 1701 (N45W28 L= 18)   412     90        0 (1417/1659-1474/1725) -59.0 W25-W29 N41-N47
      2 1471 1546 (N32W27 L= 17)   249     31        0 (1448/1528-1494/1562)   0.0 W24-W28 N31-N33
      3  726 1485 (N21E24 L=326)   258     51        0 ( 695/1469- 757/1499)  26.6 E25-E21 N19-N21
      4 1427 1356 (N19W23 L= 13)   392    136        0 (1369/1333-1496/1383)  -6.3 W18-W27 N17-N20
      5 1124 1283 (N11W04 L=354)  1131    298        3 (1005/1183-1258/1360) -26.6 E04-W12 N06-N15
      6 1670 1288 (N18W41 L= 31)   913    177        0 (1634/1188-1695/1358) -73.3 W37-W42 N11-N21
      7 1818 1134 (N10W54 L= 44)   505     99        3 (1791/1102-1846/1164)  45.0 W51-W57 N08-N12
      8 1723 1112 (N08W45 L= 35)   263     82        0 (1687/1099-1753/1126)  11.3 W42-W47 N06-N08
      9  654 1057 (S06E23 L=327)   237     84        0 ( 610/1038- 700/1075)  26.6 E25-E19 S07-S04
     10  941  951 (S10E04 L=346)   353    143        2 ( 875/ 924-1016/ 990) -14.0 E08-W00 S11-S08
     11  569  696 (S29E28 L=322)  1801    675        6 ( 351/ 561- 797/ 902)  33.0 E49-E12 S38-S14
     12 1559  507 (S31W45 L= 35)   885    221        2 (1509/ 452-1614/ 578) -28.3 W38-W51 S33-S26
  

    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