Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1147 1775 (N43E19 L= 18)  2465    435        5 (1023/1636-1293/1925) -20.4 E45-E02 N36-N54
      2  512 1708 (N21E59 L=337)  3645    292       13 ( 428/1619- 606/1791)  51.3 E68-E53 N13-N27
      3  887 1716 (N30E31 L=  5)   293     65        0 ( 871/1677- 903/1751) -56.3 E33-E29 N27-N33
      4  632 1712 (N23E49 L=348)   566     61        1 ( 603/1687- 659/1735)  76.0 E51-E46 N21-N25
      5  598 1513 (N10E40 L=357)  1359    427        2 ( 468/1375- 726/1658)  76.6 E44-E36 N00-N21
      6 1731 1441 (N41W39 L= 75)   481     44        1 (1703/1412-1753/1468) -76.0 W35-W41 N38-N43
      7  404 1070 (S20E40 L=356)   297     48        1 ( 384/1044- 424/1102)  76.0 E41-E38 S21-S17
      8  334 1043 (S22E46 L=351)   609    168        1 ( 302/ 946- 360/1099) -48.8 E49-E42 S27-S19
      9  489 1004 (S21E32 L=  4)   268     50        0 ( 455/ 988- 521/1019)  26.6 E34-E30 S23-S20
     10 1256  235 (S44W55 L= 91)   501    103        0 (1219/ 208-1305/ 265)   4.1 W48-W62 S45-S42
  

    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