Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  905 1791 (N57W26 L=184)   361     47        0 ( 876/1776- 929/1809)  -9.5 W21-W28 N55-N58
      2 1277 1775 (N41W50 L=209)   746    127        1 (1216/1741-1339/1805)  -7.6 W42-W58 N39-N42
      3 1503 1593 (N23W51 L=209)  1022    138        1 (1433/1567-1575/1626) -12.1 W45-W59 N21-N25
      4  701 1479 (N41E08 L=150)   241     71        0 ( 663/1470- 738/1490)  -9.5 E10-E04 N40-N42
      5  574 1437 (N42E19 L=139)   573    131        1 ( 515/1404- 627/1471) -41.6 E24-E15 N38-N46
      6 1201 1398 (N21W21 L=179)   264    105        0 (1162/1378-1255/1428) -45.0 W19-W23 N18-N24
      7 1018 1248 (N18W06 L=164)   854    331        1 ( 884/1201-1176/1291) -15.5 E03-W15 N14-N20
      8 1525 1120 (S03W31 L=190)   284     48        0 (1500/1100-1560/1145) -45.0 W29-W33 S04-S01
      9 1624  981 (S13W35 L=193)   800    248        1 (1588/ 869-1659/1088)  85.9 W32-W36 S20-S06
     10  405 1022 (N21E39 L=119)   311     55        0 ( 374/ 998- 439/1047) -38.7 E41-E36 N18-N22
     11  788 1020 (N12E13 L=145)   363    157        1 ( 697/ 997- 833/1043) -14.0 E18-E10 N09-N14
     12 1251  992 (S02W12 L=170)   197     81        0 (1203/ 977-1276/1010)   0.0 W08-W13 S02-N00
     13  807  968 (N08E14 L=145)   342    140        1 ( 762/ 939- 872/1003)   8.1 E16-E09 N07-N09
     14  532  859 (N09E33 L=125)   228     64        0 ( 513/ 824- 551/ 895) -90.0 E34-E32 N06-N10
     15  296  827 (N12E53 L=105)   368    111        0 ( 262/ 809- 341/ 851)  -7.1 E56-E48 N10-N13
     16 1575  718 (S28W26 L=184)   407    139        1 (1535/ 663-1619/ 783)  32.0 W21-W30 S30-S24
     17 1219  579 (S25E01 L=157)   298     88        0 (1172/ 559-1263/ 606)  -9.5 E04-W02 S25-S24
     18  893  413 (S26E27 L=132)  1488    324        5 ( 861/ 297- 932/ 580)  50.2 E35-E20 S35-S15
     19  674  324 (S27E47 L=111)   371     62        0 ( 658/ 291- 693/ 355)  23.2 E50-E43 S29-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