Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1337 1736 (N51W39 L=282)   538     59        0 (1312/1703-1355/1765)  78.7 W35-W41 N48-N53
      2 1844 1319 (N16W70 L=313)  1443    142        2 (1782/1276-1896/1363) -20.6 W62-W79 N13-N19
      3  469 1155 (N14E37 L=206)   291     46        1 ( 438/1134- 492/1176) -26.6 E39-E35 N12-N15
      4  675  872 (S04E23 L=221)   569    248        0 ( 577/ 828- 783/ 935) -24.8 E29-E15 S06-N00
      5 1081  750 (S13W02 L=246)   351    113        0 (1045/ 712-1136/ 776)  31.0 W00-W05 S15-S11
      6 1600  663 (S22W39 L=283)   478     56        1 (1564/ 641-1627/ 685) -21.8 W36-W41 S23-S20
      7  299  658 (S17E57 L=186)   549     39        1 ( 279/ 635- 330/ 682)  15.9 E60-E53 S18-S15
      8  291  594 (S22E62 L=181)   540     68        1 ( 273/ 562- 311/ 624)  76.0 E65-E59 S23-S19
      9 1478  549 (S29W31 L=274)   478    136        0 (1417/ 513-1528/ 590)  35.5 W26-W34 S31-S26
     10  715  532 (S26E24 L=219)   371    150        0 ( 673/ 479- 779/ 571)  33.7 E28-E19 S29-S23
  

    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