Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  415 1386 (N30E45 L=194)  1502    373        3 ( 328/1265- 546/1450)  24.0 E52-E33 N21-N34
      2  979 1397 (N28E01 L=239)   339    178        0 ( 893/1365-1041/1418) -12.5 E06-W03 N26-N30
      3  260 1377 (N29E64 L=176)   674     96        1 ( 220/1353- 300/1395)  10.3 E69-E58 N27-N30
      4 1409 1364 (N23W29 L=268)   235    109        0 (1367/1335-1461/1389) -33.7 W26-W32 N20-N25
      5  466 1321 (N25E39 L=200)   682    138        1 ( 403/1297- 533/1339)  10.3 E44-E33 N24-N26
      6  197 1184 (N16E65 L=174)  1612    212        7 ( 162/1089- 221/1286) -63.4 E72-E60 N10-N22
      7  849 1227 (N18E10 L=229)   294    107        2 ( 816/1191- 886/1272)  45.0 E12-E07 N16-N20
      8 1303  929 (S03W17 L=256)   503    168        0 (1242/ 882-1349/ 990)  32.0 W12-W20 S05-N00
      9 1774  882 (S10W53 L=293)  1253    210        4 (1719/ 820-1823/ 963)  41.6 W48-W58 S13-S05
     10  996  597 (S22E04 L=235)   215    128        0 ( 943/ 565-1039/ 635)  35.5 E08-E01 S24-S19
     11  767  566 (S23E20 L=219)   644    153        1 ( 747/ 503- 788/ 634)  77.5 E21-E18 S27-S18
     12  538  557 (S23E38 L=202)   417    166        0 ( 482/ 526- 605/ 604)  18.4 E43-E31 S25-S20
  

    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