Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  616 1637 (N30E37 L=192)   350     61        0 ( 586/1620- 648/1652)  -9.5 E39-E33 N28-N30
      2  770 1514 (N22E22 L=207)   381    116        1 ( 731/1478- 825/1548) -29.7 E24-E17 N20-N24
      3  446 1244 (N03E39 L=190)   442     68        0 ( 426/1198- 467/1288)  78.7 E40-E38 N00-N05
      4  678 1144 (S02E22 L=207)   194     38        0 ( 660/1120- 700/1166)  63.4 E23-E20 S02-N00
      5  755 1113 (S03E17 L=212)   193     54        0 ( 726/1091- 789/1136) -26.6 E18-E14 S03-S01
      6 1385  810 (S14W25 L=254)  1563    314        3 (1254/ 715-1530/ 941)  48.8 W17-W32 S20-S04
      7  802  917 (S14E12 L=216)   190     40        1 ( 774/ 904- 827/ 935)   0.0 E14-E10 S14-S12
      8  627  779 (S24E23 L=206)   248     70        0 ( 596/ 764- 669/ 792)   0.0 E25-E20 S24-S22
      9  785  691 (S28E12 L=217)   447    121        0 ( 718/ 649- 840/ 731)  36.9 E16-E08 S31-S24
     10 1569  656 (S21W41 L=270)   322     89        0 (1538/ 613-1591/ 699)  80.5 W40-W42 S23-S17
     11 1030  362 (S48W10 L=239)   262     54        0 ( 996/ 348-1069/ 378)  18.4 W07-W13 S49-S46
     12  647  346 (S55E27 L=202)   396     42        2 ( 622/ 327- 674/ 367)   0.0 E29-E23 S56-S52
  

    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