Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  899 1768 (N48E27 L=243)   451     46        0 ( 872/1747- 929/1788)  45.0 E30-E24 N45-N49
      2  563 1535 (N25E41 L=229)   776    177        1 ( 512/1452- 622/1623)  90.0 E43-E39 N18-N31
      3 1009 1437 (N25E07 L=263)   243     70        2 ( 990/1401-1023/1475)  90.0 E08-E06 N22-N27
      4  413 1354 (N11E47 L=224)   902    200        4 ( 379/1281- 446/1434) -41.6 E52-E42 N07-N15
      5 1145 1268 (N17W04 L=274)   421    177        1 (1066/1213-1204/1303)   8.1 W00-W07 N13-N18
      6  961 1012 (S02E04 L=267)   921    268        1 ( 903/ 898-1057/1132) -63.4 E07-W03 S07-N06
      7  729  924 (S10E17 L=254)   317    149        0 ( 675/ 880- 797/ 954)  40.6 E19-E12 S13-S07
      8  813  854 (S13E10 L=260)   196     38        1 ( 788/ 837- 839/ 877)   0.0 E12-E08 S14-S11
      9 1090  821 (S11W07 L=278)   304     92        1 (1058/ 789-1119/ 874)  80.5 W05-W08 S13-S07
     10  891  574 (S30E02 L=269)   243     72        0 ( 863/ 551- 929/ 590)  33.7 E03-W00 S31-S28
  

    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