Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1546 1290 (N16W38 L=199)   411    146        0 (1483/1269-1615/1311) -21.8 W33-W43 N14-N18
      2  716 1037 (N08E19 L=142)   391    103        1 ( 686/1004- 763/1077) -36.9 E21-E16 N05-N10
      3 1907 1019 (S05W70 L=232)  1108     80        1 (1868/1000-1948/1044) -12.5 W64-W82 S07-S03
      4 1811  829 (S15W57 L=219)   544    108        0 (1773/ 796-1845/ 853) -16.7 W52-W62 S18-S13
      5  714  704 (S13E22 L=139)   470    157        1 ( 651/ 678- 763/ 742)   7.1 E26-E18 S14-S10
      6 1352  606 (S24W19 L=180)  1942    502        5 (1194/ 437-1515/ 735)  38.0 W07-W30 S35-S16
      7 1750  634 (S28W55 L=216)   404     56        0 (1725/ 605-1777/ 659)  36.9 W52-W57 S29-S26
      8  599  340 (S39E44 L=117)   572     84        1 ( 581/ 295- 622/ 379)  37.9 E49-E39 S42-S35
      9 1072  268 (S49E04 L=157)   516     88        0 (1042/ 247-1112/ 297) -18.4 E06-W00 S51-S45
  

    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