Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  345 1103 (S16E45 L=207)  1229    248        1 ( 283/1006- 411/1202)  66.8 E49-E42 S22-S08
      2  244 1134 (S16E56 L=197)   495     83        1 ( 218/1100- 265/1171)  68.2 E57-E53 S18-S13
      3 1859 1057 (N21W56 L=309)   506    131        0 (1820/1035-1881/1106)  56.3 W52-W58 N18-N24
      4  426  983 (S21E37 L=216)   255     51        0 ( 401/ 968- 454/ 998)   0.0 E38-E34 S22-S20
      5  175  816 (S37E59 L=194)  1468    299        1 ( 128/ 780- 214/ 863)  31.0 E64-E54 S39-S33
      6  393  784 (S34E37 L=216)   818    164        0 ( 325/ 745- 457/ 826)  -8.7 E43-E30 S35-S32
      7  809  764 (S24E06 L=247)  1595    417        1 ( 597/ 688- 993/ 819)   6.3 E20-W07 S27-S21
      8 1359  672 (S15W30 L=282)  1829    502        4 (1155/ 630-1570/ 727)  17.8 W15-W43 S18-S06
      9  588  541 (S45E17 L=236)  1082    233        1 ( 505/ 511- 712/ 590)  14.7 E25-E06 S48-S40
     10 1243  382 (S35W36 L=289)   295     52        0 (1215/ 366-1270/ 402)  26.6 W34-W38 S36-S33
     11 1400  282 (S36W54 L=307)   504     59        1 (1365/ 269-1436/ 301)  53.1 W51-W55 S37-S33
  

    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