Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1405 1783 (N38W63 L=152)  1394    177        1 (1338/1727-1477/1818)  -5.0 W50-W73 N35-N39
      2 1365 1663 (N33W46 L=136)   964    146        1 (1311/1602-1410/1719)  73.3 W43-W49 N27-N37
      3 1507 1530 (N21W47 L=137)   311     53        0 (1477/1511-1538/1548) -63.4 W45-W49 N18-N22
      4 1866  958 (S22W58 L=147)   880    105        2 (1814/ 930-1902/ 998)   0.0 W50-W63 S23-S20
      5 1593  879 (S17W30 L=120)   840    197        0 (1569/ 783-1607/ 977)  63.4 W26-W33 S22-S11
      6  203  854 (N15E64 L= 26)  1930    192        4 ( 137/ 771- 250/ 925) -45.0 E73-E59 N09-N19
      7  288  814 (N11E55 L= 35)   400     50        2 ( 260/ 789- 316/ 844) -76.0 E57-E52 N09-N13
      8  269  713 (N05E60 L= 29)   447     41        0 ( 250/ 686- 286/ 740)  45.0 E62-E58 N03-N06
      9  548  608 (S06E40 L= 49)   293     50        0 ( 516/ 590- 571/ 630)  11.3 E42-E37 S06-S04
     10  926  410 (S27E24 L= 65)   532     90        0 ( 888/ 375- 961/ 449)  20.6 E28-E20 S28-S24
     11 1391  328 (S49W04 L= 94)   462     70        0 (1353/ 306-1431/ 354)   0.0 W00-W08 S51-S47
     12 1277  250 (S52E09 L= 80)   388     69        0 (1240/ 230-1312/ 266)   0.0 E13-E05 S53-S50
     13  769  240 (S35E46 L= 43)   464     78        0 ( 727/ 216- 802/ 265)  10.3 E52-E41 S36-S33
  

    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