Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  625 1526 (N35E23 L= 63)   505     81        0 ( 587/1493- 659/1550) -33.7 E26-E20 N32-N37
      2 1190 1425 (N22W15 L=101)   707    238        0 (1064/1404-1301/1439) -14.9 W06-W21 N19-N23
      3  169 1235 (N21E65 L= 21)  1192    100        2 ( 147/1183- 202/1292) -45.0 E70-E59 N17-N24
      4 1800 1220 (N03W57 L=143)   911     94        4 (1761/1171-1837/1271) -80.5 W52-W60 N00-N06
      5 1602  900 (S14W36 L=122)   312     31        0 (1585/ 873-1619/ 927)  71.6 W34-W37 S15-S12
      6  849  800 (S12E13 L= 73)   249     86        1 ( 804/ 771- 885/ 834)  18.4 E16-E10 S13-S10
      7 1570  778 (S21W34 L=120)   413     95        0 (1534/ 751-1603/ 822)  45.0 W31-W36 S22-S18
      8 1068  679 (S22E01 L= 85)   278    115        0 (1014/ 655-1119/ 712) -38.7 E03-W02 S24-S19
      9  756  584 (S25E23 L= 63)   405     64        2 ( 720/ 565- 802/ 609) -26.6 E25-E19 S26-S23
     10 1136  411 (S41W01 L= 87)   834    164        1 (1056/ 369-1219/ 441)   9.5 E05-W07 S44-S39
     11  968  354 (S44E14 L= 72)   346     76        0 ( 930/ 335- 993/ 386)  26.6 E17-E11 S45-S41
  

    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