Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  487 1537 (N51E33 L=246)   348     49        1 ( 467/1509- 510/1562)  45.0 E36-E30 N48-N52
      2 1241 1513 (N30W27 L=307)  2527    735        3 ( 936/1416-1553/1549) -28.9 W06-W45 N16-N38
      3 1869 1187 (S06W65 L=344)  1839    190        6 (1817/1135-1933/1259) -35.5 W58-W76 S10-N00
      4 1239 1183 (N11W16 L=296)   214     87        0 (1195/1149-1271/1207) -53.1 W14-W17 N08-N13
      5  328 1089 (N24E47 L=233)  1267    247        1 ( 294/1002- 377/1204) -69.0 E50-E43 N18-N31
      6  827 1089 (N15E10 L=269)   938    250        2 ( 800/ 959- 862/1189)  81.9 E11-E07 N06-N20
      7 1663 1080 (S05W41 L=321)   274     64        1 (1641/1047-1687/1104) -63.4 W39-W42 S07-S03
      8 1775  788 (S26W47 L=327)  1337    216        3 (1703/ 692-1831/ 867)  28.3 W39-W54 S31-S21
      9  423  710 (S00E46 L=234)   809    210        0 ( 394/ 613- 483/ 762)  28.6 E51-E39 S05-N02
     10  359  530 (S11E59 L=220)  7624      0       17 ( 259/ 385- 464/ 662) -75.3 E72-E50 S21-S01
     11 1304  400 (S39W04 L=283)   361     57        0 (1279/ 373-1331/ 427)  36.9 W01-W05 S41-S37
  

    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