Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  977 1634 (N32E12 L= 49)   327     81        0 ( 939/1608-1018/1655) -18.4 E14-E08 N30-N32
      2  293 1484 (N15E60 L=  1)  5312    388       11 ( 212/1330- 402/1615)  75.3 E75-E51 N05-N24
      3  380 1469 (N15E50 L= 10)  1677    263        8 ( 313/1371- 447/1559)  67.4 E53-E46 N08-N21
      4  185 1404 (N10E71 L=350)  1247     82        3 ( 156/1356- 211/1447)  90.0 E75-E66 N07-N13
      5 1721 1152 (N12W44 L=104)   883    232        2 (1684/1061-1758/1269)  74.1 W41-W47 N05-N19
      6  437 1231 (N01E39 L= 21)   301     76        0 ( 398/1207- 475/1246)   0.0 E42-E36 N00-N01
      7 1861 1121 (N12W58 L=118)  1620    203        2 (1815/1025-1891/1206) -60.9 W52-W61 N06-N17
      8  161 1072 (S10E64 L=357)  2095    173        5 ( 104/1000- 215/1149)  60.9 E72-E57 S14-S05
      9  493  818 (S23E33 L= 28)  2691    749        2 ( 288/ 615- 701/1138)  45.0 E56-E20 S37-S02
     10  248 1094 (S09E54 L=  7)   310     96        0 ( 217/1056- 291/1130)  33.7 E56-E49 S10-S06
     11  419  993 (S14E39 L= 22)   254     47        1 ( 395/ 967- 441/1015) -45.0 E40-E36 S14-S12
     12  771  923 (S15E14 L= 46)   342    139        0 ( 729/ 868- 799/ 984) -50.2 E17-E11 S17-S11
     13  240  886 (S21E55 L=  6)  2965    373        5 (  88/ 818- 417/ 979)   5.1 E83-E38 S24-S15
     14 1405  898 (S08W24 L= 85)   451    155        0 (1371/ 825-1431/ 970)  83.7 W23-W25 S12-S03
     15 1896  831 (S03W67 L=128)   533     75        2 (1874/ 789-1913/ 859) -38.7 W63-W70 S05-S01
     16 1340  708 (S20W24 L= 85)   650    222        0 (1265/ 670-1426/ 769)  -5.2 W18-W29 S22-S17
     17  999  439 (S42W07 L= 68)  2635    459        8 ( 776/ 399-1202/ 484)  10.6 E10-W23 S45-S37
  

    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