Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  485 1493 (N33E36 L= 12)  4028    557       10 ( 309/1288- 661/1644) -28.3 E72-E20 N18-N47
      2 1689 1494 (N22W55 L=103)   928    131        0 (1623/1452-1739/1515) -24.0 W48-W60 N19-N24
      3  246 1263 (N20E56 L=352)   912    169        2 ( 221/1216- 276/1321)  81.9 E58-E52 N17-N24
      4  307 1247 (N19E49 L=359)   362     80        0 ( 292/1205- 328/1279)  53.1 E50-E46 N16-N20
      5  178 1152 (N14E63 L=345)  3107    295        8 ( 116/1056- 224/1244) -36.9 E75-E56 N08-N20
      6  260 1057 (N08E53 L=355)   404     77        1 ( 230/1025- 282/1099)  68.2 E55-E50 N05-N10
      7 1694  998 (S08W44 L= 92)   820    182        1 (1660/ 914-1752/1047)  41.2 W41-W49 S12-S05
      8 1453  766 (S20W26 L= 73)   382    140        0 (1398/ 735-1525/ 790)  12.5 W21-W30 S22-S19
      9  508  564 (S25E42 L=  6)  1717    425        8 ( 442/ 404- 629/ 760)  36.9 E55-E27 S34-S13
     10 1271  691 (S24W13 L= 61)   597    217        4 (1183/ 660-1358/ 751) -26.6 W07-W19 S26-S19
     11  348  597 (S21E55 L=353)  1083    178        1 ( 312/ 510- 379/ 670)  68.2 E57-E51 S26-S16
     12 1562  472 (S42W40 L= 87)   391     57        1 (1529/ 456-1597/ 488) -12.5 W35-W44 S43-S40
  

    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