Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  580 1133 (N23E25 L=319)   550    189        0 ( 495/1107- 663/1161) -28.6 E31-E19 N20-N26
      2 1574  999 (S08W32 L= 16)   426    127        0 (1536/ 948-1625/1026)  14.0 W28-W36 S10-S06
      3  885  868 (N02E12 L=332)   721    189        2 ( 864/ 763- 909/ 957)  68.2 E14-E09 S04-N06
      4 1414  844 (S13W18 L=  2)   798    209        0 (1349/ 752-1451/ 947)  45.0 W12-W22 S16-S07
      5 1711  828 (S23W40 L= 24)   823    231        0 (1660/ 716-1737/ 945)  58.4 W34-W43 S28-S15
      6 1272  724 (S16W07 L=351)   304     83        0 (1234/ 698-1319/ 745)   0.0 W04-W10 S17-S14
      7 1382  597 (S27W12 L=356)  1197    340        0 (1223/ 477-1509/ 733)  19.2 E01-W22 S31-S21
      8 1740  673 (S34W41 L= 25)   574     84        0 (1698/ 645-1775/ 702)   8.1 W37-W44 S35-S32
      9 1043  496 (S25E12 L=332)   262     63        1 (1017/ 468-1077/ 526)  21.8 E14-E09 S26-S23
     10  980  276 (S39E26 L=318)   471    126        1 ( 958/ 235-1007/ 320)  33.7 E30-E21 S42-S36
  

    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