Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1089 1820 (N48E26 L=359)   343     53        1 (1058/1807-1117/1833)  26.6 E28-E24 N45-N49
      2 1202 1621 (N36E08 L= 18)  1378    313        3 (1131/1497-1282/1734) -39.1 E16-W00 N28-N42
      3  404 1685 (N19E68 L=318)  2457    157        3 ( 361/1639- 450/1730)  63.4 E71-E64 N14-N22
      4  583 1606 (N17E46 L=340)   596    200        0 ( 561/1533- 609/1684) -48.4 E49-E40 N13-N22
      5  448 1550 (N11E53 L=332)   715    169        2 ( 427/1487- 482/1612) -71.6 E56-E50 N06-N15
      6 1790 1408 (N43W46 L= 71)   462    103        1 (1767/1384-1816/1426)  33.7 W43-W49 N40-N45
      7 1658 1266 (N27W33 L= 59)   292     75        0 (1639/1223-1671/1301) -59.0 W31-W34 N24-N29
      8 1184 1193 (N09W04 L= 30)   221     51        1 (1164/1162-1207/1224) -76.0 W03-W05 N06-N10
      9 1821  912 (N12W56 L= 81)  1771    324        4 (1779/ 783-1867/1021) -36.3 W49-W65 N06-N17
     10 1443  983 (N04W25 L= 50)   259     78        1 (1412/ 965-1476/1003)  53.1 W22-W26 N02-N06
     11 1887  924 (N15W63 L= 88)   524    120        1 (1870/ 885-1914/ 948)  59.0 W60-W66 N12-N17
     12 1237  799 (S12W18 L= 43)   432    166        0 (1165/ 736-1324/ 833)  56.3 W15-W22 S16-S07
     13  185  802 (S39E57 L=328)   492    110        1 ( 159/ 781- 209/ 827)  26.6 E60-E54 S40-S37
     14  989  672 (S26W08 L= 33)   613    166        1 ( 900/ 654-1077/ 696)  31.0 W02-W12 S28-S21
     15 1159  521 (S29W24 L= 50)   308     51        1 (1127/ 506-1196/ 544)  26.6 W22-W26 S30-S28
     16  896  473 (S41W10 L= 35)   308     59        0 ( 870/ 455- 931/ 498)   0.0 W07-W12 S42-S39
  

    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