Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1479 1656 (N29W54 L=251)   526     81        0 (1440/1631-1522/1684) -63.4 W51-W55 N25-N31
      2 1576 1544 (N20W53 L=250)   392     65        0 (1537/1532-1615/1558) -26.6 W49-W56 N18-N21
      3 1545 1397 (N13W42 L=240)   301     55        0 (1523/1369-1571/1425)  21.8 W40-W45 N11-N14
      4  207 1386 (N47E68 L=130)   886     97        1 ( 161/1345- 244/1404)  -6.3 E85-E58 N43-N48
      5 1604 1360 (N09W45 L=243)   312     71        1 (1565/1335-1633/1389)  -9.5 W41-W47 N08-N10
      6  802 1282 (N27E07 L=191)   533    210        1 ( 706/1255- 906/1308) -30.3 E12-W00 N23-N30
      7  394 1136 (N28E39 L=159)   304     50        0 ( 365/1117- 429/1150)   0.0 E41-E35 N26-N28
      8 1052 1093 (N10W03 L=201)   373    112        1 (1011/1047-1098/1129) -74.1 W01-W04 N06-N13
      9  730 1072 (N17E16 L=182)   194    128        0 ( 716/1021- 767/1114)  51.3 E16-E12 N13-N18
     10 1644  909 (S16W35 L=232)   254     65        1 (1623/ 877-1664/ 938)  56.3 W32-W36 S17-S14
     11 1241  695 (S17W04 L=202)   631    180        1 (1186/ 623-1292/ 757)  40.6 W00-W08 S21-S14
     12  207  647 (N02E71 L=126)  1132     95        1 ( 176/ 602- 242/ 675) -63.4 E76-E66 S01-N03
     13 1616  453 (S47W28 L=226)  1286    150        2 (1551/ 419-1677/ 504)   4.1 W21-W35 S49-S44
     14  761  330 (S28E38 L=159)  2274    407        3 ( 721/ 155- 805/ 501)  39.8 E57-E27 S43-S17
     15  394  426 (S15E63 L=135)  1198    113        6 ( 365/ 364- 420/ 479)  36.0 E69-E58 S19-S11
     16  494  417 (S17E53 L=144)   749     96        3 ( 455/ 379- 538/ 458) -51.3 E55-E50 S19-S13
     17  977  407 (S28E21 L=177)   281    115        0 ( 951/ 356-1005/ 450)  35.5 E24-E17 S31-S26
     18  495  310 (S24E63 L=134)   944    105        4 ( 446/ 288- 543/ 331) -29.7 E68-E57 S26-S22
  

    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