Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  852 1540 (N36E15 L=225)   268     69        0 ( 820/1522- 892/1557)  18.4 E17-E11 N34-N37
      2  338 1505 (N31E62 L=178)   582     70        2 ( 318/1475- 366/1537) -15.3 E67-E56 N28-N32
      3 1609 1354 (N25W42 L=282)   326    114        0 (1560/1334-1664/1374) -12.5 W38-W47 N23-N25
      4  401 1254 (N15E45 L=195)   416    127        0 ( 384/1185- 428/1297) -66.8 E46-E42 N10-N17
      5  317  629 (S25E55 L=185)   559    122        0 ( 305/ 568- 343/ 677)  48.8 E58-E50 S29-S21
      6  675  533 (S30E24 L=216)   253     47        0 ( 647/ 523- 709/ 548)  11.3 E26-E21 S30-S28
      7 1474  397 (S38W41 L=281)   598    138        1 (1402/ 372-1529/ 430)  26.6 W34-W45 S40-S35
      8  987  365 (S42E01 L=239)   776    150        3 ( 937/ 327-1025/ 397)  56.3 E04-W02 S45-S39
      9  699  359 (S44E26 L=214)   319     54        0 ( 675/ 342- 730/ 380) -36.9 E28-E23 S45-S42
     10  802  328 (S46E18 L=222)   361     46        2 ( 777/ 302- 823/ 350) -76.0 E20-E15 S48-S44
  

    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