Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  794 1650 (N52E06 L=248)   348     70        2 ( 775/1617- 821/1689)  45.0 E08-E02 N48-N54
      2 1778 1350 (N11W62 L=316)  1211    213        2 (1705/1294-1843/1444) -54.0 W58-W68 N05-N17
      3  224 1075 (N19E59 L=195)   948    174        1 ( 175/1050- 294/1114)   3.8 E66-E51 N17-N21
      4 1257  916 (S04W12 L=267)   309     94        0 (1207/ 904-1303/ 929) -18.4 W09-W15 S04-S02
      5 1364  527 (S31W15 L=269)  2996    572        2 (1139/ 389-1543/ 707)  24.3 E03-W28 S37-S22
      6  299  614 (S10E61 L=194)   522     43        1 ( 277/ 584- 321/ 634)  36.9 E63-E57 S11-S08
      7  998  587 (S20E09 L=246)   289     95        0 ( 962/ 546-1024/ 629) -71.6 E10-E07 S22-S16
      8  569  539 (S17E39 L=216)   406    113        0 ( 548/ 491- 598/ 574)  45.0 E41-E35 S20-S14
      9  796  388 (S31E29 L=226)   426    104        0 ( 779/ 345- 820/ 448)  76.0 E30-E26 S34-S26
     10  570  357 (S30E49 L=206)   822    138        0 ( 546/ 294- 615/ 414)  29.4 E57-E41 S35-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