Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  677 1692 (N27E40 L=182)   476     95        2 ( 659/1640- 703/1734) -60.3 E43-E37 N23-N30
      2  571 1661 (N24E47 L=176)   457     80        1 ( 539/1638- 610/1683)  11.3 E50-E44 N21-N25
      3  925 1528 (N21E17 L=206)   755    222        0 ( 903/1433- 954/1631) -65.2 E20-E12 N15-N28
      4 1520 1514 (N35W23 L=245)  1223    338        1 (1408/1376-1606/1620) -35.2 W12-W30 N26-N39
      5  334  928 (S24E45 L=177)  1099    205        1 ( 306/ 844- 383/1026) -63.4 E48-E40 S29-S18
      6 1327  892 (S08W20 L=243)   280    134        0 (1311/ 826-1338/ 951) -56.3 W18-W22 S11-S04
      7  156  825 (S32E69 L=153)  1132    141        1 ( 129/ 790- 202/ 877) -35.5 E74-E61 S34-S29
      8  444  836 (S28E35 L=187)   262     64        0 ( 412/ 824- 481/ 856)  -9.5 E38-E32 S29-S27
      9  467  479 (S53E37 L=186)   441     58        0 ( 445/ 458- 496/ 510) -31.0 E39-E33 S54-S50
     10 1303  299 (S43W46 L=269)   482     61        0 (1264/ 277-1330/ 321)  36.9 W42-W49 S44-S41
  

    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