Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1706 1459 (N20W58 L=104)  1592    383        1 (1578/1329-1778/1598) -81.0 W52-W63 N12-N31
      2 1468 1476 (N26W39 L= 85)   394     80        0 (1429/1439-1500/1515) -71.6 W37-W40 N23-N29
      3  782 1384 (N32E11 L= 35)   239     76        0 ( 744/1370- 828/1395)  -8.1 E14-E07 N30-N32
      4  435 1294 (N30E38 L=  8)  1858    378        4 ( 300/1239- 581/1350)   4.8 E50-E26 N26-N32
      5  183 1233 (N27E63 L=343)   757    226        2 ( 145/1195- 223/1261) -26.6 E67-E58 N24-N28
      6  956 1132 (N14E03 L= 43)   556    145        2 ( 897/1093-1001/1169) -50.2 E06-W00 N10-N16
      7  358  639 (S10E53 L=353)  2871    585        5 ( 285/ 476- 413/ 811)  46.3 E69-E45 S22-N00
      8  583  729 (S07E32 L= 14)   515    114        2 ( 526/ 698- 648/ 753) -23.2 E36-E28 S09-S04
      9 1458  656 (S23W23 L= 70)   439     87        1 (1416/ 619-1500/ 695)  26.6 W20-W26 S25-S21
     10  514  661 (S10E39 L=  7)   328     78        0 ( 501/ 619- 536/ 689)  53.1 E40-E36 S12-S08
     11  474  367 (S31E54 L=352)  1733    357        1 ( 403/ 340- 552/ 412)  11.9 E63-E44 S33-S27
  

    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