Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1653 1458 (N20W53 L=186)   440    126        0 (1612/1434-1695/1472)   0.0 W48-W57 N19-N21
      2 1198 1201 (N14W13 L=146)   543    141        1 (1171/1131-1223/1270)  77.5 W11-W14 N09-N18
      3 1654 1035 (S04W41 L=174)   349     72        0 (1616/1009-1683/1061)   9.5 W37-W43 S05-S02
      4 1704  871 (S15W43 L=176)   695    145        0 (1688/ 803-1721/ 953)  77.5 W41-W44 S19-S10
      5  322  686 (S06E55 L= 78)   429     65        1 ( 304/ 661- 349/ 718)  53.1 E57-E52 S07-S03
      6 1318  660 (S21W13 L=145)   231     48        0 (1306/ 625-1333/ 691)  78.7 W11-W14 S23-S18
      7 1603  642 (S28W34 L=166)   323     57        1 (1583/ 611-1624/ 663)  56.3 W32-W35 S29-S26
      8  816  513 (S22E22 L=110)  1126    254        2 ( 712/ 460- 932/ 566)  13.2 E30-E13 S25-S19
      9  453  503 (S19E50 L= 83)   353     62        0 ( 424/ 468- 480/ 526) -63.4 E51-E47 S21-S17
     10  584  325 (S34E48 L= 85)   719    128        1 ( 552/ 294- 610/ 361)  32.0 E51-E43 S36-S30
  

    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