Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  463 1597 (N13E53 L=274)   434    117        2 ( 432/1565- 504/1645)  80.5 E55-E51 N10-N17
      2  310 1557 (N09E67 L=261)  2076    281        5 ( 257/1519- 361/1596)   4.1 E76-E61 N05-N11
      3  698 1473 (N10E31 L=297)   806    257        3 ( 629/1375- 756/1568) -32.7 E39-E24 N05-N15
      4  139 1166 (S17E68 L=260)  1113    278        1 (  85/1131- 182/1193)  14.0 E75-E62 S18-S15
      5  628 1161 (S09E26 L=301)   317     87        2 ( 607/1131- 651/1193) -36.9 E28-E24 S10-S06
      6  357 1075 (S19E44 L=284)   912    312        0 ( 256/1015- 454/1142)  -9.5 E53-E35 S21-S16
      7  189 1018 (S25E58 L=269)   754    189        0 ( 148/ 949- 235/1085)  60.9 E61-E54 S29-S20
      8  582  974 (S20E25 L=302)   625    166        0 ( 506/ 952- 652/ 998)  21.8 E30-E20 S23-S18
      9  440  773 (S36E33 L=295)   329     96        0 ( 398/ 748- 489/ 795)   0.0 E36-E28 S36-S34
     10 1385  406 (S31W44 L= 12)  2828    623        3 (1151/ 352-1626/ 463)  18.0 W26-W63 S35-S22
  

    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