Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1276 1806 (N58E08 L=218)   901     93        1 (1256/1755-1298/1848) -72.9 E11-E04 N51-N64
      2  994 1649 (N34E20 L=205)   729    237        2 ( 949/1554-1052/1739) -32.0 E28-E12 N28-N39
      3 1525 1551 (N44W22 L=248)   749     99        1 (1484/1505-1561/1590) -33.7 W18-W26 N40-N46
      4 1119 1485 (N26E06 L=220)   220     57        0 (1098/1470-1148/1508)   0.0 E07-E04 N25-N27
      5 1525 1340 (N28W24 L=250)   410    125        0 (1508/1285-1549/1398) -56.3 W22-W26 N25-N31
      6 1833 1170 (N27W56 L=281)  2535    411        5 (1723/ 987-1935/1310) -22.1 W42-W75 N18-N33
      7  202  772 (S38E63 L=162)  1019     87        1 ( 176/ 737- 235/ 808)  24.8 E70-E56 S40-S34
      8 1421  614 (S16W36 L=261)   310    121        0 (1357/ 595-1463/ 648)  50.2 W32-W37 S19-S13
      9  265  578 (S50E62 L=163)   696     83        0 ( 233/ 554- 296/ 614)  12.8 E75-E53 S52-S46
     10 1330  454 (S28W38 L=263)  2402    411        5 (1160/ 384-1512/ 517)  43.4 W28-W48 S37-S19
     11  415  434 (S59E44 L=182)   455     56        2 ( 393/ 408- 432/ 463)  33.7 E48-E38 S61-S55
     12  769  258 (S63W08 L=234)   366     85        0 ( 732/ 237- 806/ 287) -10.3 W02-W13 S63-S60
  

    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