Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  643 1638 (N30E34 L=193)   346     39        0 ( 618/1619- 675/1652)   0.0 E36-E31 N28-N30
      2  921 1463 (N20E11 L=216)   241     96        0 ( 874/1451- 982/1477)   8.1 E14-E07 N19-N20
      3  475 1249 (N03E37 L=190)   394     96        0 ( 457/1200- 499/1302) -80.5 E38-E35 N00-N06
      4 1422  814 (S13W27 L=254)  1395    336        2 (1285/ 719-1561/ 968)  54.2 W19-W34 S20-S02
      5  839  919 (S13E10 L=217)   246     61        1 ( 807/ 900- 877/ 939)  26.6 E12-E07 S14-S11
      6  654  774 (S24E21 L=205)   235     95        0 ( 617/ 752- 704/ 785)  18.4 E23-E17 S25-S23
      7  826  695 (S27E09 L=217)   414    134        2 ( 762/ 646- 882/ 732)  41.2 E13-E05 S31-S24
      8 1121  693 (S24W10 L=237)   215     74        1 (1090/ 671-1148/ 732)  63.4 W08-W12 S25-S21
      9 1592  654 (S20W43 L=270)   288     79        0 (1562/ 614-1611/ 696)  71.6 W41-W44 S23-S17
     10  659  344 (S55E25 L=201)   454     92        1 ( 631/ 321- 700/ 367)   6.3 E29-E20 S56-S52
  

    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