Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  595 1667 (N21E48 L=115)  1888    351        1 ( 564/1542- 637/1784) -54.8 E54-E42 N13-N30
      2  373 1536 (N08E60 L=103)  1849    482        2 ( 346/1381- 407/1688) -46.6 E69-E50 N00-N18
      3  239 1377 (S04E66 L= 97)   832    104        1 ( 212/1342- 271/1414)  68.2 E69-E62 S05-N00
      4  739 1378 (N07E26 L=137)   460     70        0 ( 693/1347- 780/1410)  63.4 E28-E24 N04-N10
      5 1125 1306 (N13E02 L=161)   370    135        0 (1072/1240-1163/1368) -35.5 E06-W01 N10-N15
      6  898 1341 (N09E16 L=147)   250     72        0 ( 855/1326- 932/1362)  53.1 E17-E14 N07-N11
      7 1745 1148 (N23W43 L=206)   580    111        1 (1726/1103-1761/1211) -63.4 W40-W44 N20-N26
      8  346 1030 (S22E44 L=119)   338     75        2 ( 322/ 998- 373/1058)  90.0 E45-E41 S23-S19
      9  623  861 (S25E20 L=143)  1397    457        3 ( 446/ 810- 827/ 937)  38.0 E31-E08 S32-S14
     10 1124  777 (S16W13 L=176)   463    170        0 (1041/ 742-1208/ 798)  33.7 W08-W17 S19-S13
     11 1842  735 (N04W66 L=230)   839    119        1 (1824/ 682-1857/ 777) -56.3 W62-W69 N00-N06
     12 1174  681 (S20W19 L=182)   937    258        1 (1061/ 656-1315/ 709)  28.1 W11-W26 S22-S14
     13 1513  548 (S17W45 L=208)   460    103        2 (1467/ 523-1560/ 566)  33.7 W41-W49 S19-S14
     14  808  509 (S41W02 L=165)   241     45        1 ( 787/ 484- 832/ 527) -33.7 W00-W03 S43-S39
     15  985  485 (S37W15 L=178)   239     67        0 ( 962/ 460-1017/ 501)  63.4 W13-W16 S39-S35
     16 1457  287 (S33W58 L=221)   535    103        2 (1433/ 260-1481/ 310)   0.0 W54-W60 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