Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1651 1543 (N22W55 L=221)   721    137        0 (1601/1495-1698/1589) -66.8 W52-W58 N18-N25
      2  620 1509 (N36E21 L=145)   827    216        0 ( 537/1451- 730/1539)   3.8 E27-E12 N33-N38
      3 1192  724 (S20W06 L=172)   500    142        0 (1130/ 672-1241/ 780)  32.0 W01-W09 S22-S17
      4 1482  743 (S23W25 L=192)   226     54        0 (1470/ 713-1495/ 776) -76.0 W24-W26 S24-S20
      5  984  658 (S21E08 L=158)   232     59        0 ( 947/ 641-1024/ 678) -31.0 E10-E05 S22-S19
      6  281  617 (S13E61 L=105)  1587    143        3 ( 236/ 560- 314/ 676)  81.9 E65-E56 S16-S09
      7  912  630 (S22E13 L=153)   249     52        0 ( 899/ 592- 927/ 663)  53.1 E14-E11 S23-S19
      8  342  609 (S15E55 L=111)   573    100        0 ( 324/ 563- 363/ 662)  71.6 E57-E52 S17-S11
      9  679  598 (S20E29 L=137)   303     74        0 ( 654/ 562- 707/ 633)  31.0 E31-E26 S22-S18
     10  396  383 (S30E66 L=100)   847     93        0 ( 373/ 338- 420/ 422)  60.3 E70-E61 S33-S26
  

    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