Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  454 1665 (N25E57 L=148)   449     59        0 ( 433/1641- 475/1692) -76.0 E59-E54 N23-N27
      2 1582 1586 (N41W32 L=237)   741    110        0 (1535/1544-1623/1623) -26.6 W27-W35 N38-N43
      3  301 1493 (N12E62 L=143)   609     77        2 ( 266/1463- 337/1520)  38.7 E65-E58 N10-N14
      4 1070 1434 (N19E04 L=201)   266     71        0 (1033/1415-1109/1451)  11.3 E06-E01 N17-N19
      5 1436 1422 (N25W19 L=224)   277     88        0 (1394/1398-1474/1449)  -9.5 W16-W22 N24-N26
      6 1645 1326 (N24W36 L=241)   441    111        0 (1612/1276-1681/1373) -45.0 W32-W38 N21-N26
      7  138 1109 (S12E69 L=136)  3747    165        5 (  88/1045- 196/1191)  29.7 E80-E62 S16-S07
      8  326  804 (S30E47 L=158)  1293    197        3 ( 236/ 771- 429/ 834)  10.8 E58-E37 S32-S27
      9 1839  670 (S09W69 L=274)   660     76        0 (1805/ 640-1875/ 689)  26.6 W65-W74 S11-S07
     10 1486  506 (S28W44 L=249)  2765    511        2 (1278/ 446-1724/ 577)  30.1 W29-W61 S35-S17
  

    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