Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  304 1562 (N09E68 L=260)  1594    244        6 ( 265/1525- 340/1596) -78.7 E76-E63 N06-N11
      2  708 1438 (N09E29 L=300)   450    158        0 ( 669/1385- 742/1509) -41.2 E34-E25 N05-N12
      3  157 1181 (S16E66 L=262)  1405    356        1 (  86/1134- 217/1234)  23.2 E74-E60 S18-S12
      4  612 1169 (S09E27 L=301)   282     57        2 ( 592/1140- 637/1197)  76.0 E29-E25 S10-S06
      5  346 1079 (S19E45 L=284)   861    287        1 ( 253/1021- 437/1144) -10.0 E54-E36 S21-S16
      6  163  999 (S27E61 L=268)  1198    361        1 (  90/ 937- 224/1074)  31.6 E68-E55 S30-S21
      7 1819  964 (N13W53 L= 22)   395     98        0 (1803/ 911-1829/1020) -56.3 W50-W55 N09-N15
      8  578  977 (S20E25 L=303)   700    189        0 ( 495/ 945- 668/1003)  13.0 E31-E18 S22-S18
      9  859  920 (S17E07 L=322)   309    136        0 ( 792/ 871- 902/ 972) -20.6 E11-E03 S19-S15
     10  436  774 (S36E33 L=296)   307    101        0 ( 392/ 751- 480/ 798)   0.0 E37-E29 S36-S34
     11 1334  760 (S13W25 L=354)   306     98        1 (1311/ 712-1364/ 793) -78.7 W23-W26 S15-S10
     12 1384  412 (S31W44 L= 12)  2858    605        2 (1148/ 356-1621/ 468)  18.4 W26-W62 S35-S22
  

    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