Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  634 1452 (N28E23 L=134)  1034    284        2 ( 577/1359- 725/1565) -51.3 E28-E16 N21-N36
      2  884 1525 (N31E06 L=151)   305     99        0 ( 842/1497- 928/1558) -33.7 E08-E02 N28-N33
      3  367 1394 (N27E45 L=112)   675    142        1 ( 341/1338- 395/1467) -48.4 E49-E41 N22-N31
      4  251 1278 (N20E55 L=102)  1039    148        3 ( 214/1219- 297/1316)  56.3 E60-E50 N16-N22
      5 1408 1074 (S01W24 L=181)   506     95        2 (1368/1031-1442/1116) -63.4 W21-W25 S03-N01
      6 1877  975 (S10W63 L=219)   890    106        2 (1828/ 926-1912/1021) -59.0 W56-W68 S12-S07
      7 1796  978 (S09W53 L=210)   344     56        1 (1770/ 952-1825/ 999)  18.4 W50-W56 S10-S08
      8 1801  864 (S16W55 L=211)   562    100        2 (1767/ 827-1829/ 906)  33.7 W51-W57 S18-S13
      9 1860  865 (S16W62 L=218)   410     45        3 (1833/ 844-1878/ 895)  26.6 W58-W64 S17-S14
     10 1349  800 (S17W19 L=175)   205     98        0 (1313/ 774-1394/ 833)  31.0 W16-W21 S18-S15
     11  706  716 (S18E23 L=134)   277     82        1 ( 693/ 672- 735/ 750)  68.2 E23-E20 S20-S15
     12 1206  535 (S33W09 L=165)  2002    701        2 ( 995/ 375-1488/ 725)  26.6 E08-W28 S43-S22
     13  273  703 (S15E58 L= 99)   374     66        0 ( 254/ 677- 304/ 721)   9.5 E60-E54 S15-S13
     14 1732  636 (S30W52 L=209)   873    167        2 (1694/ 604-1786/ 672) -14.0 W47-W59 S32-S27
     15  464  621 (S21E43 L=114)   299     58        2 ( 441/ 596- 489/ 652)  45.0 E44-E40 S22-S19
     16  359  557 (S24E55 L=102)   561     89        1 ( 334/ 515- 390/ 598) -80.5 E56-E53 S27-S21
     17  751  551 (S28E23 L=134)   358     72        0 ( 716/ 530- 799/ 568)   0.0 E25-E19 S29-S27
     18  922  426 (S39E13 L=143)   900    205        0 ( 855/ 365- 994/ 501) -59.7 E16-E08 S44-S32
     19 1393  246 (S60W34 L=191)  1191    107        2 (1339/ 225-1456/ 280) -10.8 W25-W46 S62-S55
  

    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