Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  361 1612 (N29E57 L=322)   795    138        1 ( 341/1562- 376/1662) -66.8 E59-E54 N25-N32
      2 1508 1533 (N32W30 L= 48)   454     42        2 (1481/1509-1537/1557) -14.0 W27-W32 N29-N33
      3  346 1487 (N20E53 L=326)   433     82        0 ( 312/1468- 380/1512)  31.0 E55-E49 N18-N21
      4 1246 1363 (N17W11 L= 29)   465    137        1 (1215/1293-1272/1435) -71.6 W08-W12 N12-N21
      5  634 1342 (N12E27 L=352)   819    227        1 ( 558/1259- 726/1394) -24.8 E32-E19 N06-N14
      6  241 1251 (N05E57 L=321)  1153    188        3 ( 203/1174- 290/1302)  40.6 E60-E52 N00-N07
      7 1483 1207 (N10W27 L= 45)   308     75        0 (1445/1179-1518/1228) -21.8 W24-W29 N08-N11
      8 1778  991 (N01W51 L= 70)   726    166        0 (1746/ 913-1798/1077) -53.1 W47-W54 S03-N05
      9  296  748 (S25E50 L=328)  1949    484        1 ( 228/ 592- 397/ 916)  46.5 E60-E40 S35-S15
     10  476  441 (S46E42 L=336)   368     84        1 ( 457/ 410- 493/ 473)  45.0 E44-E39 S48-S43
  

    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