Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  603 1401 (N35E19 L=177)   717    136        2 ( 552/1351- 663/1431) -14.0 E22-E14 N31-N37
      2  671 1348 (N30E15 L=181)   352     52        1 ( 652/1311- 695/1381)  53.1 E16-E12 N27-N31
      3 1500 1325 (N09W36 L=232)   326     71        1 (1468/1303-1532/1346) -11.3 W32-W37 N07-N10
      4 1669 1224 (S01W45 L=241)   320    116        2 (1619/1204-1701/1252)   0.0 W40-W48 S01-N00
      5  133 1075 (N24E64 L=132)   927    144        2 ( 111/1037- 161/1115) -59.0 E66-E60 N21-N26
      6 1206  965 (S04W09 L=205)   204     66        0 (1176/ 945-1245/ 981) -11.3 W06-W11 S05-S03
      7  803  733 (S08E19 L=177)  3944    895       10 ( 627/ 486- 951/ 905) -11.0 E31-E10 S21-N03
      8 1220  846 (S11W07 L=203)   306    153        0 (1154/ 825-1309/ 865) -11.3 W02-W12 S13-S10
      9  639  630 (S11E32 L=164)   360    110        0 ( 604/ 589- 673/ 673)  23.2 E35-E28 S12-S08
     10 1022  547 (S25E12 L=184)   515    150        0 ( 960/ 515-1103/ 590)   5.7 E15-E05 S25-S23
  

    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