Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1530 1611 (N33W42 L=166)   941    132        1 (1480/1558-1577/1643) -50.2 W37-W45 N29-N35
      2  713 1558 (N32E20 L=105)   882    252        0 ( 630/1462- 780/1629) -40.2 E27-E14 N25-N37
      3  988 1538 (N29E00 L=124)   901    269        4 ( 887/1482-1126/1583) -18.4 E07-W08 N25-N32
      4  863 1334 (N17E09 L=116)   367    192        0 ( 798/1260- 908/1381) -53.1 E12-E06 N11-N19
      5 1272 1314 (N14W16 L=141)   366    108        4 (1235/1268-1316/1351) -56.3 W14-W18 N10-N16
      6  850 1219 (N10E10 L=115)   333    116        1 ( 831/1160- 871/1275)  81.9 E10-E08 N06-N13
      7 1800 1077 (S01W53 L=178)  2480    540        6 (1692/ 954-1896/1158) -23.6 W43-W64 S08-N03
      8  392  970 (S03E41 L= 84)   366     54        3 ( 369/ 936- 413/1003) -63.4 E42-E39 S04-N00
      9 1343  922 (S10W19 L=144)   247    132        0 (1281/ 882-1391/ 951) -35.5 W15-W22 S12-S07
     10  453  863 (S09E37 L= 87)  1493    279       10 ( 366/ 811- 581/ 937) -21.8 E43-E28 S12-S04
     11 1646  851 (S15W41 L=165)  2532    417        5 (1487/ 744-1766/ 928)  18.4 W28-W50 S21-S10
     12  235  771 (S14E58 L= 66)  1600    346        3 ( 133/ 633- 305/ 892) -64.8 E69-E53 S22-S05
     13 1723  685 (S25W50 L=175)   697    198        1 (1644/ 611-1784/ 754)  42.0 W44-W55 S29-S20
     14  476  680 (S21E39 L= 86)   341     75        2 ( 460/ 645- 493/ 724)  78.7 E40-E36 S22-S17
     15 1012  642 (S26E02 L=122)   355    120        0 ( 978/ 605-1033/ 698)  51.3 E04-W00 S28-S22
     16  241  607 (S23E65 L= 59)   845     79        1 ( 199/ 588- 283/ 631)  -4.8 E72-E60 S24-S21
     17  661  455 (S36E31 L= 94)   315     64        1 ( 645/ 422- 679/ 491)  78.7 E32-E28 S38-S33
     18  592  380 (S42E40 L= 84)  1016    155        3 ( 530/ 324- 648/ 437)  29.7 E49-E33 S45-S37
     19  921  222 (S58E17 L=108)  1543    309        3 ( 771/ 182-1054/ 269) -12.1 E30-E02 S62-S51
  

    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