Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  366 1555 (N56E54 L=101)   664     58        0 ( 336/1530- 395/1583)  26.6 E59-E49 N53-N57
      2 1790 1378 (N07W65 L=220)   699    101        1 (1761/1332-1817/1415)  23.2 W59-W69 N05-N09
      3  529 1376 (N39E28 L=128)   284     44        1 ( 502/1358- 553/1398)  11.3 E30-E25 N37-N40
      4  946 1353 (N28W03 L=158)   707    198        3 ( 862/1315-1050/1385)  -4.4 E03-W10 N26-N30
      5 1797 1039 (S12W52 L=208)   689     87        2 (1766/ 998-1823/1075)  63.4 W49-W55 S14-S09
      6  814  632 (S11E21 L=134)   259    108        0 ( 775/ 600- 864/ 667)  15.9 E24-E17 S12-S10
      7  354  514 (S12E62 L= 94)   577     49        0 ( 331/ 486- 380/ 543) -90.0 E63-E59 S13-S09
      8  520  413 (S21E52 L=104)  1660    245        2 ( 462/ 312- 575/ 516)  29.2 E68-E41 S27-S13
      9 1067  452 (S29E11 L=144)   239    139        0 (1003/ 440-1133/ 467)  -6.3 E15-E06 S30-S27
     10  634  376 (S25E44 L=111)   675    135        2 ( 578/ 346- 681/ 422)   9.5 E50-E38 S28-S22
     11  746  319 (S31E40 L=116)   336     57        0 ( 716/ 297- 771/ 344)  18.4 E42-E36 S32-S29
  

    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