Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  617 1637 (N29E44 L=152)  2209    328        3 ( 543/1517- 709/1773) -23.8 E66-E32 N21-N37
      2  834 1474 (N22E21 L=175)   767    174        1 ( 756/1440- 931/1496)   4.8 E26-E14 N20-N23
      3 1146 1401 (N24W00 L=197)   593    214        0 (1097/1349-1241/1468) -24.4 E04-W07 N22-N27
      4  306 1318 (N03E56 L=141)   469     60        1 ( 271/1302- 344/1335)  15.9 E59-E52 N02-N04
      5 1670 1262 (N26W40 L=237)   284     60        0 (1637/1241-1696/1277)   0.0 W37-W42 N24-N26
      6 1072 1177 (N09E00 L=197)   539    181        2 ( 987/1150-1143/1215)   6.3 E05-W04 N07-N10
      7  863 1131 (N02E11 L=185)   445     94        0 ( 804/1116- 920/1144)   7.1 E15-E07 N00-N03
      8 1319 1043 (N05W17 L=214)   246    152        0 (1295/ 980-1332/1118)  90.0 W16-W18 N01-N09
      9  716  821 (S19E15 L=181)   382     79        0 ( 696/ 781- 739/ 854)  78.7 E16-E13 S21-S16
     10  277  801 (S29E52 L=145)   643    140        2 ( 259/ 729- 292/ 849)  76.0 E54-E49 S33-S25
     11  780  702 (S25E09 L=187)   237     57        0 ( 750/ 686- 807/ 719) -14.0 E11-E07 S26-S24
     12 1514  343 (S33W59 L=256)   634     49        2 (1488/ 322-1544/ 369)  21.8 W55-W62 S34-S31
  

    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