Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  456 1504 (N50E30 L= 10)   340     76        0 ( 417/1490- 499/1514) -16.7 E35-E25 N48-N51
      2 1029 1147 (N12W04 L= 44)   369    195        0 ( 999/1078-1075/1222)  54.5 W01-W08 N07-N15
      3 1378 1181 (N04W24 L= 64)   439    151        1 (1323/1135-1442/1221) -54.5 W20-W26 N00-N07
      4  744 1168 (N20E12 L= 28)   353    144        1 ( 674/1139- 808/1197) -36.9 E16-E08 N17-N23
      5  386 1084 (N24E39 L=  1)   500    115        1 ( 349/1039- 427/1116) -45.0 E42-E36 N20-N26
      6  121 1046 (N27E72 L=328)   778    119        1 (  97/1004- 150/1103) -28.6 E79-E66 N24-N30
      7 1578  859 (S19W28 L= 69)   321     46        2 (1552/ 833-1601/ 886)  14.0 W26-W30 S20-S17
      8  854  666 (S10E19 L= 21)   230     81        0 ( 814/ 652- 892/ 681) -21.8 E21-E16 S11-S09
      9 1245  429 (S36E04 L= 36)  1228    210        2 (1166/ 359-1311/ 512)  30.3 E10-W02 S39-S32
     10  754  320 (S29E41 L=360)  1740    214        4 ( 707/ 232- 803/ 410)  30.1 E52-E33 S34-S23
     11 1117  344 (S38E16 L= 25)   409     59        2 (1087/ 319-1149/ 368) -45.0 E18-E13 S39-S36
     12  994  283 (S39E28 L= 13)   677     95        1 ( 939/ 267-1048/ 305) -20.6 E32-E23 S41-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