Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  918 1780 (N60W13 L=217)   745     89        3 ( 887/1742- 960/1806)  18.4 W09-W18 N57-N63
      2  289 1329 (N36E55 L=148)   658     77        1 ( 243/1301- 319/1351)   0.0 E62-E51 N34-N37
      3  850 1208 (N21E08 L=196)   794    330        1 ( 718/1177- 985/1249) -14.0 E15-W01 N17-N24
      4 1746 1034 (S07W49 L=252)   377     83        1 (1723/1004-1766/1077) -68.2 W46-W50 S09-S04
      5  160 1024 (N17E69 L=135)   933     90        2 ( 124/ 999- 207/1053)  -8.1 E75-E61 N15-N19
      6  992  794 (S06E06 L=198)   267     54        1 ( 962/ 772-1019/ 813) -45.0 E07-E04 S08-S04
      7 1328  750 (S15W14 L=217)   318    113        0 (1292/ 703-1366/ 794)  33.7 W10-W16 S17-S13
      8  259  696 (S03E62 L=141)   960    160        2 ( 214/ 641- 303/ 728) -39.8 E68-E58 S06-S01
      9  275  616 (S09E64 L=139)   628    117        0 ( 251/ 570- 299/ 677) -90.0 E66-E62 S11-S04
     10  486  640 (S09E42 L=161)   361     53        1 ( 469/ 608- 512/ 668)  45.0 E44-E40 S10-S06
     11 1791  619 (S37W56 L=260)  1288    132        1 (1730/ 592-1852/ 654)  -2.9 W47-W67 S38-S34
     12  972  621 (S17E10 L=193)   211     50        1 ( 950/ 600- 995/ 644)  26.6 E12-E08 S18-S15
     13 1222  501 (S30W03 L=207)  1297    287        1 (1143/ 403-1327/ 608)  33.7 E04-W11 S35-S24
     14 1567  548 (S35W30 L=234)   627    101        1 (1513/ 522-1609/ 585)   0.0 W25-W34 S36-S32
     15  681  537 (S18E31 L=172)   588    167        0 ( 617/ 513- 743/ 576) -20.6 E35-E26 S19-S14
     16  892  523 (S22E18 L=186)   291     71        1 ( 872/ 495- 927/ 547)  21.8 E19-E14 S23-S20
     17  536  418 (S24E48 L=155)   347     62        0 ( 522/ 384- 557/ 447)  33.7 E51-E45 S26-S22
     18  694  325 (S33E41 L=162)  2987    318        6 ( 634/ 180- 784/ 445)  26.6 E70-E28 S47-S24
  

    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