Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  881 1847 (N49E40 L=335)  1300    194        4 ( 849/1808- 917/1892) -30.3 E46-E34 N45-N53
      2  755 1384 (N13E24 L=351)  1381    337        2 ( 645/1279- 886/1490) -20.9 E34-E13 N09-N18
      3 1698 1257 (N27W41 L= 55)   703    171        1 (1638/1187-1740/1333) -37.9 W35-W44 N23-N30
      4 1377 1240 (N18W17 L= 32)   477    109        1 (1355/1188-1410/1302) -50.2 W14-W19 N15-N21
      5  780 1095 (S03E16 L=359)   618    134        2 ( 714/1065- 854/1129)  32.0 E19-E11 S06-N00
      6  977  900 (S10E00 L= 15)   496    131        0 ( 905/ 866-1043/ 932)   0.0 E04-W04 S11-S09
      7  292  760 (S34E48 L=327)  1983    564        2 ( 178/ 651- 398/ 854)  35.8 E62-E37 S43-S25
      8 1576  651 (S12W44 L= 59)   280     75        1 (1550/ 632-1604/ 680) -26.6 W42-W46 S12-S10
      9 1424  466 (S26W41 L= 56)   962    116        1 (1362/ 431-1483/ 500)  49.4 W38-W44 S29-S22
     10 1060  459 (S36W17 L= 32)   281     62        0 (1028/ 446-1097/ 475)  36.9 W15-W19 S37-S34
  

    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