Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1120 1809 (N46E23 L=329)   387     83        0 (1086/1779-1146/1834) -14.0 E27-E19 N43-N47
      2  972 1685 (N31E24 L=328)   358    103        0 ( 921/1654-1023/1701)   8.1 E27-E20 N29-N32
      3 1485 1620 (N44W13 L=  5)   823    114        3 (1444/1581-1525/1649)  -7.1 W08-W16 N41-N47
      4 1556 1596 (N46W19 L= 12)   503    105        1 (1537/1551-1584/1628)  81.9 W17-W22 N41-N49
      5 1720 1375 (N36W36 L= 29)  3668    479        7 (1593/1197-1797/1534) -40.1 W23-W44 N25-N44
      6  407 1480 (N04E53 L=299)   508    128        0 ( 366/1442- 432/1505) -18.4 E57-E49 N01-N05
      7 1065 1250 (N07E04 L=348)   265     84        0 (1048/1210-1087/1294) -59.0 E05-E01 N04-N09
      8 1562 1159 (N16W28 L= 20)   526    133        0 (1514/1090-1591/1227) -35.5 W22-W30 N13-N18
      9 1813 1089 (N21W49 L= 41)   363     64        0 (1785/1066-1837/1117)  -9.5 W45-W51 N20-N22
     10 1480 1005 (N05W26 L= 18)   753    157        2 (1422/ 955-1553/1068) -11.3 W21-W31 N03-N07
     11 1875  977 (N17W58 L= 51)   684     91        0 (1858/ 937-1894/1013) -53.1 W56-W61 N15-N19
     12 1728  932 (N09W46 L= 38)  1224    255        1 (1638/ 901-1810/ 954)  23.2 W39-W53 N05-N12
     13 1735  750 (N00W53 L= 45)   908    164        1 (1668/ 692-1786/ 823)  60.9 W49-W55 S04-N05
     14 1241  771 (S14W19 L= 12)   307    130        0 (1178/ 757-1309/ 790)  20.6 W15-W23 S16-S12
     15  744  739 (S30E09 L=344)   268     91        0 ( 709/ 711- 776/ 777)  59.0 E10-E07 S32-S27
     16 1146  630 (S25W19 L= 11)  1115    294        3 (1018/ 526-1271/ 708)  -5.4 W09-W30 S27-S21
     17 1541  587 (S15W45 L= 37)   509    128        0 (1472/ 537-1597/ 629)  66.0 W42-W47 S19-S10
  

    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