Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  957 1642 (N48W07 L=326)   407     66        0 ( 926/1620- 993/1660)  18.4 W04-W10 N46-N49
      2 1052 1538 (N38W12 L=331)   720    177        1 (1014/1459-1101/1602) -80.5 W09-W13 N31-N43
      3  950 1442 (N33W02 L=321)   283     65        0 ( 931/1403- 962/1472) -90.0 W00-W03 N30-N35
      4  321 1121 (N20E49 L=270)  1293    326        1 ( 218/1037- 442/1207) -26.6 E62-E38 N14-N26
      5 1578  803 (S16W33 L=352)   552    213        0 (1495/ 746-1640/ 881)  26.6 W26-W38 S18-S12
      6  295  841 (N03E54 L=265)   359     37        1 ( 269/ 815- 322/ 861)  -9.5 E56-E50 N01-N04
      7  807  762 (S07E17 L=302)   785    180        4 ( 769/ 671- 848/ 834)  48.4 E21-E13 S11-S02
      8 1305  789 (S12W14 L=333)   330     75        1 (1285/ 749-1330/ 820)  45.0 W12-W15 S14-S10
      9 1681  743 (S23W42 L=  1)   726    150        1 (1660/ 681-1702/ 815)  90.0 W39-W44 S27-S18
     10 1312  551 (S28W12 L=331)   463    183        0 (1280/ 469-1351/ 630)  59.0 W08-W15 S33-S23
     11 1024  381 (S36E11 L=308)   294     36        0 (1006/ 353-1043/ 404)  63.4 E12-E09 S38-S33
  

    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