Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1628 1468 (N28W51 L=301)  1670    225        4 (1562/1399-1708/1561) -71.6 W46-W57 N22-N34
      2  323 1402 (N30E56 L=194)  1295    264        1 ( 259/1329- 420/1445)  14.0 E65-E45 N25-N33
      3  828 1392 (N29E12 L=239)   275     87        1 ( 780/1370- 866/1406)   0.0 E15-E08 N27-N30
      4  337 1326 (N26E52 L=199)   745     90        2 ( 293/1309- 393/1347)   5.2 E56-E45 N24-N26
      5  696 1227 (N19E21 L=230)   298    106        1 ( 659/1188- 731/1282)  71.6 E23-E18 N16-N22
      6 1693  902 (S08W45 L=296)   533    115        2 (1664/ 857-1721/ 963)  63.4 W42-W47 S10-S04
      7 1749  857 (S11W51 L=301)   591     63        3 (1720/ 825-1779/ 881)   0.0 W48-W54 S13-S09
      8  616  550 (S24E32 L=219)   864    195        1 ( 589/ 462- 650/ 632)  80.5 E33-E29 S30-S18
      9  414  551 (S24E50 L=201)   459     90        1 ( 376/ 524- 443/ 588)  18.4 E54-E45 S25-S21
  

    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