Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  982 1681 (N49E00 L=169)   370     55        1 ( 960/1647-1008/1707)  63.4 E02-W02 N46-N51
      2 1170 1656 (N46W16 L=185)   387     98        0 (1127/1639-1215/1675) -23.2 W12-W19 N44-N47
      3  377 1570 (N41E60 L=109)   439     77        0 ( 345/1556- 425/1581)   4.8 E65-E53 N39-N41
      4 1353 1499 (N33W27 L=196)   266     50        0 (1327/1474-1380/1520) -33.7 W24-W28 N31-N34
      5  690 1287 (N21E21 L=148)   659    191        2 ( 623/1237- 738/1326) -29.7 E26-E17 N18-N24
      6 1398 1266 (N17W26 L=195)   454    137        0 (1347/1219-1441/1321) -49.4 W22-W28 N14-N21
      7 1205  945 (S01W11 L=180)   577    251        0 (1105/ 877-1297/ 991) -26.6 W04-W16 S05-N01
      8 1711  798 (S14W48 L=217)   847    223        0 (1649/ 737-1772/ 866)  35.0 W42-W53 S17-S09
      9  860  813 (S08E11 L=158)   235    118        0 ( 841/ 776- 900/ 864) -71.6 E11-E08 S10-S04
     10  371  723 (S13E47 L=122)  1578    382        5 ( 333/ 555- 420/ 861)  61.2 E53-E41 S24-S04
     11  144  810 (S09E75 L= 94)  1131     71        2 ( 125/ 767- 168/ 840)  32.0 E81-E69 S11-S06
     12  947  714 (S15E06 L=163)  1047    402        5 ( 785/ 652-1077/ 785)   6.3 E16-W02 S18-S10
     13 1718  647 (S24W52 L=221)  2158    442        2 (1597/ 538-1847/ 771)  32.5 W42-W64 S30-S16
     14 1287  600 (S24W16 L=185)  1371    300        3 (1192/ 499-1393/ 694)  42.7 W09-W23 S30-S18
     15 1000  607 (S22E03 L=166)   850    187        4 ( 935/ 531-1053/ 659)  49.4 E07-W00 S27-S18
     16  507  314 (S44E56 L=113)   572     56        2 ( 486/ 297- 533/ 340) -71.6 E60-E52 S45-S41
     17 1109  237 (S53W04 L=173)   539     58        1 (1091/ 207-1131/ 269) -66.8 W02-W06 S57-S49
     18 1218  212 (S58W17 L=186)   676     84        1 (1190/ 183-1246/ 239)  45.0 W12-W20 S61-S54
     19 1295  189 (S62W30 L=199)   980     68        2 (1271/ 151-1313/ 232) -35.8 W24-W42 S68-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