Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1439 1815 (N36W70 L=252)  1032    114        1 (1405/1767-1462/1860)   4.4 W61-W87 N33-N37
      2 1220 1407 (N21W23 L=205)   836    245        0 (1110/1363-1325/1446)  -7.6 W15-W30 N19-N22
      3 1465 1203 (N03W30 L=212)   372     95        0 (1450/1153-1477/1252)  59.0 W28-W32 N00-N05
      4 1909 1069 (S19W65 L=248)  1951    245        5 (1867/ 970-1952/1175)  37.6 W57-W77 S23-S13
      5 1514  922 (S13W25 L=208)   368    159        0 (1491/ 853-1548/ 993)  49.4 W23-W29 S17-S10
      6 1468  592 (S32W15 L=197)  1122    316        1 (1394/ 479-1527/ 739) -87.3 W12-W19 S41-S20
      7  521  525 (S10E46 L=136)   617    143        1 ( 473/ 471- 574/ 577) -69.4 E47-E44 S14-S05
      8 1360  457 (S38W04 L=186)   548     48        2 (1329/ 422-1390/ 486)  78.7 W01-W06 S40-S35
      9  495  393 (S17E56 L=126)   663    118        0 ( 453/ 366- 523/ 433) -38.7 E60-E52 S19-S15
     10  575  397 (S19E49 L=133)   441     60        2 ( 559/ 365- 597/ 427)  23.2 E52-E45 S20-S17
     11  631  303 (S27E52 L=130)   470     47        1 ( 602/ 283- 654/ 322) -45.0 E54-E49 S28-S25
  

    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