Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  957 1709 (N53W07 L=188)   406     75        1 ( 926/1687- 989/1725)  15.9 W03-W10 N51-N54
      2 1472 1503 (N29W41 L=222)   275     70        0 (1443/1481-1503/1518) -53.1 W38-W42 N26-N30
      3 1382 1473 (N28W32 L=213)   324     67        0 (1363/1436-1403/1509)  78.7 W30-W34 N25-N30
      4  361 1415 (N37E51 L=130)   766     63        3 ( 327/1383- 397/1443) -16.7 E55-E45 N35-N39
      5  644 1322 (N29E23 L=158)   290     63        0 ( 612/1300- 669/1346)  21.8 E25-E20 N27-N30
      6  527 1298 (N28E32 L=149)   273     41        1 ( 504/1270- 552/1320)  56.3 E33-E29 N26-N29
      7  319 1136 (N19E50 L=131)   658    161        1 ( 274/1068- 354/1180) -41.2 E54-E46 N14-N22
      8  283 1065 (N15E53 L=128)   499    166        1 ( 258/ 992- 315/1146) -68.2 E56-E49 N10-N20
      9  171  926 (N07E68 L=113)  1158     96        2 ( 143/ 884- 200/ 973)  33.7 E73-E63 N03-N09
     10 1451  360 (S47W27 L=208)  1380    207        3 (1348/ 330-1530/ 403)   0.0 W16-W34 S49-S43
     11  511  302 (S38E62 L=119)  1034     79        1 ( 482/ 270- 541/ 334)  14.6 E76-E53 S41-S35
  

    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