Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1505 1510 (N18W45 L= 53)  1041    241        1 (1410/1484-1613/1530) -26.6 W38-W52 N14-N21
      2 1804 1405 (N03W66 L= 75)   883    222        2 (1765/1370-1855/1438)  -3.4 W60-W77 N01-N04
      3 1586 1166 (S02W37 L= 45)   599    134        0 (1558/1105-1621/1237)  66.8 W34-W39 S06-N01
      4 1747  902 (S22W42 L= 51)  2686    536        6 (1664/ 672-1796/1161)  56.8 W31-W50 S34-S08
      5  158  969 (N23E63 L=306)   655    123        2 ( 128/ 947- 187/ 988)  -8.1 E66-E59 N21-N23
      6  215  870 (N16E59 L=309)   408     81        1 ( 188/ 845- 242/ 902) -53.1 E61-E56 N14-N18
      7 1521  727 (S25W21 L= 30)   259     60        1 (1510/ 696-1535/ 766)  78.7 W20-W22 S27-S22
      8  191  715 (N07E69 L=299)  1100    185        2 ( 148/ 680- 237/ 749) -21.8 E74-E63 N04-N09
      9  491  541 (S08E47 L=321)   283     56        0 ( 477/ 516- 510/ 578)  78.7 E48-E45 S10-S05
     10  378  482 (S10E60 L=308)   573    119        2 ( 348/ 444- 402/ 530)  71.6 E62-E58 S12-S06
  

    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